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FOREWORD 


This study was performed by General Dynamics Convalr Division under NAS3- 
20092. Under a subcontract from November 1978 to May 1979, Science 
Applications, Inc. completed propellant manageir >nt systems comparisons 
and prepared the written draft of this report. This final report summarizes 
the technical effort from March 1976 to September 1979. Convalr program 
manager from March 1976 to August 1978 was M. H. Blatt. R. D. Bradshaw 
was program manager from Ai^st 1978 to September 1979. M. H. Blatt 
developed the refilling computer program, monitored the refilling tests, 
data reduction and modeling correlation, completed vapor Inflow testing and 
data correlation and performed the propellant management system comparison. 
F. Merino developed the vapor inflow analysis. L. E. Si den designed the 
experimental apparatus for both the liquid and v^or flow tests, R. Spencer 
and D. Uhlken were test conductors. J. A. Risberg provided technical 
support in running vapor inflow tests, reducing data from both liquid and 
vapor inflow tests and correlating the data v/ith analytical models. M. D. 
Walter developed capillary device and thermal subcooler weights. R. Makela 
determined propellant management system relative reliability. Chloe 
Bradshaw prepared the artwork and typed the final report. 

All data are presented with the International System of Units as the primary 
system and English Units as the secondary units. The English system was 
used for the basic calculations. 
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SUMMARY 


The study was concerned with three main areas (1) analytical and experimental 
investigation of refilling of start baskets with settled fluid, (2) iu<alytical and 
experimental investigation of the retention capability of wetted screens when 
subjected to vapor flow and, (3) comparison of alternative propellant manage- 
ment systems for the Centaur D-IS vehicle. 

In the refilling task, a computer program was written to include the effects of 
djmamic pressure, screen wicking, multiple screen barriers, standpipe 
(window) screens, variable vehicle mass for computing vehicle acceleration, 
calculation of tank outflow rate and vapor pullthrough height, options for wetting 
atid spilling. An experimental apparatus was designed , fabricated and tested 
to provide data for correlation with the analytical model. The test program 
was conducted in normal gravity using a scale model capillary device and 
ethanol as the test fluid. The test program provided data that was success- 
fully correlated with the analytical model. The model was then used to analyze 
Centaur D-IS LO 2 and LH 2 capillary devices for worst case mission conditions. 

Analysis of vapor inflow across initially wetted screens was performed using 
bench tests to provide a seml-emplrical model for the flow and pressure 
relationships applicable to the operation of a single screen window within a 
multiple screen capillary device. An extensive program of small scale bench 
tests was conducted using hexane, ethanol and Freon with traiisparent models 
simulating multiple screen/window configurations. The results obtained had 
lower repeatibility than was desired, creai^ing greater variability in the seml- 
empirlcal model predictions than is generally desirable for capillary device 
design. The worst case semi-empirical model predictions were used, however, 
to verify that the Centaur D-IS passive thermal conditioning system would 
have a high probability of operating successfully. A description is given of 
the recommended experimental prjgram to obtain repeatable data on vapor 
flow across wetted screen window configurations. 

Propellant management system alternatives were compared for the Centaur 
D-IS vehicle. Ten system concepts were comp.ared on the basis of payload 
weight penalty, hardware weight, relp.iive reliability, electrical power 
consumption and mission profile flexibility. Subsystem comparisons were 
made between settling rockets and start baskets; pressure feed and thermal 
subcoolers for boost pump NPSP (net positive sucticsi pressure); boost pumps, 
thenTial subcoolers and pressure feed for turbopump NPSP; cooled and uncooled 
ducts; and dumping any coolant required overboard or pumping It back into the 
tank. Comparisons indicate that the baseline Centaur D-IS system is a prime 
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propellant management system candidate based on payload weight. Depending upon 
the relative importance of hardware weight, reliability, power consumption and flight 
profile flexibility, promising subsystem alternatives worthy of consideration are 
thermally subcooled boost pumps, pressure fed turbopumps and capillary devices. 
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1 

INTRODUCTION 


The objectives of this study were to detennine the characteristics of liquid and vapor 
flow into a Centaur D-IS capillary device and to compare capillary devices with other 
propellant feed system alternatives. 

The use of capillary devices for replacing the baseline hydrogen peroxide settling 
system on Centaur was examined in detail in Reference 1-1. The analysis indicated 
that passively cooled start baskets (using layers of fine mes!i screen materials for 
wicking) coupled with thermal subcoolers (for replacing the main tank pressurization 
system) were promising for multiple burn missions. Subsequent work, reported in 
Reference 1-2, was conducted to study wicking configurations for passive cooling, 
and thermal subcoolers for replacing boost pumps. Experimental evaluation of the 
passive cooling configurations identified that closely spaced multiple layers of 
perforated plate and fine mesh screen material could satisfy all Centaur D-lS 
capillary device thermal conditioning requirements. Several important areas for 
investigation were identified as a result of the efforts undertaken in Reference 1-2. 
These areas were capillary device refilling with settled fluid and vapor flow across 
wetted screens. (Vapor must enter the capillary device to replace liquid evaporated 
by incident heat flux). A comparative evaluation of competing propellant management 
techniques for Centaur D-lS was identified as a useful method of determining the 
relative merits of promising subsystem alternatives. 

The baseline vehicle configuration for this study was the Centaur D-lS. The Centaur 
D-lS. as defined in Reference 1-3, represented a minimum change D-IT configuration 
(Centaur formerly used with Titan), modified to be compatible with the Space Shuttle 

interface, operations and safety require- 
ments. Approximately 95% of the existing 
Centaur components remain unchanged for 
the D-IS. Figure 1-1 summarizes signifi- 
cant modifications made to the existing D-IT 
to evolve the D-lS configuration. 

Mission profiles for the study were the 
planetary (1 burn), synchronous equatorial 
(2 bum) and low earth orbit (5 bum) flight 
profiles defined in References 1-1, 1-2, and 
1-3 and documented in Tables 1-1, 1-2, and 
1-3. 


MODIFICATIONS 
TO CENTAUR 
D1T 



EXPENDABLE 
(9B% EXISTING 
D IT 

COMPONENTS 

USED) 


Figure 1-1. Evolution of the Existing 
Centaur D-IT to a Centaur 
Compatible with Space Shuttle 
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Table 1-1. Planetary Mission Profile 





Event/ 
Time 
(min. ) 

Initial Maas, 
kKm (U>n,) 

Bum 

Time, 

sec 

Propellant 
Burned, 
kgm Ubm) 

Final Maas, 
kgm (Ihm) 

Initial 

Percent 

FuU* 

Initial 

Accel- 

eration, 

g 

Loading 

LO2 

11,554 (25,450) 






(T=0) 

LH2 

2,397 (5,279) 






MESl 

Vehicle 22,434 (49.413) 



8|888a9,576) 



(Tb67) 

LO2 

11.488 (25,304) 

441.4 

11.302 (24,894) 

186 (410) 1 

96 

' 0.61 


LH2 

2,344 (5,164) 


2,243 (4.941) 

101 (223) 

95 

j 


^AsBuznes 11, 946 kgg, (26,313 Ibm)* LO2 2,478 kgjQ (5,459 Ibm), LH2 for full tank. 


Main engine thrust 13,620 kgf (30,000 Ibf) 

Maximum ACS thrust * 10, 9 kgf (24 Ibf) 

Maximum ACS acceleration before last bum = 4.86 x 10*^ g*8 
Main engine flow rates, LO2 = 25.6 kg/sec (56.4 Ib/sec) 

LH2 5. 1 kg/sec (11. 2 Ib/sec) 


ISP B 443. 82 see 

Payload « 6,667 kg^ (14,465 Ib^i) 
Dry wei^ 2,015 kgm (4,439 Ibjg) 
Burnout acceleration « 1. 53 g* s 


Table 1-2. Synchronous Equatorial Mission Profile 


Event/ 
Time 
(min. ) 

Initial Mass, 

kgm (Ibjjj) 

Bum 

Time, 

sec 

Propellant 
Burned, 
kgm (“hn) 

Final Mass , 
kgm (l^m) 

Initial 

Percent 

Full* 

Initial 

Accel- 

eration, 

g 

Loading 

LO 2 

11,554 (25,450) 






(T = 0) 

LH2 

2,397 (5,279) 






MESl 

Vehicle 

21,541 (4f,447) 



12,100 (26,653) 


1 

(T = 67) 

LO 2 

n,488 (25.304) 

305.4 

7,854 (17,299) 

3,634 (8,005) 

96 

0.63 


LH 2 

2.344 (5,164) 


1,525 (3,360) 

819 (1,804) 

95 


MES2 

Vehicle 

12,162 (26,788) 



7,956 (17,525) 



(T = 384) 

LO 2 

3,621 (7,975) 

132.3 

3,403 (7,496) 

190 (419) 

30.3 

! 1.12 


LH 2 

782 (1,723) 


665 (1,465) 

117 (258) 

3L6 



^Assumes 11,946 kgm (26,313 lbm)i LO2 2,478 kgxn (5,459 lbm)f LHo for full tank. 


Main engine thrust » 13,620 kgf (30,000 Ibf) 

Maximum ACS thrust « 10. 9 kgf (24 Ibf) 

Max. ACS acceleration before last bum > 8. 96 x 10*^ g's 
Mixture ratio •= 5 . 0 

Main engine flow rates, LO2 - 25. 7 kg/sec (56. 65 Ib/sec) 
LH2 = 5.01 kg/seo (11.03 Ib/sec) 
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« 



ISP = 443.35 sec 

Payload = 5,538 kgj^ (12, 199 Ibm) 
Dry weight « 2,090 kg^^ (4,604 Ib^) 
Burnout acceleration = 1.71 g*s 


n , , ■‘in i I i>i«>rt«Illlll.lri. •‘t — ■ - — - 





Table 1-3. Low Earth Orbit Mission Profile 


Event/ 
Time 
(min. ) 

Initial Maas, 

kgua (Ibjij) 

Bum 

Time, 

sec 

Propellant 
Burned, 
kgm (Ibm) 

Final Mass 
kgm (Jb®) 

Initial 

Percent 

Full^ 

hiltlii 

Accel- 

eration, 

a 

Loading 

LO2 

xl.554 (25,450) 







0 

II 

LH2 

2.397 (5,279) 







MESl 

Vehicle 

19,090 (42,049) 


' 


16,363 (36,042) 



(T-67t 

LO2 

11,488 (25,304) 

88.6 

2,294 

(5,052) 

9,194 (20,252) 

96 

0.71 


LH2 

2,344 (5,164) 


434 

(955) 

1,911 (4,209) 

95 


MES2 

Vehicle 

16,264 (35,824) 




10,373 (22,849) 



(T = 118) 

LO2 

9,155 (20,165) 

191,32 

4,955 (10,915) 

4.200 (9.250) 

77 

0, 84 


LHg 

1,885 (4,153) 


935 

(2,060) 

950 (2,093) 

76 


MES3 

Vehicle 

10,246 (22,568) 




6.536 (14,397) 



(T = 408) 

LO2 

4,162 (9,167) 

120.51 

3.121 

(6,875) 

1,038 (2,286) 

35 

1. 33 


LH2 

913 (2,010) 


587 

(1.294) 

325 (716) 

37 


MES4 

Vehicle 

6.443 (14,192) 




5,861(12,910) 



(T = 459) 

LO2 

998 (2. 198) 

18,90 

489 

(1,078) 

509 (1,121) 

8.4 

2. 11 

1 

LH2 

295 (650) j 


93 

(204) 

207 (456) 

11.9 


MESS 

Vehicle 

5,765 (12,698) 




5.433 (11,967) 



(T = 533) 

LOj 

468 (1,031) 

10,8 

279 

(614) 

189 (417) 

3.9 

2.36 


LH2 

178 (393) 


53 

(116) 

126 (277) 

7.2 



♦Assumes 11,946 kgm (26,313 2,478 kg^j (5,459 Ibjjj), LH2 for full tank. 

Main engine thrust ■ 13,.620 kg£ (30, 000 Ibf) ISP » 443,8 sec 

Maximum ACS thrust » 10, 9 kgf (24 Ibf) Payload « 2,642 kgj^ (6,260 Ibu^) 

Max. ACS acceleration before 5th burn = 1.89x 10*5 g^g Dry weight « 2,225 kg^ (4,901 Ib^) 

Mixture ratio » 5.298 Burnout acceleration = 2.51 g's 

Main engine flow rates, LO2 = 25. 9 kg/sec (57. 05 Ib/seo) 

LH2 = 4,89 kg/sec (10,77 Ib/sec) 




CAPILLARY DEVICE REFILLING DURING ENGINE FIRING 


Partial control capillary acquisition devices (start baskets) operate by retaining liquid 
over the tank outlet between main engine bums. The start basket provides liquid flow 
to the main engines for engine restart and for a sufficient period of time to settle the 
propellants and refill the capillary device in preparation for the next bum. A 
typical start basket is shown schematically in Figure 2-1. 

Due to vehicle drag the initial position of 
propellant in the tank will be at the 
forward end of the tank. Engine thmst. 
which builds up during the start sequence 
to the steady state thrust level, will 
settle the propellant to the aft bulkhead. 

Refilling of the capillary device must be 
accomplished in the subsequent high-g 
period coincident with settling and main 
engine steady state operation. Refilling 
can be accomplished during this period 
because the retention capability of the 
screen is sized for much lower accelera- 
tion levels (e. g. , 1 x 10”3 g's) than 
occur during main engine firing (e. g. 

1 to 3 g's). For refilling to be successful, 
liquid must enter the capillary device in sufficient quantity to sustain engine firing 
without vapor pullthrough and refill the capillary in preparation for the next engine 
start. 

The process of refilling is a basic element in start basket operation. Analysis of 
refilling is required to assure that the start basket functions successfully during restart. 
An analysis was performed in Reference 2-1 for the LO 2 and LH 2 start baskets using 
conservative assumptions. Only hydrostatic pressure was assumed as the driving 
pressure. (No refilling augmentation due to liquid velocity impinging on the screen 
was used. ) Refilling was assumed to start after settling was complete. Screen wetting 
was assumed to exist during the entire refilling period. 

Refilling was computed based on pressure differences between the inside and the outside 
of the capillary device. Equations were formulated, as described in Reference 2-1, for 
capillary device geometry, and flow/pressure drop. Capillary device volumes were 
divided into shapes that permitted close form integration of the flow into each volume 
segment as a function of screen open area. 
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b the present contract the analysis of refilling was broadened to include the effects of 
dynamic pressure, screen wicking, multiple screen barriers, window (standpipe) 
sc-;.*eens that can be different mesh than the main screens, and time dependent liquid 
'X.!^tling (collection). The analysis also included the effects of vehicle mass on the 
vt^ hide acceleration as the propellant tanks are being drained. Outflow from the 
f ink was included either as an input or as a calculated value based on feed system 
pressures. Other analysis features are: calculation of vapor pullthrou^ hei^t in 
the basket based on tank outflow rate (and channel geometiy),an option for including 
a standpipe screen in the calculations, options for maintaining the standpipe in either 
a dry or wetted condition, an option for simulating liquid spilling from the capillary 
device at initiation of tank outflow (from the start basket), and an option for selecting 
t e type of multiple screen barrier to be used for the main screen in the screen wicking 
r ilculations. 

Because of the complexit A the analysis a closed form solution did not appear to be 
feasible. A finite difference solution was constructed using the basic screen flow 
equations and the continuity equation. The screen flow equations are the heart of the 
computer program developed for analyzing capillary device refilling. 

An experimental program was conducted in normal gravity using a scale model capillary 
de^dce and ethanol as the test fluid. The test program, while limited in scope due to the 
uniform gravity conditions during both retention and refilling, provided data that was 
used to successfully correlate the analytical model. The model was then used to predict 
Centaur D-IS start basket refilling. 

The followir . sections describe the refilling analysis computer program, the normal 
gravity non-cryogenic testing and the test data correlation. Documentation of the 
computer program is presented in Appendix A. 

2.1 REFILLING ANALYSIS 

An analysis of capillary device refilling with settled fluid was performed as a tool for 
design and an: I' sis of capillary acquisition device (start basket) sizing. A computer 
program was written to incorporate all significant analysis elements. Computer 
program documentation is presented in Appendix A. Input irequirements are fluid and 
screen properties, start basket and tank geometry vs hei^ above the bottom of the 
tai liquid collection vs time, screen impingement velocities, vehicle mass 
^ roperties, tank pressure and outlet geometry. 

A flow chart of the computer prognram is shown in Figure 2-2. The heart of the program 
is the flow equations governing vapor or liquid flow into or out of the start basket. The 
program elements determine the boundary conditions for calculating start basket fluid 
flo.vs uased on liquid level, system pressure and screen configuration and surface area. 
Continuity is satisfied by summing the flows into and out of the capillary device and 
adjusting the pressure difference between the inside and outside of the basket until the 
flows balance. This iterative solution is achieved for time step increments until either 
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the total bum time or the pullthrough height is reached. 

The pressure loss equations are separately delineated for each of six regions for fluid 
flow (without liquid Impingement). The six regions (for net refilling) are: an unwetted 
region where vapor flows out of the basket into vapor, a wetted region where vapor flows 
out of the basket into vapor, a region where the liquid level outside the basket covers 
the screen and vapor flows into the liquid, a region where no flow occurs across the 
screen barrier, a region where liquid enters the basket flowing into vapor and a region 
where liquid enters the basket flowing into liquid. Figure 2-3 depicts these regions. 

All six regions will not exist simultaneously. The delineation presented in Figure 2-3 
illustrates the entire range of possible conditions. 

VAPOR OITFLOW 



Figure 2-3. Schematic of Capillary Device Refilling Flow Regimes 


DPIN = Pu - Pa (2-1) 

DPS = 2. 72 a/Dgp, wetted screen retention pressure (2-2) 

where 


DPIN = pressure differential between vapor inside the basket and tanx ullage 

DPS = pressure difference that can be sustained by surface tension 

^ - liquid simface tension 

Dfip = screen bubble point 

For the equations presented in the following paragraph 
Pq = pressure in liquid outside the basket 

Pi = pressure in liquid inside the basket 

A&B are viscous and Inertial constants for vapor flow across the screen 

Y&Z are viscous and inertial constants for liquid flow across the screen 
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Liquid flow across the screen is computed from the following equation 
= YV + ZV2 

where V is the liquid velocity through the screen. 

Vapor flow across the screen is computed from 
APv - AV + BV^ 

where V is the vapor velocity across the screen. 


(2-3) 


(2-4) 


Substituting the pressure relationships for each region and solving Equations 2-3 and 
2-4 by using the quadratic formula yields the following e}q>ressions used in Subroutine 
FLOW. 

For Region I with unwetted screen and vapor flow into vi^: •. 


-A +V A^ + 4B (DPIN) 
2B 


For Region II with wetted screen and vapor flow into vapor 
-A +/a^ + 4B (DPIN-DPS) 


VVW = 


2B 


For Region in with liquid outside the basket and vapor flow Into liquid 


vvT _ - A + V A^ + 2B iPPlN-DPS> 
^ ~ 2B 


For Region IV no flow occurs. 

For Region V liquid flow into vapor occurs 

-Y +V y2 + 2Z (P^ - Pj - DPIN) 


VL = 


(2-5) 


(2-6) 


(2-7) 


2Z 


(2-8) 


For Region VI with liquid flow into liquid 

- Y +V Y^ + 4Z (P - Pi - DPIN) 

A * 


VT = 


2Z 


(2-9) 
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The outflow velocity from the start basket is computed lu subroutine OUTFLOW from 
the following equation 


V0UT = \/6 * 4DPF/(KC * CC * Pl> (2-10) 

where 

DPF is the pressure difference across the outlet (average during the time id;ep)* 
DPF is a function of DPIN 

KC is the contraction loss coefficient for the outlet 
CC is a correction factor for a rounded entrance 
p ^ is the liquid density 

The continuity equation evaluated using flow velocities and areas 

VV (AV) + VVW (AVW) + WL (AVL) + V0UT (A0UT) - VL (AL) 

- VT (AT) = 0 (2-11) 


An iterative solution is performed by varying DPIN until Equation 2-11 is satisfied. 

Screen impingement is included in the liquid refilling flow in Regions V and VI as 
applicable based on the dynamic pressure aiding refilling using the e^qpression 


V0 


-FKA + 


J 


(FKA)2 + 4 pV^ 

...kJ , 

2 (FBK + p K) 

L 


(FBK + p^ K) 
L 


( 2 - 12 ) 


for the velocity of the emergent jet incident on a screen 

F is a constant used to determine the pressure drop through the screen 
(empirically found to be equal to 2) 

K is a proportionality constant equal to the emergent Jet liquid area/impinging 
jet liquid area (empirically found to be equal to 4) 

V] is the velocity of the impinging jet 

A t B and P|^ are as previously defbsed. 

Equation 2-12 was obtained from Reference 2-2. 



Sample nms were made with the RL'FILL computer program to test program capability. 
Program runs were made in order to design an experimental refilling apparatus. 

Several of the program capabilities described in Figure 2-2 such as subroutine STAND, 
variable acceleration capability and spilling options were added based on model 
correlation requirements described In Section 2-4. 

2.2 REFILLING EXPERIMENT PREPARATION 

An experimental apparatus and normal gravity test procedure was developed with the 
following objectives: simulate actual conditions for refilling as much as possible 
considering the nonvariant ambient acceleration, simulate actual settling and collection, 
(with realistic screen impingement), run cases where both net "efilling and net draining 
of the start basket occurs, run with tank outflow at both zero and nonzero values, and 
simulate realistic standpipe and channel flow operation. 

2. 2. 1 TEST MODEL AND EXPERIMENTAI. APPARATUS. A complete set of design 
drawings are shown In Figures 2-4 through 2-6. Figure 2-4 Is a general arrangement 
of the refilling apparatus. Figure 2-5 details the multiple screen layer start basket. 
Figure 2-6 shows the screened standpipe assembly and manufacturing details. 

The basic configuration shown in Figure 2-5 consists of a start basket with multiple 
screen barriers for thermal conditioning and a screened window on top of the basket 
to allow vapor penetration during thermal conditioning and vapor removal during 
refilling. Fine mesh screened channels within the basket were designed to provide 
liquid outflow from within the start basket and will be maintained full of liquid during 
testing. The channel screen was 325 v 2300 mash. The start basket screen was 
multiple barriers of 200 x 600 screen. 'Fv/o of the basket walls were made of Lexan 
to allow backlighting, visual observation, and photographs to be employed. The top 
screen window on the basket (view M-M) was either 200 600 or 50 x 250 mesh 

screen. With 50 x 250 mesh screen, 23. 45 cm (4. 9 in) of liquid head could be 
retained in the basket when it is completely wetted, (similarly 12. 45 cm (4. 9 in) 
of vapor will be trapped In the completely wetted basket during refilling). With 
200 X 600 screen the basket could hold its own head of liquid or vapor. A single 
layer removable standpipe of 50 x 250 screen (View Q) was used for some of the 
runs. The start basket cuter envelope was a square prism, 16. 5 cm high x 45. 72 
cm wide (6.5 in x 18 in wide). Channels were 22 cm (8.5 in) long x 15.2 cm (6 in) 
wide and 0. 95 cm (3/8 In) thick. The outlet tube was 5. 08 cm (2 in) in diameter. 

A square tank was fabricated with front and back sides of Lexan sheet to serve as the 
outer tank. Within this tank a movable cover over the basket controlled the collection 
rate of liquid. This cover consisted of a square cross section box with an open frame- 
work on top. The cover has front and back sides constructed of Lexan. The cover 
keeps liquid betv.’een the tank wall and the cover until the cover is lifted. Lifting the 
cover allows liquid collection around the basket. Splashing of liquid against the screen 
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was controlled by deploying a fence for some of the runs to cause liquid spilling over 
the fence when the cover was removed. The cover was lifted at different rates to 
simulate different collection times. A piston/stop arrangement was used to control 
cover motion. This approach was selected because the collection rate is controllable 
and the geometry can be made to simulate the wall flow that would occur during an 
actual axlsymmetrlc reorientation case. 

A plumbing schematic of the capillary device refilling apparatus is shown in Figure 2-7. 

A schematic showing instrumentation is presented in Figure 2-8. 

2. 2. 2 TEST VARIABLES AND INSTRUMENTA TION. Ethanol was used as the test fluid 
because of its good wettability, low toxicity and low vapor pressure at room temperature. 
The properties of ethanol allowed reasonable start basket heights and flow rates to be 
used. 

Surface geometry was recorded photographically on 16 mm color movie film. A scale 
with 0 . 1 mm divisions was mounted on the start basket Lexan surface. Liquid level 
Inside and outside the basket was measured using General Dynamics Convair fabricated 
parallel plate capacitance probes. Variable reluctance type pressure transducers were 
used to measure the pressure difference across both the start basket and channel. 

Flow rate out of the basket was measured with a turbine flow meter in the outlet line. 
Temperatures were measured with chromel/constantan thermocouples. The tempera- 
tures measured were Tj, the temperature at the start basket AP transducer; T 2 , the 
temperature at the channel AP transducer; T 3 , the temperature at the top of the start 
basket outer screen; T 4 , the temperature In the outflow line; T5, the fluid temperature 
In the outer annulus; and T 0 , the supply tank temperature. 

Analog recordings were obtained for each run. Data channels recorded were outer 
chamber liquid level, start basket liquid level, start basket pressure drop, channel 
pressure drop, channel outflow rate, and the six selected temperatures. All quantitative 
data was obtained from the analog recordings. The motion picture runs were used for 
qualitative observations and for assistance in determining boundary conditions such as 
settling flow pattern and screen wetting during refill. 

Test variables were outflow rate, top screen (standpipe) configuration, initial start 
basket wetting, liquid collection rate and refilling geometry. 

After assembly of the start basket, bubble point testing was conducted preparatory to 
c^^mpleting test set up assembly. Leakage occurring at the curved scarfed Lexan 
spacers between the 200 x GOO screen layers was sealed with fast curing epoxy. Shown 
in Figures 2-9 and 2-10 is the assembled test apparatus. The start basket shown 
does not have a standpipe. Figure 2-10 shows the outlet region under the test model. 
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Figure 2-5. Inner Basket Assembly/Refilling Test Apparatus (Sheet 1 of 5) 





Figure 2-5. biner Basket Assembly/Refllling Test Apparatus (Sheet 2 of 5) 




Figure 2-5. Inoer Basket Assembly/Refilltng Test Apparatus (Sheet 3 of 5) 



Figure 2-5. bmer Basket Assembly/Refllling Test Apparatus (Sheet 4 of 6) 
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Figure 2-5. Buaer Basket Assembly/Refilling Test Apparatus (Sheet 5 of 5) 
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Filcure 2-6* Standpipe Assembly/RefillinfiC Test Apparatus (Sheet 1 of 2) 






Figure 2-7. Capillary Device Refilling Apparatus Plumbing Schematic 




Figure 2-8. Schematic Showing Test Apparatus hastrumentatlon 
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Figure 2-9. Refilling Test Model 
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2.3 REFILLING TESTS 




A series of runs were made with the test configuration shown in Figures 2-9 and 2-10 
in order to check out instrumentationt valvest and a]q>aratus controls. After checkout 
a series of twenty nine runs were made with data taken as described in Section 2. 2. 2. 
Results were obtained over a wide spectrum of iv^fllllng times ranging from no ref'Uing, 
with the 200 x 600 mesh screens, to rapid refilling with the 50 x 250 mesh standpipe 
and initially dry screens. 

The matrix of test runs is presented in Table 2-1 showing the screens tested, the 
outflow rate, initial wetting condition, liquid collection rate, start basket initial 
liquid level at start of collection and start of outflow. R<ms recorded on motion 
picture film are noted with an "M. " 

For the start basket top screen, three configurations were used; 50 x 250 mesh and 
200 X 600 mesh flat screens and a 50 x 250 mesh standpipe 3.81 cm (1.5 In) high and 
5. 08 cm (2 in) in diameter. 

The majority of runs were made with the start basket screen wet. Several runs were 
made with the basket initially empty and no outflow (Runs No. 10,11,19,23 and 24) 
with the screen therefore initially dry. Other runs were made with dry top screens. 

In these runs the basket was filled to a low level and the liquid level between the basket 
and the cover was initially equal to the liquid level in the basket. 

Zero time for each run was at the timing mark on the analog recording just preceding 
the lifting of the cover to allow liquid to collect over the basket. Timing marks were 
made at one second intervals. Liquid levels and flow rates are approximately linear 
between the time values shown in Table 2-1. 

A typical test sequence occurred as follows 

1. Dry top screen with GN 2 purge and drain residual liquid from test tank in 
preparaticm for start basket filling. 

2. Fill basket and inner annulus to oeslred level with ethanol. 

3. Drain inner annulus (between start basket and cover) as required. 

4. Fill outer annulus (cover is down) to desired level. 

5. Open outflow valve (if flow is desired) and drain start basket. 

6. When start basket is at desired level, lift cover (cover lift time is as desired) to 
collect liquid over the start basket. 
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Table 2-1. Refilling Test Data 
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Table 2-1. Refilling Test Data (Continued) 


itun 

No. 

Baakct 

Top 

Scraafi 

initial 

Top 

Scteen 

Wettlns 

Tint* 

aeo. 

Stan Baakec 
Internal 
Liquid Level, 
cm (In) 

Start Basket 
External 
Liquid Level, 
cm (in) 

FloerCbanael, 
OP, kS/m^ 
(in lifO) 

Start Basket • 
OP, kti/m^ 
((nH20) 

Flow Rate, 

em^/$€c 

(CPM) 

Start 1 

Basket 1 

Top j 

Temp. , 
K(F) j 

14 

50x250 

Dnr 

0 

3.8 

(l.S) 

1.5 

(0.0) 

1.17 

(4.7) 

-0.12 

(-0.5) 

1 

100.9 

(1-9) 

292 

(66) 




1 

3.8 

(1.5) 

10.4 

(4.1) 

1.47 

(8.9) 

-0.12 

(-0.5) 

i 100.9 

(U«| 

292 

(66) 



1 

2 

4.6 

(1.8) 

16.3 

(6.4) 

1.17 

(4.7) 

0 

(0) 

1 94.6 

(1.51 

292 

(66) 




4 

8.1 

(3.2) 

18.8 

(».«) 

0.6 

(3.4) 

-0.5 

(-2.0) 

94.6 

(1.5k 

292 

(66)1 




10 

9.4 

(3.7) 

23.4 

(0.3) 

0.37 

(1.S) 

-0.37 

(-1.8) 

107.2 

(1.1) 

292 

(66)1 

i 

) 



47 

13.5 

(5.3) 

21.6 

(8.5) 

0 

(0) 

-0.25 

(-1.0) 

100.9 

(1.6) 

292 

(66)1 

1 



315 

14.2 

(5.6) 

16.3 

(«.«) 

0.17 

(0.7) 

-0.12 

(-0.5) 

107.2 

(1.7) 

292 

(66)1 

1 



3S0 

i 14.0 

(5.5) 

15.2 

(6.0) 

0.17 

(0.7) 

-0.12 

(-0.5) 

107.2 

(1.7) 

.292 

(66)5 

15 

50x250 

Dry 

0 

2.5 

(1.0) 

2.39 

(0.9) 

2.79 

(U.2). 

0.27 

(1.1) 

126.2 

(2.0) 

287 

(57) 




1 

2.5 

(1.0) 

3.1 

(3.3) 

2.79 

(11.2) 

0.54 

(2.2) 

126.2 

(2.0) 

287 

(57) 




2 

3.6 

(1.4) 

17.3 

(6.8) 

1.87 

(7.5) 

-0.12 

(-0. 5) 

126.2 

(2.0) 

287 

(67) 




7 

5.3 

(2.1) 

24.4 

(9.6) 

1.17 

(4.T) 

-0.25 

(-1.0) 

145.1 

(2.3) 

287 

(57) 




35 

8.4 

(3.3) 

23.6 

(9.3) 

0.87 

(3.5) 

-0.25 

(-1.0) 

126.2 

(2.0) 

287 

(57) 




135 

10.9 

(4.3) 

21.1 

(7.9) 

0.60 

(2.4) 

0 

(0) 

132.5 

(2. 1) 

287 

(57) 

, 



335 

13.5 

(5.3) 

15.5 

(6.1) 

0.37 

(1.5) 

0.1 

(0.4) 

j 126.2 

(2.0) 

287 

(57) 

i 



435 

13.5 

(5.3) 

13.2 

(S.2) 

0,37 

(1.51, 

0.25 

(1.0) 

■126.2 

(2.0) 

287 

(57) 

: 

50x250 

Wet 

0 

2.3 

(0.9) 

1.3 

(0.5) 

4.0 

(16.1) 

1.69 

(6.8) 

88.3 

(1.4) 

289 

(60) 




2 

2.5 

(1.0) 

17.5 

(6.9) 

4.0 

(16.1) 

-0.50 

(-2.0) 

‘100.9 

(1.6) 

289 

(61) 




7 

2.8 

(1.1) 

23.9 

(9.4) 

3.9 

(15. 5) 

-O.SO 

(-2.0) 

j 107.2 

(1.7) 

289 

(61) 




1357 

3.8 

(1.5) 

3.3 

(1.3) 

2.2 

(9.8) 

0.37 

(1.5) 

82.0 

(1.3) 

289 

(61) 




1407 

3.8 

-L1.S) 

2.8 

11^1) 

2.1 

(8.5) 

0.70 

(2.8) 

82.0 

(1.3) 

289 

f61i 

17 No Data ! 





j 


— 

! 1* 

50x250 

Wet 

0 

2.5 

(1.0) 

1.3 

(0.5) 

— 

0.87 

(3.5) , 

1.04 

(4.2) 

0 

(0) 

290 

(63) 




2 

2.8 

(1.1) 

18.8 

(7.«) 

1.0 

(4. 1) ' 

-0.62 

(-2.5) 

0 

(0) 

290 

(63) ; 

j 



$ 

3.0 

(1.2) 

25.7 

(10.1) 

0.87 

(3.5) 

-0. 75 

(-3.0) 

0 
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2 
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Table 2-1, 


Refilling Test I>'t (Continued) 




|Run 

No. 

i Iniciai 

' Basket | Top 
. op J Screen 
Screen i Wetting 

Time. 

sec. 

Stan Basket 
Internal 
Liquid Level, 
cm (in) 

Start Basket 
External 
Liquid Level, 
cm (in) 



Plow Channel. 
AP. kN/m2 
(laHjO) 

Start Basket, 
AP, kN/m* 
(in HjO) 

Flow Rate. 

cmVaee 

(GPM) 

Start 
Basket 
Top 
Temp. . 
K (F) 

! 23 

30x250 

1 Dry 

0 

0 

(0) 

0 

(0) 

— 


0 

(0) 


(0) 

292 

(66) 

M 

Stand- 


2 

0 

(0) 

17.0 

(6.7) 

~ 


0. 25 

(1.0) 

0 

(0) 

292 

(66) 


Pipe 


4 

4.1 

(1.6) 

20.6 

(8.1) 



-0. 50 

(-2.0) 

0 

(0) 

292 

(66) 

i 



7 

4.8 

(1.9) 

23.6 

(9.3) 



-0.50 

(-2.0) 

0 

(0) 

292 

(66) 




108 

15.5 

(6.1) 

20.6 

(8. 1) 

- 


-0. 50 

(-2.0) 

0 

(0) 

292 

(66) 



i 

13S 

15.5 

L6. 1). 1 

20.1 

(7.9) 

- 


-0. 62 

(-2.5) 

0 

(0) 

292 

(66) 

I 23A 

5(h250 

! Wet 

0 

1.7 

(0.8» 

1.3 

(0.5) 

_ 


-0. 12 

(-0.5) 

0 

(0) 

292 

(66) 

! M 

Stand- 


2 

2.5 

(1.0) 

18.8 

(7.4) 

- 


-0.50 

(-2.0) 

0 

(0) 

292 

(65) 

i 

Pipe 

! 

3 

3.7 

(1.4) 

21.3 

(8.4) 

- 


-0.75 

(-3.0) 

0 

(0) 

292 

(65) 




7 

3.8 

(1.5) 

25.7 

(10.1) 

- 


-0.87 

(-3.5) 

0 

(0) 

292 

(65) 



; 

25 

4.2 

(1.6) 

25.7 

(10.1) 

- 


-0. 92 

(-3.7) 

0 

(0) 

292 

(65) 

1 



55 

4.1 

(1.6) 

24.4 

(9.6) 

- 


-0. 95 

(-3.8) 

0 

LO) 

292 

(65) 

24 

200><600 

1 Dry 

0 

0 

(0) 

1.8 

(0.7) 

0 

(0. 1) 

-0. 12 

(-0.5) 

0 

(0) 

294 

(69) 

M 


1 

2 

0 

(0) 

17.0 

(6.7) 

0 

(0) 

0.27 

(1.1) 

0 

(0) 

294 

(69) 



1 

3 

4.1 

(1.6) 

19.6 

(7.7) 

-0.1 

(-0.5) 

-0. 50 

(-2.0) 

0 

(0) 

294 

(69) 



i 

! 

6 

4.8 

(1.9) 

25.4 

(10.0) 

-0.2 

(-0.9) 

-0.47 

(-1.9) 

0 

(0) 

294 

(69) 



1 

35 

10.9 

(4.3) 

25.4 

(le. 0) 

-0.5 

(-2. 0) 

-0. 10 

(-0.4) 

0 

(0) 

294 

(69) 




75 

12.7 

(5.0) 

25.1 

(9.9) 

-0. 6 

(-2.5) 

-0. 07 

(-0.3) 

0 

(0) 

294 

(69) 




1 - - 

300 

14.2 

(5.6) 

24.6 

19.7) 

-0.5 

(-2.0) 

0.05 

0.2 

0 

(0) 

294 

(69)! 

25 

1 

200x600 

1 Dry 

0 

3.6 

(1.4) 

3.0 

(1.2) 

2.8 

(11.2) 

— 

0.1 

(0.5) 

157.7 

(2.5) 

296 

(73)1 

1 M 


t 

2 

4. 1 

(1.6) 1 

16.6 

(6. 5) 

2.3 

(9.11 

! 0.4 

(1.6) 1 

157.7 

(2.5) 

296 

(73)1 

1 


t 

! 

3 

5.6 

(2.2) 

19.8 

(7.8) 

1.7 

(6.9) 

-0.15 

(-0.6) i 

i 164.0 

(2. 6) 

296 

(73)1 

i 



6 

6.1 

(2.4) 

1 26.4 

(10.4) 

1 1.4 

(5.8) 

! -0.1 

(-0.5) 

170.3 

(2.7) 

' 296 

(73)1 



J 

! 61 

11.6 

(4.5) 

26.2 

(10.3) 

1 0.6 

(2.4) 

! 0 

(0) 1 

1 170,3 

(2.7) 1 

1 296 

(73)1 

i 



422 1 

1 14.2 

(5.6) 

16.3 

(6.1) 

0.3 

(1.2) 

0.1 

(0.5) 

157.7 

(2.5) 

' 29b 

(73)! 

1 

! 26 

200x600 

Dry 

1 

0 ! 

1 


(3.5) 

10.4 

1 

(4. 1) j 

1 . .. 

1 7.4 

(29.6) 

0.3 

(1.1) ' 

' 2523 

(40. 0) 

! 297 

l76)| 



, i 

1 ! 

* 3.5 

(3.3) 

17.5 

(6. 9) 1 

1 6.7 

(27.0) 

0.5 

(2. 2) 

i 2523 

(40. 0) 

1 297 

(76)| 

1 


. 1 

* 1 

10.4 

(4.1) 

26.9 

(10.6) 

1 

(25.9) : 

! 0.5 

(2.2) 

i 2965 

(47.0) 

! 297 

(75) 

j 


. 1 

34 I 

12.7 

(5.0) 

1 13.2 

(5.2) 1 

1 6.2 

(24.8) 

' 0.8 

(3.39 

' 2744 

(43.5) 

297 

(76) 




! 44 1 

1 12.7 

(5.0) 

' 5.S 

12.3) 

i 6-0 

(24.0) 

2.1 

(8.3) ^ 

' 2681 

(42.5) 

298 

(77) 

1 

: 27 

200x600 

1 Wet 

! 0 1 

j S.6 

(3.4) ! 

' 6 

(0) 1 

7..) 

l29.2) 

. 1. 1 

(4.6) 

2397 

(38.0) 

I 

1 297 

(75) 

M 


i 

1 - ' 

! 3.6 

(3.4) 

! 19.8 

(7.4) 

6.7 

(27.0) 

0.4 

(1.8) 

' 2523 

(40. 0) 

j 297 

(75) 

I 


j 

1 « i 

1 9.4 

(3.7) 

25.9 

(10.2) 

[ 6.6 

(26.5) 

0.5 

(2.1i ] 

1 2649 

(42.0) 

j 297 

(75) 



j 

i i 

! 12.7 

(5.0) 

, 10.9 

(4.3) ! 

6.3 

(25.4) 

' 1.0 

(4.2) 

i 2681 

(42.5) 

1 297 

(76)1 

i 
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I 6. 1 
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2.4 

; 
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7. Continue outflow from the basket (throu|^ the channels) until pullthroug^ occurs. 


Test data from Table 2-1 was studied for selection of test rues to be compared with 
REFILL program computer predictions. 

2.4 REFILLING MODEL CORRELATION AND VEHICLE PREDICTION 

The REFILL computer program was successfully correlated with test data. The 
computer program was then utilized to predict refilling for Centaur D-lS LH 2 and LO 2 
start baskets. 

2. 4. 1 REFILLING COMPUTER MODEL CORRELATION. Five test runs were 
selected for computer program correlatitm: test runs 13 to 16 and 22 as listed in 
Table 2-1. Several computer program Improvements were made at this point* as 
described in Section 2. 1* in order to accurately simulj^e the test runs. These modifica- 
tions were: addition of a table for pullthrough hei^t vs flow rate, an option to eliminate 
spilling at initiation of outflow, and an option to allow the standpipe to remain dry 
during refilling independent of the liquid level outside the basket. 

Computer runs were made with both air and ethanol vapor as the gas inside and 
surrounding the basket. Tie data was best correlated with air properties for vapor 
flow across the start basket. Tank liquid level vs time obtained from test data was 
an input for each computer run. Given the initial liquid level in the basket and the tank 
liquid level versus time the major variables used in obtaining the correlation were 
AST, NDR, NW' and DBF2. AST is the screen area of the standpipe. NDR is a flag 
that keeps the standpipe screen dry if not set equal to zero. If NDR = 0 the standpipe 
screen is wet. NW is a flag that determines wetting conditions of the uncovered 
screen excluding the standpipe. If NW = 1 wetting of the screen does not occur tdiove 
the liquid level outside the basket. If NW = 2 the basket is wetted to the base of the 
standpipe. (If NW - 3, Subroutine SWET is used to calculate screen wetting based on 
screen/plate wicking. ) DBP2 is the bubble point of the top or standpipe screen. (The 
top screen consisting of a single screen layer is treated as a "standpipe" in the 
computer simulations. ) Comparisons of test data and post test computer simulaticai 
are shown in Figures 2-11 to 2-13. For runs no. 14 and 15 the screen is initially 
dry, thus NDR = 1. Minimal screen wetting is assumed (NW = 1). The top screen 
bubble point is 65 microns for 50 x 250 screen. The standpipe area that remained 
dry during refilling was determined parametrically by fixing the initial basket level 
and varying the standpipe unwetted area until the basket level at the final time (46 
seconds for run no. 14 and 34 seconds for run no. 15) was matched. Results are 
plotted in Figures 2-12 and 2-13 with straight lines connecting the data points and 
simulation points. 

For runs no. 13, 16 and 22 the standpipe screen was initially wet (NDR= 0) and the basket 
screen was initially wet (NW=2). For runs no. 13 and 16 a 50x250 flat 5. 08 cm (2. 0 in) 
dia top screen was used as the "standpipe. " For run no. 22 a 3. 81 cm (1.5 in) high by 
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Figure 2-11. REFILL Computer Program 
Simulation of Test Run No. 13 



Figure 2-12. REFILL Computer Program 
Simulation of Test Run No. 14 
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REFILL Computer Program Figure 2-14. REFILL Computer Program 
Simulation of Test Run No. 15 Simulation of Test Run No. 16 
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Figure 2-15. REFILL Computer Program Simulation of Test Run No. 22 

5. 08 cm (2. 0 in) diameter cylindrical standpipe was used. Initial computer runs made 
for test runs no. 13, 16 and 22 indicated that no refilling should occur. Inspection of 
test notes and movies indicated that some leakage had occurred around the standpipe 
screen/main screen attachment ring. Computer runs were then made to determine 
this leakage area. From the initial runs with no Riling the minimum size leak that 
would produce refilling was determined to be approximately 130 microns. Parametric 
evaluation of the standpipe area, using a bubble point of 130 microns, determined the 
areas required to produce the desired refilling for each run. This evaluation was 
similar to that used for runs no. 14 and 15 with the initial basket level fixed and the 
area of the leak varied until the "Rnal" liquid level was matched. 

Although only five runs were evaluated in detail because of time limitations, these 
runs represent a wide range of geometry, wetting and outflow conditions. He 
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usefulness of the program In predicting test conditions qualified the program as a 
versatile and apparently accurate tool for predicting actusd start basket refilling 
under mission conditions. 

2.4.2 CENTAUR D-IS START BASKET REFILLING. Refilling analysis of the Centaur 
D-lS LO 2 and LH 2 start baskets was completed using the REFILL computer program. 
Three reference missions described in detail in Tables 1-1, 1-2 and 1-3 were used 
for capillary device design and performance evaluation. The missions were a single- 
bum planetary mission, a two-bum S3mchronous equatorial mission and a five-bum 
low earth orbit mission. Each bum for these three missions was examined as a 
possible worst case condition for refilling. Three burns were selected for analysis 
based on the refilling gravity level and engine bum time. The first bum of the two- 
bum and five-bum missions were selected because of their relatively low gravity refill 
conditions. The fourth bum of the five-bum mission was selected because it required 
refilling to be accomplished in the shortest time period. Refilling time was computed 
by subtracting the settling time under main engine thrust (after settling produced by the 
start sequence) from the main engine bum time. Settling time was computed using 
five times the free fall time with the start sequence and main engine thrust levels. For 
both of the first bum conditions examined (the two bum synchronous equatorial mission 
and the five-burn low earth orbit mission) the start sequence thrust was predicted to be 
sufficient to settle the propellants. This allowed the full bom time to be used for 
refilling. For bum four of the five-bum low earth orbit mission (bum time = 18. 9 sec) 
some settling with main engine thrust will be required. This reduces the allowable 
refilling time to 16 seconds for LH 2 and 16. 9 seconds for LO 2 . 

No liquid collection was assvuned to occur around the basket during settling of the tank 
contents. Liquid collection was assumed to occur instantaneously after settling was 
complete. Liquid refilling calculations commenced at this point. (The computer program 
can handle refilling calculations during settling and collection, however the conservative 
assumption of instantaneous collection simplified the calculations required to generate wiie 
input for the simulation. ) Schematics of the LO 2 and LH 2 baskets are shown in Figures 
2-16 and 2-17. The LO 2 basket was modelled as two layers of 50 x 250 screen around 
the entire screen surface with the exception of a 5. 5 cm standpipe which is a single 
layer of 50 x 250. The LH 2 basket was modelled as a double layer 50 x 250 surface 
with a 40 X 200 mesh 10. 7 cm standpipe. The standpipe was input with an 84 micron 
bubble point and the flow properties of 50 x 250 screen. The actual configuration 
consists of a 50 X 250 double screened lower surface and 40 x 200 double screened upper 
surface with a 40 v 200 single screened standpipe. An 18 x 18 mesh cross screen 
separates the 50 x 250 lower and 40 x 200 tqjper compartments. The assumptions 
used in modelling the LH 2 configurations were made because no data is available for 
40 X 200 screen flow/pressure drop and the REFILL program does not have the 
capability to handle multiple compartments. Initial basket liquid volume was based 
on subtracting the thermal conditioning usage betw'een bums and the liquid outflow 
from the basket prior to settling from *he start basket volume. 

Results of the analysis are shown in Figures 2-18, 2-19 and 2-20. The volume versus 
height relationships used for both baskets are shown in Figures 2-21 and 2-22. 

2-29 



BASKET LIQUID LEVEL, cxn (ft» 


i 



Figure 2-17. Schematic of LH£ Start Basket 
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Figure 2-19, Centaur D-IS LO 2 and LH 2 Basket Refilling, Bum One • 
Low Earth Orbit Mission 
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Figure 2-21. Centaur D-IS LO 2 Start Basket, Volume Figure 2-22. Centaur D-IS LH 2 Start Basket, Volume 
Versus Hei^t Versus Height 



The refill computer program with variable g-level capability was used for this 
analysis. For the three bums considered, the results indicate that start basket 
refilling will occur at the imposed g-level and within the refilling time constraints. 

2. 4. 3 CONCLUSIONS AND RECOMMENDATIONS. The computer program REFILL is 
an important capillary device design tool that permits reRlling time for start baskets 
(using settled fluid) to be predicted. The program as documented in Appendix A has 
been correlated vith ground test data and demonstrated to be a useful predictive tool. 
The program is recommended for use as a design tool to determine start basket volume 
and geometry. In fact, the program has already been used extensively at General 
Dynamics Convair to evaluate the capability of a start basket designed for ground 
testing with LH2 to be refilled with LN£. Test data to be obtained from the on-going 
program (NAS8-31778) will provide additional verification of program accuracy in 
predicting start basket refilling. 
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CAPILLARY DEVICE RETENTION A A POR INFLOW 


In order for cryogenic capillary devices to function properly between engine restarts, 
screen retention capability must be maintained by keeping screens in a wetted condition. 
For partial acquisition devices, such as the LO 2 and LH 2 start baskets designed for 
the Centaur D-IS in Reference 3-1 and 3-2, contact between the screen surfaces and 
main liquid pool will not be continuously maintained between bums. The screen 
surfaces will dry out due to heat input to the surface if liquid is not continuously 
resupplied to the screen. Supply of thermal conditioning liquid to the screen 
surface by use of wicking between multiple screen/plate combinations was 
discussed in detail in Reference 3-2. Analysis and testing performed showed that 
liquid could be supplied to the screen surfaces in sufficient quantity to prevent screen 
dry out during all possible heat flux, adverse acceleration and wicking distance 
conditions. 

One concern that was not within the scope of the NAS3-19693 effort reported in Reference 
3-2 was how to handle vapor inflow. Vapor must enter the start basket volume to 
replace the liquid being evaporated from the screen surface. Unless this vapor can 
be directed to specific areas of the start basket, the vapor will detrimentally affect 
wicking flow and screen retention capability. 

If vapor enters the multiple screen/plate wicking barriers, the path for wicking flow 
will be reduced to the point where liquid retention is lost. Figure 3-1 illustrates a 
typical multiple barrier wicking configuration. The gap between the innermost layers 
will be in the order of 0. 064 cm (0. 025 in) for low gravity application. This approach 
is satisfactory as long as sufficient liquid remains within the gap to provide wicking 
flow. However, eventually the gap will become filled with enough vapor that wicking 
flow rate is insufflcient to intercept heat input to the outer screen. Outer screen dry- 
out results and liquid in the channel is lost. Because the gap between the wicking 
barriers is small, only a small amount of liquid will be retained in the gap and the time 
required for screen dryout will not be significant compared to the time between bums. 

A concept which successfully avoids the adverse effects of vapor penetration is a 
multiple screen liner device with a window, shown schematically in Figure 3-1. The 
so-called ’’window" is a single screen that enables communication between the ambient 
and the inner volume of the surface tension device. This communication is established 
by providing a latter bubble point diameter screen for the window than for the multiple 
screens, to guarantee that vapor penetration will occur at the window screen and enter 
the device inner volume. The advantage of allowing vapor penetration at the window 
screen is that the multiple screens will remain wetted and not dry out. Consequently, 
heat interception capability of the device is maintained. 
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MULTIPLE SCREEN AND WINDOW BUBBLE POINT 
aP't. 


Figure 3-1. Multiple Screen Liner With Window 

Vapor inflow aorosa wetted screened windows was examined in detail during this study. 
An extensive program of small scale bench tests was conducted in order to verify and 
correlate the equations. Persistent difftculties were experienced in obtaining repeatable 
results. 

The remainder of this section discusses the equations believed to govern the oper^ion of 
a single window screen within a multiple screen barrier cimfiguration. Descrif^ionst 
drawings and photographs of the experimental apparatus and test specimens used are 
presented with a discussion of losuUs. analytical model correlations and recommenda- 
tions for additional experimentation. 

3. 1 VAPOR INFLOW ANALYTICAL MODEL DEVELOPMENT 

Analytical evaluation of the vapor inflow process was performed in conjunction with 
some simple bench tests to obtain empirical factors for developing the governing 
equations. A small transparent box, representing a multiple screen/window configura- 
tion, was tested in hexane for this purpose. The box is shown schematteally In Figure 
3-2 along with the nomenclature used in the analysis. Figure 3-3 Illustrates the 
interface conditions and pressures at the screened window. Using the nomenclature 
of Figures 3-2 and 3-3, the driving pressure for wieking across the window screen, 
^^wlck determined to be 
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‘'’wick = ■’Ljj ■ Pi. 

1P<t, - Pa - Pi. '»- 2» 

at tht oolnt of vapor breakthrough across the 
window screen. 

APwlck-PL2-<Pa*^V 

where “^vvlck** Equation 3-3 Is a 

maximum due to the minimum value of the 
liquid film pressure Pl being u.sed, APct 2 
Is the bubble point of the window screen. 

For wlcklng flow to occur along the window 
screen PL 2 ^ ~ ^^ 02 * 

The heat intercept capability of a single 
horizontal screen is found from Reference 
3-3 to be 

Q.' As - mX/Ag - pXf (C/l?) (o 


uqiiiii ‘’^wmixiw maw 

bVAPt^llATION SCIIEEN HCHKEN 



Figure 3-2. Multiple Screen Liner 
With Window 
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where 


Q is the applied heat rate 

As is the screen surface area 

p is the liquid density 

X is the latent heat of vaporization 
6 is the screen thickness 

C is a correlation constant 

L is the distance from the liquid source to the wicking front 
a is the liquid surface tension 

H is the liquid viscosity 

The derivation of Equation 3-4 is based on the assumption that A^Vick screen 
bubble point AP are equivalent. For the window screen application, however, these 
terms ai'e not equivalent and Equation 3-4 must be adjusted accordingly. The adjust- 
ment has been to introduce a correction factor, , which results in 


Q/Ag = m X/As = P X 6 (C/l2) (a/^) (F^ ) 

(3-5) 

where = APwlck^^^cg 

(3-6) 

i'rom Equation 3-5 one can show that 


Q/A oc F^/l2 

(3-7) 


The window screen will remain wetted as long as its wicking capability is equal to, or 
exceeds the incident heat flux, Q/A. Because liquid wicking capability is decreased as 
F(j (i.e. , APwick) reduced, the screen will begin to dry once F^ drops below a 
minimum level. Screen drying can be explained from Equation 3-7 which shows that 
as Fg. drops in value, L must also be reduced in order to maintain a constant Q/A. 

A reduction in L is equivalent to having a portion of the screen become dry because of 
the unavailability of liquid. Vapor flow will occur through the screen at this point. 

The flow of vapor into the inner volume will increase Py, PLg and F^. The increase 
in F(j will increase window screen wicking flow rate capability. Eventually the dry 
portion of the screen will re-wet, and vapor penetrrtion will cease. Thus the pressure 
difference across the window screen cycles between minimum and maximum values. 
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Using the nomenclature of Figure 3-2, where H3 is the maximum liquid head that 
can he retained by the window screen and is^ls value Just after vapor break- 

through 



PL2'PA-»g«3_^^ (3-9) 

’Pwlok ' Pa - P BH3„,„ - <Pa - : J - 0 , (8-10) 


and F. is found experimentally from 


F- 


(H3 - H3 . » ' H3 

’max min '’max 


(3-11) 


F^ is a quantity to be determined from testing and is a direct measure of the ability of 
the window screen to rewet and remain wetted. Substituting back Into Equation 3-3 
allows the heat Interception capability of the configuration to be determined. 


Also to be determined from the testing is a factor that measures when the screen 
cannot rewet. As the liquid level in the capillarv device under the window screen 
drops, the capability for wicking across the window screen is reduced. Experimentally, 
a relationship was found that predicted the limit of the window /main screen system. 

This equation 


^'^(max) 




(3-12) 


is a pressure balance at the time of screen breakdown. (The eiqiloratory tests conducted 
prior to formulation of the analysis and test plan had shown that the quantity was a 
constant equal to 0. 44 for the apparatus shown schematically in Figure 3-2 using 
hexane as the test fluid.) The following equations transform Equation 3-12 into an 
expression that can be evaluated from experimentally measured quantities. 

Combining Equation 3-3 with the relationship Ptj Pl^ -t- pgHj (from Figure 3-2) gives 

Pu - OgHi = (Pa - aPag) ^ APwick (3-13) 

Substituting Equation 3-12 into 3-13 results in 

l^A - l^u ^ ~ l^l ~ (1*A - APct 2^ APwlck (3“14) 

or 

Pg Hi = (1 - K„) (3-15) 
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Substituting Equations 3-8 and 3-10 into 3-15 


Hl= W'HW-H3„ta' 

or 

H, = Ho . - K„ Ho 
^ ’’min “ ’’max 

Dividing Equation 3-17 by 


(3-16) 

(3-17) 


(3-18) 


Since Equation 3-11 shows Ho . /Ho = 1 - Fa 

’’min -^max 

Equation 3-19 Identifies the maximum head difference that can exist between the inner 
volume and the outer volume liquid levels before window screen drying occurs. At this 
point Hi is replaced by Hiqj. to designate a critical head difference and Equation 3-19 
becomes 


HiCr/«3 - 1 - F„ - Ku (3-20) 

max 

The heat interception capability and Hicj. can be related by combining Equations 3-5 
and 3-20 


H 


ICr 


H, 


= 1 - K. 


’max 


- 

-pX6(C/L^) a/n . 


(3-21) 


The single unknown in Equation 3-21 is K^. An extensive program of normal gravity 
bench tests was conducted to determine values of Fa and Ku (using Equations 3-11 
and 3-20) for a variety of screen configurations and test fluids. This information 
would then be extr^olated to cryogenic fluids and low gravity start basket configurations 
based on the physical property and geometric relationships found. A corollary objective 
of the testing was to increase the basic understanding of the physical processes 
controlling the wetting of screens subjected to transverse vapor flow. 

3. 2 VAPOR INFLOW APPARATUS AND TESTING 

A series of tests were run in order to evaluate the equations presented in Section 3. 1. 
Initial testing with a single test article using hexane provided repeatable results that 
were used to formulate the analytical model. Additional tests using several test fluids 
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and test articles were planned based on the initial results. These results, however, 
were found to be disappointing both in repeatability of runs and rewetting performance. 
Considerable effort was expended in refining both test conditions and test article 
geometry in order to achieve results consistent with that obtained with the initial 
test article. Only limited success was achieved in this effort. 

The following section chronologically describes the test configurations and test 
ccmdltions employed during the study. Testing was conducted in five test series. The 
first test series used a single small test article and produced repeatable results. 

Based on these results six new articles were built and tested with hexane. The 
configurations were designed to yield parametric data on the screen wetting process 
when subjected to vapor flow. Results obtained were less repeatable than the first 
test series and showed poorer wicking performance. Some of the problems ware felt 
to be a result of (I) lack of control of wicking between the main screen and window- 
screen, (2) reduced volume under the window and (3) increased screen deflection 
because of the increased screen span between supports. New articles were fabricated 
with reduced span, improved main screen/window screen wicking and increased volume 
under the window. The third test series was conducted with hexane using a bell jar for 
improved environmental control with the new test article as well as the original article. 
Somewhat better results were obtained under these conditions. Two additional boxes 
were then fabricated with different window lengths and tested in the fourth test series 
using hexane, ethanol and Freon TF. The test data was not as consistent as that 
vybtalned in the third test series. In the fifth test series two of the boxes were modified 
and tested with hexane and Freon TF. One of the boxes used Teflon dams to promote 
unidirectional wicking. The other box used screens with higher retention capability 
than previously tested. 

Quantitative experimental results and the analytical model correlation are presented in 
Section 3. 2. Recommendations for additional testing are presented in Section 3. 3. 

3.2. 1 VAPOR INFLOW TESTING - INITIAL TEST ARTICLE 10/76 TO 11/76. Tests 
were initially conducted using the device illustrated In Figure 3-4, using hexane as 
the test fluid. Test observations for the first test (test no. 1) were run as follows. 

1. After the inner and outer volumes were filled with hexane, the fill valve was closed 
and the model oriented as shown in Figure 3-5. 

2. The vapor pressure of hexane was sufficiently high that evaporation readily 
occurred at ambient pressure and temperature conditions. Consequently, 
evaporation commenced immediately at the screen. 

3. Bubble penetration at the 50 x 250 n window screen was observed to occur at 20 
second intervals for flow periods of about 2 seconds (Figure 3-5a). ft was 
believed that liquid evaporation at the screen caused liquid flow from the inner 
volume to the outer volume throu^ a slot at the bottom of the device. This 
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-U) 165 ^ 500 U TWTLLBD-WEAVE 
DUTCH SCaEEN (BUBBLE 
POINT aP » 2. SI kN/m2 
(60 psf) 


(2) 50 X 250 M plain-weave 
SCREEN (WINDOW) 
(BUBBLE POINT AP . 1.44 
lcN/m2 (30 ptD 


NOTES: 

1. CONTAINER CONSTRUCTED OF 
0.32 cm (1/8 In. ) THICK PLEXIGLAS 
EXCEPT FOR 1. 91 cm (3/4 In. ) 
DIAMETER USED FOR WINDOW SCREEN. 

2. ALL DIMENSIONS ARE IN cm (in.) 


0.32 (1/5) X 2.54 (1) SLOT 


5.6(2.2) 

\* 3.8 (1.5) • 




Figure 3-4« Multiple Screen Liner With Window Test Device 



(a) Liquid eviqwrated at screen ts replaced by (b) Liquid leaka^ begins once level drops below 

liquid from Inner volume, which is replaoed base of screen* This Is evidence that screen 

bj' vapor penetration at window dry-out has occurred. 


Figure 3-5. 


Multiple Screen Liner With Window - Initial Test Series - 
Test No. 1 
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liquid flow reduced the Inner volume pressure suffloiently to cause bubble 
> pen^ration (about ten to twelve 0. 32 cm (1/8 In) diameter bubbles in two seconds). 

I Apparently the eighteen seconds between bursts of bubble flow was the liquid flow 

I time required to exceed the bubble point ^ capability of the window screen* 

4* Surface evaporation at the screen continued until the inner volume liquid level 
dropped below the screen at the base of device (Figure 3-5b). Once the screen 
dried out due to surface evaporaticsi. liquid retention was lost and liquid leakage 
began. The test was terminated at this time. 

5. By visual observation bubbling rates and bubble sizes were used to determine an 
ai^roximate liquid evaporation rate of 0. 008 om^/sec (1. 8 in^/hour)« which 
translates to a heat flux at the screen of 419. S watts/m^ (133 Btu/hr-ft^). During 
testing an air stream was directed at the screen surface for a short period. It 
was estimated that liquid evaporation rate was increased by a factor of five to ten. 

Screen drying did not occur in the time period observed. 

The next test run (test no* 2) provided the following observations. 

1. After the inner and outer volumes were filled with hexane, the fill valve was 
closed and model oriented as shown in Figure 3-6. 

2. While the inner volume liquid level remained above the window (Figure 3-6a)» 
bubble penetration at the SO x 250 H window screen was obsen'ed to occur at 20 
second intervals. This was an identical observation to that in test no. 1 



a) Small vapor pock«t U formod next to 
winciow aa vapor penetratea window 
screen. 


b) Li^itd fllni La maintained at 
window by wtoklng along 
window aoreenfrom channel. 


0) LUmid leakage beglna onoe level 
drops below base of screen. This 
U evidence that screen diyout has 
occurred. 


Figure 3-6, Multiple Screen Liner With Window - Initial Test Series - Test No. 2 
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3. Surface evaporation at the screen continued even after the liquid level dropped 
below the window screen. Bey(»id this time* a liquid film was maintained (m 
the window screen by wicking from the outer liquid volume (Figure 3-6b). 

4. The test was terminated after the liquid level had dropped below the screen at the 
base of the device (Figure 3*6c). Liquid leakage began at this time, indicating 
that screen dry-out occurred. 

5. Liquid ev^oration rates wera the same as for test no. 1, approximately 0. 008 cmV 
sec (1. 8 inVhr). As with test no. 1, evaporation rates were increased for short 
time periods by directing an air stream at the screen suiiace. 

These tests were strong evidence that use of a window with a multiple screen liner 
could provide high heat flux interception capability, minimize wicking distance to 
prevent screen dryout and preferentially allow vapor to penetrate the capillary device 
while maintaining screen retention ciq>ability. 

Additi<»al tests in this initial test series were conducted to gain further insight into 
the mechanism of vapor penetration. Tests were performed to determine internal 
fluid pressure during liquid flow from the outer volume (main screen) to the window 
screen. A manometer was connected to the inner volume under the window to allow 
visual determination of the article internal pressures. Test observations for the 
subsequent test (test no. 3) were as follows. 

1. After the inner and outer volumes were filled with hexane the model was oriented 
as shown in Figure 3-2. 

2. Evaporation commenced immediately at the screen surfaces, which resulted in 
liquid flow from the inner volume to the outer volume. This flow demand resulted 
in manometer level drop as the inner volume pressure decreased. 

3. Once steadj^state conditions were established, bubble penetration of the 50 x 250 
mesh window screen occurred at about 180 second intervals for flow periods of 
about 7 seconds. 

The manometer level increased during vapor penetration to reflect increasing 
inner volume pressure. 

4. Figure 3-7 describes the observed liquid head variations with time. Note that 
vapor penetration of the window screen commenced each time that liquid in the 
manometer dropped 11. 5 cm (4. 53 inches) below the fine mesh screen. Note 
further that vapor penetration ceased each time that liquid in the manometer 
rose to within 4. 95 cm (1. 95 inches) of the inner volume liquid level. 
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Figure 3-7. Manometer Liquid Head Readings During Test 

5. Figure 3-8 describes inner volume ullage pressure, P^, and outer volume liquid 
pressure at the screen, Pl 2 * These pressures are referenced to ambient 
pressure. Pa, and are given in terms of window screen bubble point AP, ^ 02 * 
which is equivalent to 11. 5 cm (4. 53 inches) of liquid hexane. 

6. The test was terminated when reached 4.32 cm (1.7 inches). The screen device 
was still performing satisfactorily. Liquid evaporation rates were the same as 
found in tests no. 1 and 2, approximately 0.008 cm^/sec (1.8 in^/hr). 

Several conclusions were drawn from Figure 3-8. Fluid pumping from the outer voluir j 
to the window is a function of Pu - Px^g defined in Figure 3-3. This quantity is a 
function of the window screen bubble point APa 2 shown in Section 3. 1. Window 
screen retention capability must be less than the main screen retention capability. 

Based on these initial tests, the analysis presented in Section 3. 1 was developed. A 
comprehensive test plan was prepared (Reference 3-4) based on this analysis. The 
operation of the test article was demonstrated to NASA/LeRC personnel. After 
presentation to NASA/LeRC, the test plan was approved and the test configurations 
were fabricate^. 

3.2.2 APPARATLS FABRICATION AND VAPOR INFLOW TESTING OF SECOND SET 
OF TEST ARTICLES - JUNE-AUGUST, 1977. Five test articles were fabricated 
according to the drawings of Figures 3-9 and 3-10. Article characteristics are 
tabulated in Table 3-1. A schematic of the test apparatus used is shown in Figure 3-11. 
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TIME (NO SCALE) 

Figure 3-8. Liquid Internal Pressures During Test 
Table 3-1, Surface Tension Device Configurations 


(1) 

Config- 

uration 

Upper 

Box 

Ass’y 

Window 
Length 
cm (in.) 

Window 

Screen 

Mesh 

Top 

Screen 

Mesh 

Window Screen Wink- 
ing Direction in 
Relation to Warp Wire 

A 

-1 

1.27 (0.50) 

50 X 250 

165 X 800 

Perpendicular 

B 

-2 

2. 54 (1. 00) 

50 X 250 

165 X 800 

Perpendicular 

C 

-3 

S. 08 (2. 00) 

50 X 250 

165 X 800 

Peipendicular 

D 

-4 

5.08 (2.00) 

50 X 260 

165 X 800 

Parallel 

£ 

-5 

1.27 (0.50) 

165 X 800 

200 X 1400 

Perpendicular 

F 

-6 

2.54 (1.0) 

166 X 800 

200 X 1400 

Perpendicular 

(1) Only one lower box will be required, 

The bottom screen will 

be 200 X 1400 mesh. 
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Figure 3-10 Surface Tension Device/Lower Box 




BBLLJAn 



Figure 3* U. Vapor Inflow and Test Schematic 





The test plan of Reference 3>4 describes the test conditions planned for evaluating 
screen dryout under heating conditions and determination of Ku under zero heat flux 
conditions. 

3. 2. 2. 1 Test Apparatus. In order to achieve near uniform evapoi'ation at the window 
screen a bell Jar is used to contain the test article and a vacuum is pumped on the 
bell Jar atmosphere at controlled flow rates. The flow rates and pressures can be 
used to compute equivalent incident heat flux conditions to the screen surface based 
on test fluid vapor pressure. The test article assembly consists of an upper box and 
a lower box which are Joined together at the inlets and outlets by tubing. Each box 

is constructed of Lexan and metallic Dutch weave screen. Details of the upper box 
are given in Figure 3-9. The six configurations of Table 3-1 are identical in every 
respect excei^ for the top screen and window screen mesh and window screen length. 

A single lower box configuration common to all vq;>per boxes was assembled according 
to the details of Figure 3-10. The box is constructed of plexiglas with 200 x 1400 
mesh screen over the bottom port. 

3. 2. 2. 2 Instrumentation. The liquid heads (Hi and H3) were determined by visual 
observation and were measured from the top plane of the upper box to the respective 
levels in the sight tube and manometer. An etched scale was included in the 
apparatus for each liquid head measurement. 

The test setup incorporated two methods of measuring liquid evaporation. One 
method used a Flow Technology Omniflo turbine flow transducer. This type of transducer 
features interchangeable orifices which allows use of the same transducer for up to six 
over-lapping flow ranges. For the tests using hexane and Freon as the test fluids, the 
range of evaporation rates was expected to vary between 0. 01 and 0. 10 ACFM (actual 
cubic feet per minute). Two transducer orifice sizes were selected; one for the 0. 01 
to 0. 05 ACFM range and one for the 0. 05 to 0. 10 range. 

The flow transducer installation used is shown in Figure 3-11. In order to avoid problems 
associated with saturated vapor flow rate measurement, gas temperature could be 
increased by passing it through a hot water heat exchanger as illustrated in the test 
schematic. The water temperature could be maintained at approximately lOOF. The 
flow transducer output was monitored on a Flow Technology MdL PRl-102 Flow Rate 
Indicator and also recorded on an oscillographic recorder. The downstream temperature 
and pressure was also oscillographlcally recorded for flow correction to volume ^CFK'<) 
or mass flow units. Because most tests were run without using the vacuum pumps, this 
measurement method was limited in use and the method below became the primary method. 

The primary method for determining evaporation rate at the screens was by measuring 
the inner volume liquid level. Referring to Figure 3-11, both the inner volume of the 
upper box and the sight tube were calibrated sxich that mass or volume can be determined 
directly from inner volume liquid levels. By noting time intervals between level measure- 
ments average mass evaporation rates were known for specified intervals during the tost. 
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The bell Jar chamber pressure was monitored using a 9tatham strain $^ge type absolute 
pressure transducer. This pressure and the flowmeter outlet pressure were osclllo- 
graphically recorded. 

Chromel/Constantan thermocoiq)lea were used to monitor temperatures In the test 
article as shown in Figure 3>9. Humidifier reservoir and bell Jar atmosphere 
temperature were also measured. 

3. 2. 2. 3 Test Procedure. The oonflguiatlcms shown In Table 3-1 were tested. A large 
matrix of tests were planned as documented In Reference 3-4. Only a small fraction of 
these tests were actually run as a result of difficulties with repeating the data. Zero 
heat flux tests were planned to be run by manually increasing the separation distance 
between filled )x>xes to as much as 25. 4 cm (10 In) until vapor breakthrough at the upper 
box window occurs. The test specimens were then to be reassembled with rigid glas*) 
cubes, filled with test fluid and tested In the bell Jar using the vacuum pump. Liquid 
head data was to be recorded until loss of retention at the lower box screen. 

The testing actually proceeded as described In the following narrative. Assembly of all 
seven boxes was completed and leaks repaired as required. Two of the boxes (the lower 
box and uj^r box no. 4 from Table 3-1) are shown in Figure 3-12. The two boxes 
mounted in the test setup are shown In Figure 3-13. The overall test setup is shown In 
Figure S-14. 

Initial vapor Inflow checkout runs were made with ethancl, hexane and Freon TF using 
several of the boxes listed in Table 3-1. These runs illustrated that the concept funct- 
ioned in a reasonable manner, but the results differend from those obtained in the initial 
tests (see Section 3. 2. 1). Pressure recovery within the window area was much less 
than previously experienced. Figure 3-15 shows results obtainod with the small box 
used in initial testing and the results obtained in this second series of testo. The 
difference in results is believed to be partially due to the greater screen deflection with 
the boxes described in Figure 3-9 and 3-10. 

To compare these results with those obtained in the initial test series, the small article 
tested in the initial te t series was returned by NASA/LeRC. The article had been 
damaged and repaired at NASA/LeRC. The window screen as received at General 
Dynamics Convair was clogged and was replaced with new 60 x 250 screen. A series 
of tests were run with both technical grade and reagent grade hexane. Screens were 
cleaned and tested. A set of tests was also conducted with both a new top screen 
(154 A 800) and window screen (50 x 250). Results indicated that original observations 
could be duplicated but pressure fluctuations and screen breakdown were erratic. This 
was felt to he caused by the sensitivity of the rewettlng process to the space for wicking 
between the window screen and the top screen. This spacing can change during a test 
due to pressure fluctuations and heating variations. The spacing also changes between 
test articles due to the difficulty in exactly repeating the screen geometry during 
attachment. 
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Figure 3“ 12, Lower and Upper Box (Typical) 
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Figure 3-13, Vapor Flow Test Model Setup 
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RESULTS OBTAINED IN INITIAL TEST SERIES 
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Figure 3-15. Vapor Inflow Test Results 

Tests were then run using top boxes nos. 2 and 3 with the lower box (similar to those 
shown in Table 3-1 and Figures 3-9 and 3-10). The top boxes had 165 x 800 main 
screen and 50 x 250 window screen. The lower box had 200 x 1400 mesh screen. Runs 
were made using hexane with either small or large tubing connecting the upper and 
lower boxes as shown in Figure 3-16. A total of 17 test runs were made. All runs 
were made with the boxes exposed to ambient conditions. 

Several observations were Immediately evident from the runs; con jtlon of water on 
the apparatus affected test results (results were Improved with Incr. aed condensation). 
This is probably due to the water reducing heat input. Fewer pressure cycles were 
obtained per test using the small tubing compared to using the large tubing between 
boxes. This is because of the increased volume in the window section using large tubing. 
Resealing of uppc .■ box no. 3 was less rapid than the box no. 2 tests. This is because 
the window length for box no. 3 is 5. 08 cm (2. 0 in) compared to a length of 2. 54 cm 
(1. 0 in) for box no. 2. Additionally, tests run with both technical grade and reagent 
grade hexane gave similar results. This Indicated that technical grade was sufficiently 
pure for testing purposes. 
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h* - 10.2 cm (4 la.) 

A) SMALL WINDOW VOLUME B) LARGER WINDOW VOLUME 

(MANOMETER PORT AND PORTS CONNECTING MAIN SCREEN TO BOTTOM BOX ARE NOT SHOWN.) 


BOXES ARE LEXAN WITH STAINLESS STEEL 50 x 250 MESH WINDOW SCREE»^165 x 800 MESH 
-5CREEN FOR THE REMAINDER OF THE UPPER BOX AND 200 x 1400 MESH SCREEN FOR THE 
LOWER BOX. 


Figure 3-16. Large Box Test Configurations 

The data obtained was reduced and compared to the analytical model presented in Section 
3. 1. Evaluation of the data centered on solving Equations 3-11 and 3-19 for the quantities 
Fa and Ky. Data was reduced tc determine HiCr» H3max ^3min each run and is 
presented in Table 3-2. Considerable variation was found in for a given 

run. (In runs made in the initial test series, was found to be fairly 

constant). Values of Hjcr did not approach the value anticipated based on Task n testing. 
(Hicr was generally in the 1. 5 to 2 cm range where it had been projected to exceed 4. 3 cm). 

Examining the test data and test procedures used, several possible explanations were 
found for the anomalous results. The variation in during a run was 

likely due to environmental conditions such as changes in heat flux, water condensation 
and fluid temperature. This could be eliminated using the bell Jar and other environ- 
mental controls such as outlined previously and shown in Figure 3-11. The observation 
that Hi did not reach its full value with the larger boxes was possibly due to changes in 
internal pressure differential causing wicking potential to decrease. The lower 
pressure differential across the screens causes screen spacing to increase when 
occurs. This resulted in a longer wicking distance and lower driving pressure for 
wicking. The recommended solution to this problem was to systematically spot weld the 
window screen to the main screen in the region of overlap and wicking. 

The larger boxes 12. 7 cm x 12. 7 cm (5 in x 5 in) exhibited considerable screen 
deflection. Volume under the window was small, allowing only a few vapor break- 
through cycles to occur before reached its critical value. In the third test series 
wicking variations were controlled by controlling screen deflections, screen spacing 
variations and wicking length. Environmental conditions such as temperature, 
humidity (water), test fluid partial pressure, vapor velocity around the boxes and test 
fluid purity were controlled by using a bell jar, heaters and other equipment described 
in Figure 3-11. The test fluid was filtered prior to use. 
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Table 3-2. Summary of Runs Made During Later Portion of Second 
Test Series 


Run Number 

1 

1 

a 

n 

K 

'a 

a 

n 

X 

1 

u 

X 

No. of Cycles 

's 

cl 

X 

's s 

n 1 

a 

o 

f-H 

X 

d 

X 

Comments 

lA 

11.7 

7.54 

4.7 

19 

0.34 

0.402 

0.26 

Original Box 



7.7 

1.6 

5 

0.36 

0.132 

0.51 

Box No. 2 



9.9 

3.4 

13 

0.20 

0.276 

0.52 

Large Tubing 



12.3 

2.2 

26 

0.09 

0.159 

0.75 


1C 

12.7 

6.4 

2.0 

1 

0.5 

0. 16 

0.34 

Box No, 2 

2C 

12.7 

9.8 


1 

0.23 

0.126 

0.65 

Small Tubing 

3C 

13.6 

9.5 


2 

0.30 

0. 093 

0.61 


4C 

13.7 

8.4 


1 

0.39 

0.120 

0.49 


5C 

13.5 

9.9 


7 

0.26 

0. 113 

0.63 


6C 

13.7 


1.8 

5 

0.21 

0.130 

0. 66 


ID 

13.7 



1 

0.333 

0. 120 

0.55 

Box No. 3 

2D 

11.6 

7.0 


2 

0.396 

0.143 

0, 47 

Small Tubing 

3D 

13.5 

10.0 


1 

0.254 

0. 132 

0.62 


IE 

12.8 

8.1 


3 

0.37 

0.267 

0.37 

Box No. 3 

2E 

13.3 

9.1 

1.8 

1 

0.314 

0.133 

0.56 

Large Tubing 

3E 

13.3 

6.4 

2.0 

3 

0. 52 

0.152 

0. 33 



3.2. 3 VAPOR INFLOW TESTING - THIRD TEST SERIES, AUGUST 1977 TO OCTOBER 
1977. A new test article was fabricated to the configuration shown in Figure 3-17. 

While these articles were being fabricated tests were run with the box used in the 
initial testing to illustrate the effect of the bell jar on controlling test conditions. The 
test article was filled with hexane, mounted on the test stand and covered with the bell 
jar. Four test runs were made. Cycling time was significantly longer for runs within 
the bell jar than with ambient runs, increasing up to 15 minutes per cycle as the bell 
jar became saturated with hexane. 

Vacuum pumping of the bell jar was used to increase the heat flux to the screen surfaces. 
Rapid pressure decay produced bulk boiling under the main screen. Pumping down 
caused the pressure within the box to be above the pressure In the bell jar. The place- 
ment of the top of the manometer below the level of the top of the box was used to relieve 
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HEATING ELEMENT 



Figure 3-17. Box Design to Limit Deflection and Control Screen Spacing 

the box internal pressure by overflowing from the manometer. Reducing the pressure 
slowly to 68. 9 kN/m2 (10 psi) allowed runs to be made at an Increased heat flux without 
causing bulk boiling in the main chamber. Data was fairly consistent at this pressure, 
however,H 3 minimum was below the level of the bottom of the box. (This caused spilling 
from the box). The tests indicated that the bell jar was satisfactory in controlling 
environmental conditions and that increased heat fluxes could be obtained by vacuum 
pumping, however, careful control of system pressures is required to obtain acceptable 
data for simulated heating conditions. 

The new box configuration for the third test serles,shown In Figures 3-17 and 3-18, 
was fabricated and tested using hexane under ambient conditions and in the bell jar 
under both atmospheric pressure and reduced pressure conditions. Atmospheric 
pressure data was more consistent than data obtained in the second test series. 
stayed above 7. 62 cm (3 In) until fell below 1. 91 cm (0, 75 in), (Since the box is 
7. 62 cm (3 in) high, dripping from the 200 x 1400 screen occurred if was less 

than 7. 62 cm (3 in). ) Subsequent tests showed the window screen remained wetteu 
imtil the outer and inner compartments were almost completely empty. This box had 
a 2. 54 cm window length. 



Test Configuration 






Pump down tests were run at low pump down rates at pressures ranging from 101 
kN/m^ (14. 7 psla) to 55 kN/m^ (8 psla). During all runs, vapor appeared under the 
main screen when the pressure was reduced below 90 kN/m^ (13 psla). Some vapor 
may have been attributable to dissolved air coming out of solution, and some may have 
formed by bulk boiling. 

Based on the favorable results oMahied the npw b<%, twp additional boxes with 
window lengths of 1. 27 cm (0. 5 in) and 3.81 cm (1.5 in) were fabricated using an 
ethanol compatible adhesive. Testing was performed using ethanol. Freon and hexane. 

3.2.4 VAPOR INFLOW TESTING -FOURTH TEST SERIES, DECEMBER 1977. Two 
additional boxes were fabricated similar to the configuration shown in Figure 3-17. 

The three configurations tested in the fourth test series are summarized in Table 3-3. 
Tests were run using hexane, ethanol and Freon- TF. 


Table 3-3. Vapor Inflow Test Configurations - Fourth Test Series 


Config. 

Window Length 
cm (in) 

Window Screen 
Mesh 

Main Screen 
Mesh 

Bottom Screen 
Mesh 

I 

1.27 (0.5) 

50 X 250 

165 X 800 

200 X 1400 

n 

2.54 (1.0) 

50 X 250 

165 X 800 

200 X 1400 

1 ^ 

3.81 (1.5) 

50 X 250 

165 X 800 

200 X 1400 


Data was obtained using the bell jar for tests run without heating and for those tests run 
using a 150 watt lamp for heating the top screen. One test was run using a blow dryer 
for convective heating without the bell jar. All tests run with an external heat input 
resulted in fluid dripping out of the bottom screen or complete screen breakdowr. As 
described in the third test series, dripping of fluid out of the bottom screen will occur 
if the pressure at the box bottom exceeds the ambient pressure at vapor breakthrough. 
A chronological summary of the tests run during the fourth test series is shown in 
Table 3-4. 

3.2. f) VAPOR INFLOW TESTING - FIFTH TEST SERIES, FEBRUARY 1978. Boxes 
no. 1 and no. Ill were altered in preparation for additional testing. The screens of box 
no. Ill (see Table 3-3)were replaced with 200 x 1400 window screen, 250 x 1400 main 
screen and 325 x 2300 bottom screen, (lliese screens are similar to those typically 
used for start basket ground testing. ) For box no. I, an attempt was made to restrict 
wlcklng to the window length (from the window/main screen overlap towards the front 
of the box) by eliminating other sources of wicking. A pair of Teflon strips were 
treated and bonded along the main screen/ window screen as shown in Figure 3-19. 

A series of tests were run with the new boxes. The boxes are designated as shown 
in Table 3-5. 
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Table 3-4. Test Lo& - Fourth Test Series 
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Bell jar and 12/14/77 Results similar to those described in run 21 for Box I. 

ambient air 




Table 3-4. Test Log - Fourth Test Series (Continued) 



3.81 cm (1.5 in) 200 x 1400 250 x 1400 325 x 2300 




Figure 3-19. Bux Design to Control Direction of Wicking 
(Box No. IV) 

Initially box no. V was run with Freoti TF in the bell Jar under ambient heat flux 
conditions. The general observation for the six runs made was that only one full cycle 
down to the window screen bubble point of 24 cm (9. 4 in) could be obtained. Vapor 
breakthrourh and pressure recoverj' was not reproducible. Twice the pressure hovered 
at 1 8 cn« (fi. 3 in) before large amounts of vapor broke through. In all cases there was 
an intemier'iaii- cycle from close to the top of the Dux to about 6 . C cm (2. G in), then 
both the wii.dow ar-^ main screen dried out. 

Box no. / »>. ac ul'^o teeted with hexane. Since the 200 x 1400 window screen could hold 
more than GO cm of he nd, the boll Jar was too small to contain the apparatus during 
these tests. The single ;c.s( ren was therefore accomplished in ambient air. For this 
test, the Irst cycle went down to an HSmax ^ cycling then occurring 

over a narrow ran^je betw^^en 50 cm (19.5 in) and 57 cm (22.2.^ in). The first cycle took 
25 minutes. Js;;bseqccnt cycling between the levels noted occurred for more than half an 
hour until was greater than 3. 3 cm (1. 3 in). At Hj = 4. 1 cm (1 . 6 in) the pressure 
(H 3 ) was down to 11. 7 cm (4 . 6 in). Breakdown occurred at an Hx level greater than 4. 1 
cm ( 1. 6 in). 

A series of tests was then run using box no. IV with Teflon "dams'* installed as shown in 
Fig^tre 3-19, using hexane in the bell Jar under ambient heat flux conditions. Two tests 
were run to breakdown. For the first test, which lasted four days, complete breakdown 
occurred at Hx (the level in the window area below the top of the box) « 4. 1 cm (1 . 6 in). 
Dripping commenced when Hx was approximately 2. 3 cm ( 0 . 9 in). For the second 
complete run (which ran eight days) screen breakdown occurred when Hx was 3.0 cm 
(1. 2 in). Dripping occurred when Hx was below 1. 27 cm (0. 5 in). 
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Several tests were run with box no. n in preparation for heat flux testing with hexane. 
Box no. 2 is the box with the 50 x 250 window screen, 165 x 800 main screen and 2.54 
cm (1 In) window length. A test was run in the bell Jar under ambient conditions. 

Cycling without dripping occurred until Kx was below 3. 7 cm (1.45 In). Breakdown 
occurred at Hx below 3.9 cm (1.55 In). (Breakdown level occurred during the evening 
and was not observed. ) Tests were then run with a vortex tube, using the hot side air 
to heat the screen. Ambient temperature was 298K. Hot side vortex tube temperature 
was in excess of 333K. The vortex tube was attached to a ring stand directly over 
the screen. Tests were conducted In ambient air under a vented hood. Temperature 
Just above the surface of the screen was measured with a thermometer. Two tests 
were craiducted with the vortex tube. The first test started with the hot end of the vortex 
tube approximately 2. 54 cm (1 in) above the screen. (The outlet of the tube hot side was 
a 2.54 cm (1 in) AN elbow.) After ten minutes, with the Hx at 1.8 cm (0.7 in) and 
cycling occurring between 7. 1 cm (2.8 In) and 8. 1 cm (3.2 in), the vo 'ex tube outlet 
was moved up to 7 . 6 cm (3 In) away from the surface of the screen. Temperature of 
the impinging gas above the screen surface was 311K. After 5 minutes, with rapid 
cycling continuing between 7. 1 cm (2 . 8 in) and 8 . 1 cm (3. 2 in), the vortex tube was 
moved up to 2. 54 cm (1 in) above the screen with Hx - 2. 3 cm ( 0 . i) in) . Ambient 
temperature at the screen was 61C. Cycling continiied between the same limits for 
about 2 minutes, then cycling occurred between 6. 6 cm (2. 6 in) and 7. 1 cm (2.8 in) 
for 7 minutes from t = 18 minutes until t - 25 minutes. At t = 25 minutes, Hx was 

5.3 cm ( 2.1 in). Cycling continued with H 3 hovering around 7 . 6 cm (3 in) as the entire 
window area drained. Vapor could then enter the main compartment. Window screen 
dryout occurred with the liquid level a fev tenths of an inch below the top of the main 
compartment. 

Another run was mad.n vMth hexMe and 'he vortex tube with box no. V. The run was 
continued until the window screen dried out. Two cycles were run under ambient 
conditions. The .ortex tube was then placed 2.i--i cm (1 In) above the window area. 

Rapid cycling occurred during the test with heing as low as 5 . 6 cm (2. 2 in). 

Generally, however, performance was similar to that ot the p* evlous test with the 
window area draining completely without breakdown and liquid level under the main 
screen dropping to below 1. 0 cm ( 0 . 4 in) of the bottom of the box without breakdown. 

This concluded the vapor inflow testing perfomed under this contract. 

3.3 ANALYTICAL MODEL CORRELATION 

Test results obtained in the five test series were examined, reduced and con'.pared 
to the analytical models presented itt Section 3. 1 . Unfortunately,for this data analysis 
no quantitative data was available to evaluate K^ or Fj for the first test series and the 
early runs of the second test series. Data from test series 2,3,4 and 5 'vas avail- 
able for analysis and correlation. Evaluation of the data centered on solving Equations 
3-11 and 3-19 for the quantities Fa and Ku. Data evaluation performed for the second 
test series has already been presented In Table 3-2. Similar data rcductloi,- was 
performed tor test series 3, 4 and 5 and is presented In Table 3-6. 
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Each line of Table 3-6 represents the results of a number of runs, as shown In Column 
7. For each run, many cycles were run and was determined for each cycle and 
averaged to produce an average value of F^ for that run. The data from each run was 
then taken and grouped with other similar runs in that series. For example, with box 
no. II using hexane and 50 x 250 windcw screen under ambient conditions in test serf es 
No. 3, six runs were made. The data shown in the table represents the minimum, 
maximum and nominal values of Fq and K^ for these six runs. 

A set of runs was made with a heat lamp for heating the screen in series no. 4 with 
box no. I in the bell jar. In series no. 5 a vortex tube was used to heat the screen 
directly. 

Inspection of the data shows that F^ and K^ vary significantly. There does not appear 
to be any correlation between window length and F<j or K^. For example, for bell jar 
runs made with hexane and 50 x 250 window screen, for the 1. 27 cm window, = 

0. 35; for the 2. 54 cm window, F^ = 0. 34; and for the 3. 05 cm window, Fj = 0. 39. 

For the ambient runs, however, F(j was 0.38, 0.21 and 0.45 for the 1.27, 2.54 and 
3. 05 cm window lengths respectively. 

Using Freon as the test fluiu, higher values of Ku and Fa were obtained for the 200 > 
1400 screen (box no. V) than f u- the 50 x 250 screen (box no. HI). Results with 50 x 
250 screen were consistent from box no. I to box no. III. Tests with 50 x 250 screen 
were not run in a bell jar because breakdown occurred too q. ickly for the boxes to be 
placed in a bell jar. These results are therefore not directly comparable because 
the 200 X 1400 screen/Freon tests were run in a bell jar. 

The results shown in Table 3-6 give a general range of Fq and Ku that can be used fo^ 
preli ninar:* de Ign purposes. Because of the lack of clear trends, this information 
should be used conservatively. 

Heat flux was applied to box no. I using a 150 watt spotlight with the box in the bell jar. 
Calculations indicated that a heat flux of approximately 0. 02 watt/cm^ was incident on 
the screen surfaces. Wicking calculations showed that the main screen wicking 
capability was less than 0. 01 watt/cm^ and the window screen capability was greater 
than 1 watt/cm^. Theory would indicate that the box should remain wetted and 
function properly. In fact, the window screen dried out shortly after the lamp was 
turned on. This may have been due to radiant beat transfer directly into the sides of 
the transparent box causing vapor to form in the box at screen/box intersections or in 
other areas that would cause precipitous failure. 

Heat flux runs were made with box no. n in test series no. 5 using a vortex tube. The 
vortex tube flowed warm air at 333K directly onto the screen surface. Heat flux 
incident on the screen averaged 0. 07 watt/cm2 and peaked at 0.39 watt/cm^. Pressure 
cycling of the box occurred on a regular basis with no apparent degradation due to heat 
flux. In two runs, draining of the entire window volume was possible before the window 
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Table 3-6. Vapor Inflow Test Summary 
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screen dried out. The wlcklng capability of the window screen was equivalent to 64 
watt/cm^. 

The nonrepeatability of the data as well as the importance of understanding vapor inflow 
for start basket thennal conditioning require that more work be done in this area. 

Work done under the present contract has produced some test results that appear 
encouraging for future efforts in this area. Better control should be sought in 
providing cleanliness of the apparatus and test fluid, controlling screen geometry, and 
controlling environxr^ntal variables such as temperature and humidity. A description 
of recommended action is presented in Section 3. 4. 

3. 3. 1 CENTAUR D-IS CAPILLARY DEVICE THERMAL CONDITIONING. Thermal 
conditioning calculations were performed to determine whether the window screen on 
the LO 2 and LH 2 basket would dry out when subjected to heating. Calculations were 
made to determine the values of the parameters, F^, and Ku (defined in 

Section 3. 1) that would be required to permit the • dndow screen to function properly 
over the entire set of mission conditions. Values of Fg. were determined by calculating 
the heat flux interception capability of wicks in low gravity using the equations formulated 
for screens in Table 2-2 of Reference 3-2. Screen properties were obtained from 
Reference 3-3. The maximum wicking distance for the standpipe windows was taken 
to be the standpipe height plus one half the standpipe diameter. For the LO 2 start basket, 
standpipe height is 5. 52 cm (2. 17 in) and standpipe diameter is 2. 10 cm (0. 828 in). For 
the LH 2 start basket, standpipe height is 10. 98 cm (4. 32 in) and standpipe diameter 
is 4. 14 cm (1. 63 in). Heat flux incident on the standpipe region was obtained from 
Table 2-1 of Reference 3-2. The maximum heat flux interception requirement was 
4.41 watts/m^ (1.4 Btu/hr ft2) for LO 2 and 2.52 watts/m^ (0.8 Btu/hr ft^) for LH 2 . 
Required values of Fa were found from 


<5/A incident 

F = ;*-T- (3-22) 

^'■^wicking capability 

For the LO 2 standpipe the worst case F^ was found to be a maximum of 0, 0076. For 
the LH 2 standpipe the worst case F,j was found to be a maximum of 0. 058. Both these 
values are weh below the majority of F^j's found in Table 3-6. 

The critical height at which breakdown could occur was evaluated by determining the 
maximum vap )r volume that could be generated by incident heating to the start basket. 
Calci-latlons were performed for the period between the first and second bums of the 
two bum S 3 mchronous equatorial mission. Maximum liquid volume expended was 0. 15 
m*^ (5.37 ft^) for LH 2 and 9.36 x lO""! m^ (0.026 ft^) for LO 2 . The vapor volume was 
converted into a liquid level below the top of the standpipe. Hjcr ’^'tst be at least as 
great as this h^ig it. For LH 2 , Hicr^as 13.5 cm (5.3 in) and Hicr/HSmax was 
0. 0022 where H0|^^ is the height of liquid that can be supported by surface tension. 
Similarly, for LO 2 , Hicr was 7.32 cm (2.88 in) and was 0.00149. Using 
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Equation 3-20 values of Ku were determined to be a minimum of 0. 92 for LH 2 and 
0. 991 for LO 2 . In order for the HiCr/HSmax requirements to be satisfied, at the Fj 
values given above, Ku values must be less than these quantities. Only one value out or 
30 rims had a K^ greater than that required to make the window screen operate properly. 

Based on these calculations, the window screens should be capable of remaining wet 
during the complete set of mission conditions. The main screen multiple barriers 
should perform successfully based on the testing reported in Reference 3-2. 

3. 3 . 2 CONCLUSIONS. The data presented in Table 3-6 can be used as a rough estim^e 
to find the limits of the thermal conditioning capability of a given configuration. For 
detailed thermal conditioning design data, additional experimental work will be required. 
This woidi is described in Section 3. 4. 

3. 4 VAPOR INFLOW TESTING RECOMMENDATIONS 

As described in previous sections, difficulties in obtaining repeatable data were 
exfierienced during the vapor inflow testing. Because of funding limitations many 
of these difficulties could not be overcome within the scope of this contract. The 
following pages describe the recommended experimental program to be implemented 
in order to resolve the operating characteristics of screen window configurations 
required for allowing vapor to pass into a cryogenic capillary device subjected to 
incident heating. The specific problems uncovered during testing are listed in Table 
3-7 along with their proposed soluticm. 

The general plan recommended is to separately evaluate wicking of screens fed by 
overlapping screens and vapor flow across wetted screens. Following these tests an 
overall test of the window concept is recon.mended, with an objective similar to that 
of the testing described in Section 3.2. All tests should be conducted with a controlled 
environment, minimized handling during filling and testing, and automatic recording 
of data. 

The basic chamber to be used is shown schematically in Figure 3-20. This will be 
used for all three types of tests: wicking across an overlap joint, and v^or flow across 
wetted screens with both small specimens and with a simulated capillary device and 
configuration. The chamber incorporates the recommendations of Table 3-7. Visual 
data will only be required during wicking tests. 

Wicking across an overlap screen joint will be studied first. Fabrication of overlap 
joints that produce high wicking rates that are repeatable from test to test will be 
studied. Primary configurations will be spot welded screen sandwiches. Use of 
spacers and perforated plates will also be considered. The overlap joints will be 
tested using techniques developed in References 3-2 and 3-3. The wick will be 
stretched across knife edges with the innermost section of the overlap immersed in 
the test fluid in the chamber shown in Figure 3-20. For each specimen, wicking 
results with the overlap joint will be compared to that of screen alone in order to obtain 
a measure of the fraction of wicking that is still available with the overlr.n .wint. 
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Table 3-7* Vapor Inflow Test Improvements 


Problem 


Solution 


1. Basic evaluation of the apeclllc mechan- 
isms involved in aucceaaful operation 
of the window configuration have not been 
separately evaluated. 


I, Analyse each specific mechanism (e. g. , screen wicking from overlap 
screen )oints, and vapor flow across a wetted barrier) separately by 
develq;>lng a model for each mechanism and testing it separately after 
determining the individual riialytical relationships 


2. Apparatus is moved during each test in 2. Remotely fill the test article that is mounted in the test chamber. Use 
order to fill and mount the test article. valves and lines to control filling, venting, etc. Use an enclosed test 

Filling and mounting are time consuming, chamber to control environment, 
nonrepeatable due to the uncontrolled 
environment. Operations during filling 
are hazardous. 


3. Limited measurements were made due 
to time restrictions and paucity of 
instrumentation. 


3. Monitor and record all significant variables to determine the significance 
of all events that cculd have an effect on performance. Monitor temperatures 
in aU areas where heating Is occurring, (particularly in the heater area). 
Monitor pressures in all significant areas. Pressure measurements will be 
sensitive enough to detect pressure changes of 0. 05" of test fluid head. 
Monitor liquid ^evels (of window and main screen volumes) versus time. 


4. Heat flux to the screen was difficult to 
quantify. 


5. Net heat flux (or evaporation) v^ies as 
liquid is boiled off and saturates 
environment. 


4. Apply heat flux directly using guarded heaters. Perhaps use screen as 
the heater element or use an overlay type beater with holes for venting 
and insulation on the top and side, protecting the box from heating away 
from the screen and concenlrating all the heat to the screen. 

5. Maintain atmosphere partial pressure by venting the container enclosing 
the apparatus over a narrow pressure band (perhaps use the baratron cr 
barastat). 


6. Temperature variations will cause 
changes in property values that 
could a^ect results. 


6. Control the temperature of the apparatus by using thermocouples and 
and heaters or coolers (air conditioning). Run tests in an air conditioned 
room. 


7. Cleanliness of apparatus and measuring 
tools was not carefully maintained. 

8. Water condenses on the screen suiiace 
from the atmosphere during handling. 
No control was maintained on the 
humidity (HnO) of tl.** atmosphere 
around the test specimen. 

9. Test fluid purity is not controlled and 
could contribute to poor repeatability 
of results. 

10. Cont.'ol iLe geometry more precisely 
than can be done with glue joints 
betv'een screens and siq^port. 

11. Control the wicking from the window 
screen to the main screen transverse 
to the desired wicking path. 

12. Control spacing between the screens in 
the window' /main sc teen overlap area. 


7. Clean the screen and other portions of the apparatus and the Instrumentation 
carefully before the tests and preceding each test after the apparatus has 
been disassembled. 

8. Control the atmosphere using boiled off bubbled through test fluid 
through fittings Into the transparent test enclosure. 


9. Use clean test fluid on an open loop basis or use filtration and purification 
processes to alloa* reuse of fluid. (If necessary, use reagent grade fluid). 

10. Consider an all welded configuration for controlling spacing, minimi zlrg 
geometric variations and preventing leakage. (If a transparent apparatus 
is desired some glue joints will be required between the supports and 
the transparent material). 

11. Use DODwettlng Btr*ps such as Teflon to prevent wicking from other than 
the desired direction. 

12. Use spot welding to join materials. Conduct separate fabrication evaluation 
to determine the optimum welding pattern and specimen overlap 
configuration. 
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TO KESEItVOIK 


Figure 3-20. Schematic of Basic Test Chamber 

Vapor flow across wetted screens would be tested in the basic chamber with the 
apparatus shown in Figure 3-21. The circular screen is wetted from an annular 
reservoir. Pressure will be cc .’trolled Independently on both sides of the screen. 

A typical test run would reduce pressure in one of the chambers until vapor broke 
through the screen from the higher pressure chamber. The mechanics of rewettlng 
the screen would be studied by manipulating the pressures in each chamber. Many 
screens should be tested with at least two test fluids. The apparatus would function 
as follows. The bottom box is the low pressure side of the apparatus and is held at 
constant pressure with a relief valve. The high pressure side would have a regulator 
that could be adjusted to regulate to a pressure above the bottom box pressure. The 


ADJUSTABLE 
REGULA TOR 



Figure 3-21. Schematic of Test Article lor Vapor Flow Across Wetted Screens 
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pressure In the upper chamber would be raised to exceed the bubble point across 
the screen. Flow across the screen would be noted by flow across the relief valve. 
Pressure in both chambers will be careMly measured. Pressure will then be reduced 
to a low€ r level until the screen reseals and flow is stopped. Cycles will be repeated 
to assure reliability. Other variables will be steady state flow rate across the screen 
and time at this flow rate. Analysis should be performed, prior to testing, to predict 
the wetting capability of the screen subjected to vapor flow. This analysis will be used 
to correlate tost results. 

Combine fabrication and test information to design an apparatus to test vapor flow 
across a window of coarser screen within a main screen surface. The conl^guration 
may be similar to that used in the test articles de'^-oiibed in Section 3. 2. Changes 
that ^pear possible will be barriers to wicking between the window and main screen 
transverse to the wicking path being measured. The use of larger specimens more 
closely ^proximating a real configuration is made possible by the automatic recording 
of data making the resulting long test luns more feasible. Tests would be run with several 
screens and test fluids. An analytical model, combining the basic tiieory developed in 
Section 3.1 with the correlated models for screen wicking with an overly joint and vapor 
flow across a wetted screen, will be compared to ttie test data. The result of the total 
test program will be information that can be used for design of capillary devices for 
propulsion stages. 



4 

FEED SYSTEM COMPARISONS 


Feed system alternatives were developed for pressure fed engines and combinations of 
feed system components in order to determine the optimum propellant feed system for 
the Centaur D-lS. Comparisons were made between prc^ellant settling and capillary 
acquisition, thermal subcooling and pressux'lzation for boost pump NPSP (net positive 
suction pressure), boost pumps, thermal subcooling and pressurization for turbopump 
NPSP, uncooled and cooled propellant ducts, and pumping coolant back into the tank or 
dumping coolant overboard. Capillary device designs used for the Centaur D-IS 
reflect the refilling and vapor Inflow analysis described In Sections 2 and 3, 

Comparisons were made on the basis of payload penalty, hardware weight, reliability, 
electrical power consumption and mission profile flexibility. The comparisons were 
made for three engine candidates: (1) the existing RLlOA-3-3, and two lower NPSP 
alternatives, (2) the RL10A-3-3A and (3) the RLIO Category 1. Characteristics for 
these two advanced engines are described in Reference 4-1. Three missions were 
considered for each engine candidate: (1) a one-bum planetary mission, ( 2 ) a two^ium 
synchronous equatorial mission, and (3) a five-bum low earth orbit mission (Tables 
1-1, 1-2 and 1-3, respectively). A total of ten feed system concepts were compared 
for each of the three engines and three missions. 

4. 1 SELECTION OF FEED SYSTEM ALTERNATIVES 

In order to select the most promising feed system alternatives for comparison, a matrix 
of system candidates was constructed by selecting an alternative from each of the 
following subsystems or operations. 

A. Acquisition 1. 

2 . 

B. Boost Pump NPSP 1. 

2 . 

3. 

C. Turbopump NPSP 1. 

2 . 

3. 

D. Propellant Duct Cooling 1. 

2 . 


Propellant Settling 
Capillf iy Device 

Boost Pump Pressurization 
Boost Pump Thermal Subcooling 
No Boost Pump 

Boost Pump 
Thermal &ibcooler 
Pressurization 

Uncooled 

Cooled 



E, Coolant Handling 




1. No Coolant Required 

2. Coolant Pumped Back Into the Tank 

3. Coolant Dumped Overboard 

Selecting all the possible combinations produced 108 alternatives. Infeasible combina- 
tions were then identified. For example, boost pumps cannot be used with cooled 
propellant ducts because boost pump cooling is impractical. Pressurization and 
capillary acquisition were shown to be infeasible in NASS- 17802. Other incompatibilities 
are turbopump subcooler or pressurization with a boost pump; no coolant required and 
a cooled feedline, turbopump thermal subcooler, or boost pump thermal subcooler; 
turbopump NPSP from a boost pump and no boost pump; coolant pumped or coolant 
dumped with uncooled duct and a pressure fed system; no cooling required with coolant 
pumped or dumped; and coolant pumped or dumped with no subcoolers or cooled duct. 
Candidates using cooled propellant ducts and not having capillary devices are also 
infeasible because liquid could not positively be contained within the duc<- or supplied 
for cooling the propellant duct without a capillary device. 

The screening process resulted in twelve combinations that were feasible. These are 
listed in Table 4-1. Candidates selected for additional analyses are circled. The reasons 
for selecting these ten concepts are described in the following paragraphs. Comparisons 
performed on previous studies (NAS3-17802 and NAS3- 19693) have shown that boost pump 
thermal subcoolers have weight advantages compared to boost pump pressurization, 
pumping cooling fluid back into the tank is advantageous over dumping fluid overboard, 
boost pumps are lighter than turbopump thermal subcooling and uncooled propellant ducts 
are lower in weight penalty than cooled propellant ducts. These comparisons were used 
as guidelines in seleccing the feed system alternatives to consider. Several key 
comparisons were desirable in order to determine the best feed system candidate. 

A critical comparison to be made was between capillary acquisition and settling. Several 
pairs of concepts in Table 4-1 could have been compared for this purpose (i.e. , B and K, 
C and L, D and M, or E and N). Based cm comparisons done on previous contracts, 
concepts B and K should be the lowest weight of all these pairs and, therefore, were 
selected as the primary comparison pair for propellant settling versus capillary 
acquisition. 

Concepts A and B were selected for comparing pressure fed boost pumps to thermally 
subcc'oled boost pumps. 

Concepts B, D and H compared the use of a booLt pump, thermal subcooler, or pressure 
feed for supplying turbopump NPSP. 

Another comparison involved the propellant duct cooling. Uncooled propellant ducts, 
cooled propellant ducts with cooling fluid dumped overboard and cooled propellant 
ducts with cooling fluid pumped back into the tank were compared by considering 
concepts O, N and P. 
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Table 4-1. Feed System Candidates 




Concei 

5t 




m 

c 


(D 

(D 


2 ) 

M 




A. 

Acquisition 

1 


1 

1 

1 









1. Propellant Settling 

2. Capillary Device 

1 

X 

1 

1 

1 

X 

X 

X 

X 

X 

X 

X 

B. 

Boost Pump NPSP 



1 











1. Boost Pump Pressurization 

2. Boost Pump Thermal 
Subcooling 

3. No Boost Pump 

1 

X 


1 

1 

X 

X 

X 

X 

X 

X 

X 

C. 

Turbopump NPSP 














1. Turbopump NPSP - Boost 
Pump 

2. Turbc^ump NPSP - Thermal 
Subcooler 

3. Turbopump NPSP - 
Pressurization 

X 

X 


X 

X 

X 

X 

X 

X 

i 

1 

X 

X 

X 

D. 

Propellant Duct Cooling 














1. Uncooled Propellant Duet 

2. Cooled Propellant Duct 

X 

X 

X 

X 

X 

X 

X 

1 

X 

X 

X 

1 

X ' 

X 

E. 

Coolant Handling 








1 

1 


1 



1 

i 

1 

1. No Coolant Required 

2. Coolant Pumped Back Into 
The Tank 

3. Coolant Dumped Overboard 

X 

X 

X 

X 

X 

X 

X 

1 

X 

X 





Comparison of subcooling coolant flow being dumped or pumped could have been made 
for boost pump subcooling by comparing concepts B and C or concepts K and L. 
Concepts K and L were evaluated for this purpose since most of the data for this 
comparison had already been generated on NASS- 17 802. 

For comparing pumping or dumping subcooler flow that provides turbopump NPSP, 
concepts D and E and concepts O and P were compared. 

Based on the considerations described in the preceding paragraphs, analysis of concepts 
A, B, D, E, H, K, L, N, O and P will be required to provide the necessary 
comparisons. (This is the same number of concepts originally proposed for analysis. ) 
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Candidates that were analyzed in Reference 4-2, NAS3-17802, and Reference 4-3, 
NAS3-19693, have their letters underlined. 

Each of the feed system candidates circled in Table 4-1 was analyzed for the ono'burn, 
two-bum and five-bum missions described in Reference 4-2. The matrix of engine 
alternatives analyzed for each candidate feed system is shown in Table 4-2. Concepts 
using boost pumps were only examined for the existing RLlOA-3-3 engine since new 
boost pump designs were required for the other two engines. 


Table 4-2. Matrix of Feed System Concepts and Engine Alternatives 
Compared 


Feed System Concept 
(See Table 4-1) 

Alternatives 

RLlOA-3-3 

RL10A-3-3A 

RLIO Category 1 

A 

X 



B 

X 



D 

X 

X 

X 

E 

X 

X 

X 

H 


X 

X 

K 

X 



L 

X 



N 

X 

X 

V 

0 

X 

X 

X 

P 

X 

X 

X 


4.2 PAYLOAD PENALTY 

Payload penalty calculations were performed using payload sensitivity factors given in 
Table 4-3. These factors were applied to the differential proi -llant and Jettison weights 
for each of the ten concepts. Information generated for the systems comparisons 
presented in References 4-2 and 4-3 were used as a foundation for conducting this study. 

Capillary device weight calculations were performed to incorporate changes in capillary 
device volume and comer construction. The increased capillary device volume was 
associated with concepts N, O and P. These concepts considered turbopump thermal 
subcoolers with greater volume than the boost pump thermal subcoolers previously 
used in the weight calculations of Reference 4-2. Comer construction techniques that 
would aid in wlcklng and ease fabrication requirements were incorporated into capillary 
device designs. These changes both resulted in increased capillary device weiv’"’’t. 
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Table 4-3. Centaur D-lS Payload Sensitivity Factors 


Criteria 

Mission 

Planetary 

Sync. 

Equatorial 

Low 

Altitude 

Jettison Weight 

-1. 000 

-1. 000 

-1. 000 

Propellant Weight 

+0. 646 

+0. 587 

+0. 420 

LH 2 and LO 2 Loss Before Bum No. 1 

-0. 646 

-0. 587 

-0.420 

(Without Igp Effeet) No. 2 


-0. 946 

-0.610 

No. 3 



-0. 771 

No. 4 



-1.215 

No. 5 



-1.337 

LH 2 and LO 2 Lost After Last Bum 

-1.646 

-1.587 

-1.419 

(Residual or FPR) 




Auxiliary Propellant Used No. 1 

0 

0 

0 

Prior to Bum No. 2 


-0. 459 

-0. 072 

No. 3 i 



-0. 352 

No. 4 i 



-0. 796 

No. 5 



-0. 917 

Auxiliary Propellant Used After Last 

-1. 000 

-1. 000 

-1. 000 

1 Bum (Residual) 





Pressurization system weights for all systems were computed using a detailed inventory 
of Centaur pressurization system weights. These computations revetied that the 
original pressurization system weights documented in Table 2-6, Reference 4-2, 
included excessive vent system weight for the baseline system and insufficient support 
system weight for the system not requiring main tank pressurization. This is a 
significant revision, resulting in a weight increase of 13. 5 kg (29. 8 lbs) for the 
pressurization system used with thermal subcoulers and a reduction in weight of 
30.3 kg (66. 8 lbs) for the baseline system (Concept A). 

Lower tanking density ^ malties were computed based on increased tank pressure 
requirements for "blowdown" options not using boost pumps and using the RLlOA-3-3 
or RL10A-3-3A engines. Pressure fed option with -3 engine also requires higdiui^ tank 
pressure. 

Propellant supply duct weights were computed using w eights taken from Reference 4-3. 
Recirculation weights, included in sipply duct weights in Reference 4-3 are separated 
out into "other hardware" in this study. 
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Supply duct chilldown weights reflect the use of the recirculation system to chllldown 
the propellant ducts when boost pumps are used. Cooled ducts do not have any duct 
chilldown requirements. CXher chilldown weight data was taken from Reference 4-3. 

Supply duct cooling fluid dumped overboard was changed considerably from Reference 
4-3. The cooling of the sump area was neglected in Reference 4-3. Calculations were 
performed to include the sump area. 

Hardware for keeping supply ducts wet was taken directly from Reference 4-3, Table 

3- 12 for all concepts using cooled feedlines. Excess LO2 in the duct over the amount 
required for chilldown of these options and engine chilldown requirements were also 
taken directly from Reference 4-3, Table 3-12. 

Subcooler weights were calculated based on schematics described in Reference 4-3. 
Subrooler cooling fluid dumped overboard was taken directly from Reference 4-3, 

Table 3-12. Boost pump weights, peroxide usage and settling system penalty were 
taken directly from Reference 4-3, Table 3-10. 

Pumping system weights to return coolant to the tank were computed for feedline 
coolant, thermal subcooler coolant (turbopump or boost pump NPSP) and combinations 
of both. System weight penalties included boiloff, battery and hardware (sump, pump, 
’Ines and and accumulator) weights. 

Sump assembly weights are a function of whether a boost pump is used. Sump assembly 
weights with and without a boost pump are taken from Reference 4-3. The volume penalty 
of the added hardware is calculated as the volume of propellant which cannot be tanked 
because of this additional tank hardware. Calculations included capillary device, sub- 
cooler, pressure bottle, and thrust barrel revisions. The volume penalty of increased 
tank skin weight for concepts requiring tank pressures in excess of the baseline levels 
IS taken from Reference 4-3. Pumping options are assumed to have hardware outside 
the tank. 

Subcooler cold side residuals were computed based on cold side volume and appropriate 
payload partials. Thrust barrel revisions required to aid in refilling the capillary 
devices were determined from Reference 4-2. Other hardware weight Included in the 
payload weight comparison is that of the recirculation system for options containing 
boost pumps and dump valves, propellant utilization probe revisions and bypass lines 
for options containing capillary devices. Weight data was obtained from Reference 

4- 2. 
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All payload penalties are Included in Table 4-4 through 4-12 for the three engines and 
three missions of interest. 


Tabulations were made of the feed system alternatives to illustrate the relative merits 
of the subsystem alternatives. These comparisons are shown in Tables 4-13 through 


4-21. 


In Tables 4-13 and 4-14 a capillary acquisition system is compared to a propulsive settling 
system. Capillary devices are heavier than propulsive settling for <he missions consid- 
ered. For missions with greater than five bums, capillary device payload weight penalties 
will increase at a lower rate than propulsive settling payload weight penalties. 

Table 4-1.^ compares boost pump thermal subcooling to pressure fed boost pumps. It 
appears, for the cases considered, that thermal subcooling is lower in weight penalty 
for missions having more than two bums. 

Table 4-16 compares the three methods of supplying turbopump NPSP. It appears that 
using a boost pump to supply turbopump NPSP is the best option. 

Table 4-17 compares propellant duct cooling options. The best option is to not cool the 
duct and to pump any other coolant back into tne tank. < x'his option, concept M, was 
assessed by adding the weight of the pumping system of concept O to concept N and 
deleting the coolant los. from concept N). After this, the next best option is to cool 
the duct and pump the coolant back into the tank. The next best option is not to cool 
the duct. The worst, from a payload penalty standpoint, is to cool the duct and to 
dump the coolant overboard. 

Tables 4-18, 4-19 and 4-20 compare pemplng coolant back into the tank versus dumping 
the c<x)lant overboard. In all cases the coolar.*^ pumping system is lower in payload 
weight penalty than dumping the coolant overboard. 

Pressuri itlon system options are compared in Taole 4-21. Options considered v/ere 
cryostored or ambient stored helium for an autogeneously pressurized LH 2 tank and 
cryostored gaseous heliun> for the entire LH 2 tank pressurization requirement. 

Aulogeneous LH 2 tank pressurization with cryogenic helium storage is the lowest 
weight alternative. Ambient storage of Glle with aulogeneous ldl2 tank pressurization 
is heavier than cryogenic storage. The option with the highest payload penalty is the 
one using GHe pressurization and cryogenic storage. (Because of (he added complexity 
of cryogenic storage, ambiently stored systems are baselined for the pressure fed 
syftem, II, wherever possible.) 

Summarizing the payload weight comparisons, the following conclusions are true. 

Settling is superior to capillary devices for the missions considered. Subcooling of 
boost pumps by thermal or pressure means is a close trade-off, with thermal 
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Table 4-4. Payload Weight Penalties for System Comparisons* kg (lbm)» One Bum 
Mission, RLlOA-3-3 Engine 
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Table 4-S. Payload Weight Penalties for System Comparisons, kg (Ibm), Two Bum 
Mission, RLlOA-3-3 Engine 


rSKD 

smkiu 

coNC^rr 


weiotrr 

PENALTY 

LLCMCNT 


,! li 

ii 


Ml {,■ 


iil SiJ 

51? 5|| 

ih 2 ia 
\V II* 


il lil ill 111 III 111 111 Hi 


I 1 ll 

!i ill 


Caplllari Otvtea 

Lowrar Taaktam 
DapBlly 

PrBBBurltBtloa 

Syst«iii 

St^iply 

Duel 

Supply Duct CbltldOifB 
PropcUai4 

Supply Duct CiKillog 
Dumped Overboard 

Hardware to Keep 
Duels Wat 

Exoast LOg In Ltoaa 
Ovar Cbllldown 

Engine CbtUtknm 
Propellant 


IS. BooM Pump 

15, H^^j Hsaga, Uooat 
Pump 

14, Sattitng Syateio 

Including Panalty 

16, Sump Assembly 


Pumping Sys. to Hatum Ltl^ 
Coolant to Tank tnol. 

BolloH, Batt.Alldirr. 


Ul| 

LO, 

. 

- 

“ 

- 

Lll, 

- 

- 

69 (191) 

09 (161) 

W>S 

- 

- 

438 (946) 

499 (906) 


jioo (241) 

49 (94) 

43 (H) 

43 (94) 

Lllg 

10 (S3) 

10 (23) 

10 (23) 

10 (23) 

LO^ 

10 (S3) 

10 (23) 

10 (23) 

10 (23) 

til. 

- 

- 

3 (7) 

9 0) ! 

1.0a 

- 

- 

6 (10) 

6 (10) 

Ulg 

- 

- 

- 

- 


“ 

- 

- 

** 

i.iij 

- 

- 

- 

- 

LOa 


- 

- 



I 

- 

• 

- 

uu 

1) (26) 

~i“^6T 

U (25) 

IT'cm)" 


2 (6) 

2 (6) 

2 (6) 

9 (6) 

uia 

- 

9 (19) 

49 (107) 

49 (107) 

LOa 

- 

11 (24) 

20 (43) 

SO (43) 

"lh7 

- 

__ 

- 

177(3947 

,4 

- 

- 

- 

100(1261) 

Lila 

99 (05) 

39 (86) 

- 

- 


29 (04) 

29 (04) 

“ 

- 


17, Voluma Panally Due 

Lila 

to Hardware Aiklad 

loa 

10. Fluid Healduals 

Liu 

Cold b.Ja SubcooUr 

J^2_ 

19. Thruat Darral 


Ravlaloot 


20, Otitar llantwara 

21. Tank Skin A Walglit 


Ovar baaallna 


TarAidi 


fi 04) e (14) 

46 (lOS) 49 (109) 49 (109) 49 (109) 

4 (14) 6 ()4) 3 (7) S (7) 

8 07) 9 (17) 7 (16) 7 (16) 

0 (13) 10 (S2) 

6 (12) 11 (26) 


1 ( 2 ) 1 ( 3 ) 

4 (10) 4 (10) 


4 (8) 21 (47) SI (47) 

39 (89) 61 (112) ei (112) 




82 (837) 292 (646) 929(2048) 1666(3647 



86 (187) 86 (187) 102 (226) lOj (226) 

SO (461 20 (46) 21 (17) 21 (47) 

69 (161) 09 061) 
436 (906) 438 (966) 

43 (94) 43 (94) 43 (94) 43 (94) 

10 (23) 10 (23) "l~(b) 10 (S3) 

10 (23) 10 (23) 10 (23) 10 (23) 


102 (226) 
SI (47) 


80 (178) 90 (170) 


11 

(26) 

2 

(6) 

t 

(19) 

11 

(24) 

- 

99 

(06)1 

29 

(04) 

6 

(14) 

- 

8 

(14) 

i • 

(17) 

8 

(13) 

6 

(12) 



1 

0) 

9 

()») 

4 

(0) 

99 

(03) 

6 

(11) 

0 

(10) 


13 (29)> 179 (394) 

40 (88)1668(1261) 


0 *)\ 8 iV 3 (tT] 3 (7) 



• ( 10 ) • ( 10 ) 6 ( 11 ) 6 01 ) 6 ( 11 ) 


905 (000)1407 (809)|l747(S848|m6(24«l|l999(4107 


3 

(6) 

11 

(23) 

21 

(47) 

61 

(118) 

6 

(111 

6 

(111 

7F 

(1601 

40 

(091 








































Figure 4-6. Payload Weight Penalties for System Comparisonst kg (Ibm)t Five Burn 
Mission, BLlOA-S-3 Engine 
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Table 4-7. Payload Weight Penalties for System Comparison, kg (Ibjn), One Bum 
Mission. RL10A-3-3A 
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Table 4-8. Payload Weight Penalties for System Comparison, kg (Ibm), Two Bum 
Mission, RL10A-3-3A Engine 
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Figure 4-9. Payload Weight Penalties for System Comparisons, kg (Ibm ), Five 
Burn Mission, RL10A-3-3A Engine 
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Table 4-10. Payload Weight Penalties for ^stem Comparisonst One Bum Mission, 
HLIO Category I Engine 
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Table 4-11. Payload Wel^t Penalties for System Comparisons, kg (Ibm) 
Two Bum Mission, RLIO Category I Engine 
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Table 4-12. 
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Payload Weight Penalties for System Comparisons » kg (Ibj^) Five Bum 
Mission» BLIO Category I Engine 
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Tabic 4-13. Settling Versus Capillary 
Device Payload Penalty Comparison 
(Systems B & K) 
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Table 4-14. Settling Versus Capillary 
Device Payload Penalty Comparison 
(Systems E & N) 
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Table 4-15. Boost Pump NSPS, 
Pressure Fed Vs Thermal 
Subcooling Payload Penalty 
Comparison, (Systems A & B) 


Table 4-16. Methods for Supplying 
Turbopump NPSP Payload Penalty 
Comparison (Systems B, D & H) 
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Table 4-17. Propellant Duot Cooling Options Payload 

Penalty Comparison (Systems N, O, P & Mk 
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Table 4-18. Boost Pump Suboooler 

Coolant Disposal Payload 
Penalty Comparistm 
(Systems K & L) 


Table 4-19. Turbopump Thermal Subcooler 
Coolant Disposal Payload 
Penalty Comparison Settling 
System (Systems O & E) 
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Table 4-20. Turbopump Thermal 
Subcooling Disposal Payload 
Penalty Comparison Capillary 
Dovices Systems O and P) 


Table 4*21 . Presaurisation System 
Options Payload Penalty Compar '^on 
(Systems H| 6 H 2 ) 
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subcooling preferred for the multibum missions. Pumping coolant back into the tank is 
better than dumping overboard. Unoooled propellant duets are better than cooled 
propellant ducts and boost pumps are the best method of siq>plving turbopump NPSP. 
based on these comparisons of relative payload weight penalties, feed system concepts 
A or B appear best for the Centaur D>1S application. Feed system concept A is the 
baseline system. Feed system concept B Is the baseline system with the substitution of 
thermal subcooling for main tank pressurization. 

4. 3 HA RDWARE WEIGHT COMPARISONS 

Hardware weight penalties are shown in Table 4-22 tbrou^ 4-30. Hardware weight 
comparisons are shown in Tables 4-31 through 4-39. Hardware weight comparisons are 
not as meaningful as the payload weight penalty comparisons shown in Section 4.1 
because dicy do not reflect differences in fluid weight penalty. 

Table 4-31 and 4-32 compare settling to capillary devices. Capillary devices are 
generally about 80 to 90 kg heavier than settling in terms of hardware weight for one, two 
and five burn missions. 

Table 4-33 shows that thermally subcooled boost pumps have lower hardware weight than 
pressure fed boost pumps. Table 4-34 shows that boost pumps have the lightest hardware 
weight of the turbopump NPSP supply methods (boost pump, pressure, thermal subcoolor). 

If pressurant is cryogonically stored then concept Hi is comparable in weight to boost pumps. 
For the Category I engine thermal subcooling is lighter than pressurization for suppl}ring 
turbopump NPSP while for the RL10A-3-3A engine, pressurization is lighter than thermal 
subcooling for supplying lurbopump NPSP. 
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Table 4-22. Hardware Weight PenaUies for %stem Comparisons, kg (Ibm). One Bum Mission, RLlOA-3-3 Engine 






p*d«iaim|«09 *)oaap«|M3 

a, 

it 

1 

IP 

ia 

o o 

p* m 

ai 

p p 

|a 

58 

1 1 

1 

A A 

tg 

m »• 

1 1 

pm 

pm 

P 

A 

fi4 

P 

SI X 

T 

i 

ptdanwi 901(000 
•Iona pnooo ‘lonoooqiit 0 
4nndoqjnx *90\A9q ^jvuidfo 

it 

1 

p* 

•P 

fa 

O 9 
m m 

la 

p p 

58 

1 1 

t 

£§ 
« «• 

la 

C 

mm 

P 

s* 

p 

r§ • 

a ^ 

S 4 

1 

8 

podtana 901(000 
*pna poiooooa *donoooqiiB X 
dwndoqmi '•o\jl9q ^i^udio 

ft 

8 » 
m 

1 

w 

§§ 

o o 

m* 

1 

S5 

58 

1 


i 

1 

P 

P 

8 8 
0-^ 

S X 

1 

i 

9011000 pod«a *90oa poiooooa 

*dama 9«oos •Sinooom j 
*00(000 ^»Uldio 

82 

1 

•la 

§a 

m S 



ii 

88 


i 

■ 

P 

s 

p 



Hi(ooo podvn^ *9000 
po|oooin *dt8nd 9000 Q ‘>i(ioooqiig iic 
(viuoqx *oo(Aoa iCiinKlto 

82 

1 

a§ 

e e 

m mm 


A A 

c5 

p 

ii 

p p 
« N 


P4 pm 

mp 

P P 

« eJ 

p p 

pm 

P 

p 

p 


p 

X 

8 

M 

pojpiboii 900(000 OM 

•p«a potoooon X 

*oqOQ 3 poj ojntiojici 

■ 

■ 

i 

5 

■■ 



■ 

9 

■ 

i 

1 


1 

podona 910(000 ‘(ona po(oosoo 
*«onoooqns domdoqjox *»in»OB 

1 

1 

’e* 

§a 

P o 

pp m 

1 

s$ 

pm - 

p e 

1 

§ 

8 

S 5 

« !• 

1 

1 

■ 

|s 

8 at 

poiinbod 910(000 
*9300 poiooooa O 
•luuoooqns diandoqjni *^ 1 | 9 los 

■ 

1 

♦ 

1 

1 

£ « 
e ^ 

p e 

1 

8 

9 

§a 

S £ 

■ 

■ 

e» p 

C 

P u 
p n 

9 

g 

X 

ei 

podton^ 910(000 
1 *pnn po(oooif) *dtim.{ (ooog 

•»u»(oooqns (ouuoqx '^OUIIIS 

1 

i 

1 

■ 

P •» 
— W 

P •- 

JP- * 

i 1 
88 


i 


■ 

£ 

n 

i ( 

s 

mp 

P 

N 

poitobOH 9i«iooo OM 
•vnq p»(003un *dum.; looog < 
*uonc*(40tOOJ,| *9(ut(9906 

1 1 

s 

X 

aa 

O P 

< 1 

1 t 

8 S 

* i 

P P 

rs ?i 

1- 

e 

w e- 

M PP 
mp w 

P P 

1 1 

1 

r» 

n 

I t 

8 

• 1 
8 

- y 
£ 

s > S 

^ ^ h fc 

§ S 

li 

u u 

5 >• 

is 

?§* 

a 

ff 

J 

1 

•i 

i 

B 

ii 

1 

li 

ff 

f 

1 

1 

1 

K* 

iS 

f 

< 

r 

aTi* 

1 

li 

m 

li 

1 

i 

< 

9 

1 


4-20 


^ e 










































Fable 4-23. Hardware Weight Penalties for Svstem Comparisons, kg (Ibm)* Two Bum Missi<m, 
RLlOA— 3*3 EnginG 






































































































































Table 4-25. Hardware Weight Penalties for System Comparisons, kg (Ibjn), One Burn 
Mission, RL10A-3-3A Engine 
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Table 4-26. Hardware Weight Penalties for System Comparisons, kg (Ibnj), 
Bum Mission, RL10A-3-3A Engine 
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Table 4-29. Hardware Weigbt Penalties for System Comparisons, kg <lbm)» Two Bum 
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Table 4-30. Hardware Weight Penalties for System Comparisons, kg (Ibm), Five Bum 
Mission, RLIO Category I En,^ioe 
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Table 4-31. Settling and Capil- 
lary Device Hardware Wei^ 
Compaxiscm (Systems B&K) 


EDftlM 
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PantUlMt kKlb^) 

Food Syotom 

Pood Syotom 
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Bonf) 
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S«Mlng 

Ciptuorr 

Dovioo 

HL10AO*3 

NoUt Botl 
pun 

1 
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5 

1 contfiii \ 
pedL 

1 

212 (468) 
212 (468) 
212 (468) 

Moatpumpa, 

i 

292 (64S> 
292 (649) 
292 (645) 

coolant 
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Table 4-32. Settling and Capillary 
Device Hardware Weight 
Comparison (Systems E 6 N) 
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Buma) 

Hardwar* Watglp 
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167 (STT) 

256 (562) 
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225 (502) 
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255 (627) 
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305 (673) 

3A 
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323 (712) 
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25S (862) 

342 (784). 

RUOA-S-3 

1 

266 (886) 

364 (808) 
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287 ...2) 

388 (848) 


$ 

821 (707) 

418 (824) 

Not / 1 ' 3opui&p flulicooUiig, coolant dtimpod. 


4-33. Boost Pump NPSP, Pressure 
Fed System and Thermal Sub- 
cooling System Hardware Weights 
Comparison (Systems A & B) 


Table 4-34. Turbc^ump NPSP, Boost 
Pump System; Pressure Fed System 
and Thermal Subcooling System 
Hardware Weight Comparison 
(Systems B« D & H) 


Eoidno 

Option 

Mlooton 
(No. of 
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Hardwar# W#lght 
Penattlaa , k| (Ibg,) 

Feed Syftom 

Feed Syatom 
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1 1 
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Note: Booet pump thennal oubcooUaf { coolant 

pumped back into the tank* 

i- .1 
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PMia}tt.a. kg (Ibg,) 
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Boost pump 
for Turbo- 
pump 
NPSP 

Thermal 
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for Turbo- 
pump NPSP 

Preofuriea> 
tlon for 
Turbopump 
NPSP 

RLlOCat 1 

1 1 


152 (399) 

186 (411) 
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155 (355) 

211 (466) 


5 


213 (465) 

208 (488)* 

RLlOA-3- 
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235 (517) 

515 (1135) 

3A 
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253 (555) 

483 (1065)* 
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273 (602) 

502 (1106)* 

RLlOA-3-3 
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212 (465) ' 

289 (630) 
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308 (679) 
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Propellaot duot cooling options M, N, O and P, shown in Table 4-35 « illustrate that 
the lowest weigM option, N, has an unoooled duot with other coolant dumped oveiboard. 
The highest hardware weight has cooled ducts with coolant pumped back into the tank. 
Concepts M and P are very close in hardware weight (unoooled duot with other coolant 
pumped versus cooled duot with coolant dumped). 

Table 4-36 shows lower hardware weight for dumping boost pump themal suboooler 
coolant versus pumping it back into the tank. Tables 4-37 and 4-38 show the same 
result for turbopump thermal subcooler coolant disposal. The weight of the coolant 
pumping system accounts for the added weight. 

Table 4-39 shows that oryogenioally stored pressurant with autogeneous steady state 
LHg pressurization is the lightest pressure fed system because of the lower weight 
cryo-bottles. (All increments are due to bottles and si^poit hardware.) 

4.4 RELATIVE RELIABILITY 

Relative reliability for each of the ten concepts identifled in Table 4-1 was determined 
by analyzing each major subsystem component to determine mean missions between 
failures (MMBF). Component features are similar to those presented in Reference 4-2, 
Appendix A. For pressure fed systems, not previously analyzed, component character- 
istics are as shown in Table 4-40 and Figure 4-1. 

Elements used in determining the relative reliability of each candidate concept are 
identified in Figures 4-2 to 4-11. 


Table 4-35. Propellant Duot Cooling Options, 
Hardware Weight Comparison (Systems 
N, O, P & M) 


Enidno 

OpMofi 

Mliiiiloii 
(No. of 
lluma) 

llardwtre Wetffht (Ibm) | 

Food Syatofn 

1 

1 

Feed Syetem 

Feed System 

N 

O 

P 

M 

Unoooled 

Dttct 

Coolftnt 

Dumped 

Tooled 

Duot 

Coolent 

l\itfiped 

Cooled 

Duot 

Coobuit 

Dumped 

Unoooled 

Duot 

Coolant 

Pumped 

RMOCtt 1 

1 

250 (502) 

207 (051) 

281 (010) 

271 (597) 


2 

228 (502) 

271 (505) 

253 (550) 

246 (541) 


S 

285 (027) 

328 (722) 

310(081) 

303 (Ml) 

IILlOA-3- 

1 

305 (073) 

943 (757) 

326 (718) 

333 (T12) 

3A 
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323 (712) 

902 (799) 

343 (767) 

342 (754) 


5 

342 (734) 

983 (845) 

362 (799) 

363 (800) 

RtlOA-n-H 

1 

304 (803) 

400 (898) 

385 (861) 

380 (850) 


2 

385 (840) 

429 (949) 

400 (897) 

409 (901) 
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419 (924) 

405 (1028) 

440 (972) 

446 (980) 

ipPPffP 

Hi 

eoollfiif t»pt|o 

mpip 

iplllarv deviec 

i» aed 
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Table 4-36. Boost Pump Suboooler 
Coolant Disposal Hardware 
Weight Comparison (Systems 
K & L) 


Rnplne 

Option 


llardurare Weifht 
renalUea , kg (lh,y|) 

Feed Syetem 

Feed Syetem 

K 

h 

(No. of 
liiirtia) 

i'f»o1ant 

ritmpotl 

Coolant 

Dumped 

RMOA 3 3 

Note: Cooli 
nuluv 

1 

292 (015) 

292 (045) 

292 (015) 

nl from boost 
plltarv devleej 

281 <020) 

281 (620) 

281 (020) 

pump 
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Table 4-37* Turbopump Thermal Sub- 
oooler Coolant Dlaposal Hardware 
Wei^ Comparison* Settling 
System (Systems D 6 E) 


Cnirtne 

Option 


Hsrdtrftre Wolgbl 
Peneltlee, kg iVb^) 

Feed System 

Feed Syetesi 

D 

£ 

(No* of 
Burse) 

Coolant 

Riamped 

Coolant 

Dumped 

RLlOCat 1 

RLIOA-S- 

3A 

RLl0A-3*3 

Note; Coob 
s\ibo< 

1 

1 

2 

5 

1 

2 

S 

1 ! 
2 
S 

mt dieposi 
^llne for 

Its (333) 
m (3S8) 
213 (468) 
238 (517) 
253 (558) 
273 (502) 
283(630) 1 
308 (679) 
343 (757) 1 

il flrom turbo| 

•eltUiur syete 

1 — 1 

l«7 (377) 
133 (307) 
198 (432) 
218 (481) 
238 (520) 
255 (582) 
288 (588) 
287 (832) 
321 (707) 

^ump tbennal 
me* 
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Table 4-38* Turbopump Thermal Sub- 
cooling Coolant Pumped Versus 
Dumped Hardware Wel^t 
(Systems 0 & P) 
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Mteeton 
(No* of 
Bumi) 

Hardvaro Waif 18 
PonaltioiR kg (tbg|) 

feed Syetam 

r«*d8nt«i 

0 
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Coolant 

Pumped 

Coolmt 

Dumped 

RLlOCtt 1 

1 

287 (661) 

281 (616) 


2 

271 (888) 

283 (568) 


5 

3M (712) 

310 (881) 

RLIOA-S- 

1 

343 (757) 

326 (718) 

3A 

2 

382 (729) 

343 (757) 


3 

363 (>845*) 

352 (799) 

RLIOA'3'3 

1 

406 ^898) 

385 (851) 


2 

429 (948) 1 

406 (597) 


8 

468 (1028) 1 

440 (972) 

Notet Coolant dlapoeal for turb<n>ump 

suboooU&f with oiu>UUry davtooo* 


Table 4-39. Pressurization System Cations Hardware 
Weight Comparisons (Systems H* Hi and H2) 


Engine 

Option 


1 

Har<>*ar* PanaltiM (kf db^> | 

Mieeioo 
(No. nf 
Bures) 

Feed Syetem 

Feed System 

Feed System 
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«2 

Autofer.oua 
U?: fenk 
Ambient 
Stored OHe 

Autogenoue 
LH 2 Tank 
Cryoetored 
GHe 

CM* for 
LH; Ttnk 
Cryoatorad 
GHa 

RLIOCU t 

1 

188 (411) 

175 (380) 



2 

211 (466) 

159 (374) 



* 

• 

227 (901) 



1 

1 

illilJIJI.llttliijlllliiJilLI.iiMI 

RLlOCat 1 

1 

109 (241) 

99 (209) 



2 

134 (2M) 

89 (197) 



$ 

« 

113 (249) 



* IniufflolMit room for ambient stored bottles* 


Table 4-40. Pressurization System Bottles for Pressure Fed Turbopumps 


MISSION 

Ul*10A-:i -:iA l%NCUNK 

IU,I« (.ATKtSOUY 1 KNUINJ.; 

oNfc: iu;uN 

3 SMAI.I. AMHIKNT UUITl.KS 

1 LAtU'iK AMmKNT UolTI.I'; 
► 1 SMALL AMIHIiN'r lutm.l-: 

'I'WO UIIICN 

1 «MAI.I. (;UY«Mik;NK:AU.Y STOUi:i) IU»TTI.K 
1 ] SMAI.I. AMDIL'N'I' liU'm.fc: 

1 lUJK A M HI KNT m » m .K 

12 SMALL AMHIKNT I«» I 'I'LL.S 

I'lVK UUItN 

2 l<AliUK (?HyotilONICAM.Y SI‘OIU;i) UOTTI.KS 
H SMAM. AMlMEN i' UOITI.I-: 

1 SMALL fUYtMU:NU;AI.LY .SrOlti:i» lio'm.L 
Il SMALL AMIJIKN r IIUI TI.K 


4-31 
















TablA 4-41. Mission Profile and Enviroiiinent Summary 




HELtUM PRESSURE BOTTLES ARE STORED W THE 
LIQUID HYDROGEN TANK. A HEAT EXCHANGER 
IS USED TO RAISE PRB88URANT TEMPERATURE TO 
USAGE LEVELS IN THE LHg TANK. COLO HELIUM 
IS USED FOR PRES8URIZINO THE LO 2 TANK. 

FOR SYSTEt* K” CRYOGENIC HELIUM IS USED FOR 
ALL MAW TANK PRESSURIZATKM. FOR SYSTEM 
H AND H' THE HYDROGEN TANK IS AUTO<fflNEOUSLY 
PRESSURIZED AFTER THE WITIAL BURP WITH 
HELIUM. 


Figure 4-1. Cryogenlcally Stored Helium Pressurant 
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EruP^ffij 






Fftilar* RatMs 

Bottto (HjOa, H«) 
Sattlinf Motor 
Sotonotd V«lv* 
OrUUM 


X 

0.6 X 10*« 

8.0x 10"« 
4. 0 * 10“< 
0,6x 10*« 


Chook Volvo 6.6 y 10-* 

Bubbtor« DUflpotor 1. 0 x 10*^ 

Booot Pump iS.Ox'O*^ 

Proonir* Traaaduoor 10 ^ i 


Figure 4*2. Sohematio and Failure Rates for Concept A (Settli 
Pres8urizatlon« Boost Pump» Uncooled Duct, No 
Coolant Required) 



Foliar* Rotooi X 


Throttlo Volvo 3. 5 X 10*^ 

Thonaol SttboooUr 1.8 x lO'O 

Surg* Took 0. 6 X 10“* 

Voouam Pump 9. 0 x lo**^ 


Figure 4-3. Schematic and Failure Rates for Concept B (Settling, 
Thermal Suboooling, Boost Pumps, Unoooled Duct, 
Coolant Pumped) 



Figure 4-4. Schematic for Concept D 
(Settling, Turbopump Suboooling, 
Uncooled Duct, Coolant Pumped) 


CoDOCpI 0 

(-) 

r Sorgo 


VAOOtttt 

PlUDO 


Figure 4-S. Schematic for Concept E 
(Settling, Turbopump, Subcooling, 
Unoooled Duct, Coolant Dumped) 
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RL 10A-3-3A ENGINE 
ONE BURN AMBIENT SYSTEM 


2 OF 3 WORKINO SAT. 



RL10A-3-3A ENGINE 

TWO BURN; ONE CYROGENIC 

<^E AMB. BOTTLE 

RL10A-3-3A ENGINE 

FIVE BURN 5 TWO CRYOGENIC *" DlMloitor 
ONE AMB. BOTTLE 

FAILURE RATE! 

HEAT EXCHANGER X -1. 8 x 10’^ 


LH, Prasi. 
Sraior 

LO 2 Pmii. 
" Senior 



1 LHj Prast. 
1 Smor _ 

LO 2 Prill. 
Senior 



LHj Praia. 
Senior 

^ LO 2 Preii. 
"" Senior 


2 OF 3 WORKING SAT. 


RLIO CATEGORY I ENGINE - SAME LOGIC SYSTEM AS ABOVE EXCEPT NO. BOTTLES 
ONE BURNS TWO AMBIENT BOTTLES 
TWO BURNS THREE AMBIENT BOTTLES 
FIVE BURNS ONE CRYOGENIC, ONE AMBIENT BOTTLE 


Figure 4-6* Schematic for Concept H (Settling, Pressure Fed Engine, Uncooled 
Duct, No Coolant Required) 
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Fsllure Hales: X 
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SUitBukat 1.4x10-^ 

Screened Cluuuiele 1.4 x 10*^ 

Sidicooler Heat Exchanger 1.8 x 10 *^ 

Orifice 0 . IS x l(r« 

Dump Valve/Shifioff Valve 5 . 0 x l 0 -« 

By-Pnea Valve 4 . 0 x lo-* 


Figure 4-7. Schematic and Failure Rates for Concept K (Cig>illary Device, 
Thermal Subcooling, Boost Pump, Uncooled Duct, Pumped 
Coolant) 


Concept K 


(-) 


Surge 


Vacuum 

Tank 


Pump 


Concept L 


H 


Boost 
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Figure 4-8. Schematic for Concept L 
(Capillary Device. Thermal Sub- 
cooling, Boost Pump, Uncooted 
Duct, Coolant Dumped) 


Figure 4-9. Schematic for Concept N 
(Capillary Device, Turbopump 
Subcooling, Uncooled Duct, 
Coolant Dumped) 



Failure Rates: ^ 

Throttle Valve 15. 0 x W® 


Shutott Valve 5. 0 x 10 “® 

Cooling Coil 0.5x10-6 


Boost 
Pumps ( 2 ) 


Figure 4-lOe Schematic and Failure Rates for Concept O (Capillary Device, 
Turbopump Subcooling, Cooled Duct, Coolant Pumped) 


Vacuum 
Pump 

Figure 4-11. SchematJc for Concept P (Capillary Device, Turbopump 
Subcooling, Cooled Duct, Coolant Dumped) 
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The reliability of any system is expressed by the equation 


R = e 


-XKt 


where 


X is the failure rate for a system comp<nient 
Kt is the time**environment factor product 

For the reliability of propellant feed systems the values of X for the components are 
shown in Figure 4 *2 to 4-11 below the concept schematics. The mission profiles 
considered were the one» two, and five bum missions for the Centaur D-IS vehicle. 

The Kt factors are shown in Table 4-41 where the environments of boost, main engine 
bum and coast are c<msidered. The environmental factor K, a measure of the severity 
of the environment, was determined in an earlier Titan/Cmtaur-Viking (Ref. 4-5). 

The analysis used Concept A as the baseline with other c<mcepts evaluated as modificatiois 
to this ccmcept. The results of the reliability analysis are summarized in Table 4-42 and 
indicate the reliability of the eleven concepts relative to each other. A second indicator is 
also shown, the mean number of missions between failures (MMBF) defined by iX/Sm. R). 
Concept H has the highest reliability rating; this is achieved by replacing both boost 
pumps with a pressurization system for turbopump NPSP. The pressurization system 
iiSS many components, as indicated in Figure 4-6, however reliability achieved by 

Table 4-42. Comparison of Relative Reliability for Concept IMder Study 


Concept 

One Bum 

Two Burn 

Five Bum 

R* 

MMBF** 

R* 

MMBF** 

R. 

MMBF- 

A 

0.999271 

1371 

0.998957 

958 

0. 998792 

827 

B 

0.999025 

1025 

0.998605 

716 

0. 998384 

618 

D 

0.999466 

1872 

0.999236 

1308 

0. 999114 

1128 

E 

0.999582 

2392 

0. 999403 

1675 

0.999307 

1443 

H (RL10A-3-3A) 

0.999670 

3030 

0.999503 

2012 

0.999415 

1709 

H (RL10*Cat I) 

0.999676 

3086 

0.999524 

2100 

0. 999425 

1739 

K 

0.998853 

871 

0.998360 

609 

0.998101 

526 

L 

0.998969 

969 

i 

0.998527 

678 

0.998294 

586 

N 

0.998528 

679 

0.997897 

475 

0.997564 

410 

0 

0.997910 

478 

0.997013 

334 

0.996542 

289 

P 

0.997927 

482 

0.997205 

357 

0. 996734 

306 


* ReUabtUty 

** MMBF 3 mean mlselODs between failure defined as 1/in R. 
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parallel paths exceeds that of two boost pumps. The reliability of all acquisition 
systems is lower than that of all settling systems by relative reliabilities of 0. 999712 
versus 0. 999540 for settling and acquisition components, respectively (single bum 
mission, Kt - 12. 26). This is further illustrated by the comparison of systems B and 
K where the difference in settling versus acquistion accounts for a delta of 154 MMBF 
for a one-bum mission. 

Other comparisons can be made on reliability of subsystems. Components in boost 
pump pressurization give it a reliability of 0. 999559 versus boost pump thermal 
subcooling 0. 999429 indloating the higher reliability without the throttle valve and 
thermal subcooler components. When considering turbopump NPSP, the boost pump 
alternative has a reliability of 0. 999428 versus 0. 999869 from the thermal subcooler 
without boost pump versus 0. 999958 for pressurization, reflecting the hig^ reliability 
of the latter Sj stem. Components added for a cooled propellant duct reduce the 
reliability by a factor of 0. 999497. For coolant handling, pumping the coolant back 
into the tank reduces the reliability by the Victor 1. . 999883 while dumping coolant 
overboard incurs no penalty. 

Referring again to Table 4-42, the analysis shows that c<mcept H using the Category 1 
engine gives the highest reliability and MMBF, while Concept O will give the lowest. 

Concept H is a significant improvement over the baseline Concept A. Failure rates 
used in the analysis are from RADC Notebook RADC-TR-75-22 (Reference 4-6), and 
from Reference 4-2. 

4. 5 ELECTRICAL POWER CONSUMPTION 

Electrical power consumption calculations were performed, and the results of this analysis 
are summarized in Tables 4-43, 4-44 and 4-45. Power consumption for the valves and 
sensors was neglected. The two main power requirements are for pumping cooling fluid 
back into the tank and for warming cryostored helium to usage temperature. Heat ex- 
changer power requirements are relatively high. Table 4-43, for the RLlOA-3-3 engine, 
shows that the highest power requirement is for the option that pumps turbopump thermal 
subcooling fluid and propellant duct cooling fluid back into the tanks (Concept O). This 
consumption is only slightly greater than Concepts D and M which only pump turbopump 
subcooler fluid back into the tank. Concepts B and K, pumping boost pump subcooler 
fluid back into the tank, have lower power consumption than the other concepts requiring 
electrical power. Concept H was not designed for the RLlOA-3-3 engine. Concepts A, 

E, L, O and P do not require electrical power. 

Table 4-44 for the RL10A-3-3A engine, shows that the highest power requirement is for 
the option that uses cryostored helitun Option Hx) and the associated heat exchanger to 
warm the pressurant. Option O, pumping turbopump thermal subcooling fluid and propellant 
duct cooling fluid back into the tank uses the next highest power requirement. Thi.' 
consumption is slightly greater than Concepts D and M which only pump turbopump sub- 
cooler fluid back into the tank. Concepts A, B, H2» K and L were not designed for the 
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Table 4-43, 


System Power Consumption With RLlOA-3-3 Engine 


I 

■| 

j 


(Power tn watt-hourt^ 


Power Conmimiitloa Elenoiit 


B 

D 

E 

H 

HX 

■a 

n 

n 

M 


0 

1 

1 ONE-BiniN MISSION | 

P'.unptng Power • CooUat itelurr 

Bl 

28.9a 

503.04 

- 

B 

B 

B 

28.98 

• 

503.04 

B 

607.80 

B 

Battery Power > Preemreiit 
Heat Exeluaiter 

■ 

- 

B 

- 

B 

Ki 

Ki 


B 

- 

- 


- 

Total Power Coesumptloe 

- 

28.98 

503.04 

B 

m 

m 

B 

m 

- 

503.04 

- 

mm 

- 

TWO-BUttN MISSION * | 

Pumptag Power * Coolant Return 

- 

30 « a 

649.82 

• 

B 

B 

w 

30. a 

B 

049.92 

B 

692.30 


Battery Power - Preaauraot 
Heat ExelMiicer 

m 

B 

- 

B 

WM 

Ki 

Ki 

- 

B 

- 

B 

- 

Bl 

Total Power Conaumptlon 


30.11 

649.92 

B 

m 

B 

B 

30.11 

B 

649.92 


692.30 

Bi 

1 FIVE-BimK MISSION | 

Pumplof Power * Coolant Retem 

- 

33.44 

782.88 

B 

m 

B 

B 

33.44 

- 

782.88 

B 

841.92 


Battery Power - Preaaurant 
Heat Exchanger 

■ 

- 

• 

B 

Ki 

WM 

WM 

- 

B 

- 

- 

- 

B 

Total Power Consumption 

- 

33.44 

782.88 

B 

B 

B 

B 

33.44 

- 

782. 88 

- 

841.92 



Table 4-44. System Power Consumption With RL10A-3-3A Engine 


(Power Id watt-houre) 


Power Consumption Element 

A 

mQiii 

D E 

H 

«»J 

■9 

Bl 

mm 

M 

N 

0 

■Bl 

ONE-BURN MISSION 

P imping Power - Coolant Return 

B 

B 

250. 08 

Bl 

B 

- 


B 

B 

250. «)« 

- 

254.38 

m 

Battery Power - Pressurant 
llrat CxehaiiKer 

X 

B 

- 

- 

B 

336 

Hi 

B 

B 

- 

B 

1 

B| 

Total Power Consumption 

/ \ 

B 

250.08 

- 

B 

336 

- 

B 

B 

250. 08 


254.38 

- 

1 TWO-BURN MISSION * | 

Pumping Power - Coolant Return 

w 

\/ 

311.04 

- 

- 

. 

NoT^ 

Analyz. 

\ / 

V 

311.04 


142.38 

- 

Battery Power - PressuranC 
Heat Exchanger 

mi 

B 

- 

B 

- 

864 

- 

B 

X' 

- 

“ 

- 

- 

Total Power Consumption 

III 

B 

311.04 

- 

- 

864 

- 

B 


311.04 

B 

342.38 

- 

Fn^E-BURN MISSION 

Pumping Power - Coolant Return 

B 

B 

406.06 

B 

B 

- 

rjQjjH 

B 

B 

400. 08 


449.76 

- 

Battery Power - Preasurant 
Heat Exchanger 

B 

WM 

B 

B 

B 

2016 

B 

B 

R| 

- 

B 

1 

- 

Total Power Consumption 

B 

T\ 

400.08 

B 

- 

2016 

B 

B 

B 

too. 04 

B 

449.76 

- 
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Table 4-45. System Power Consumption With RLIO-Category I Engine 


(Power In watt-hours) 


1 Powsr Coiuumptiun Element 

A 

D 

— 


■B 

Hi 

II 9 

IT 


l>!_ 

N 


P 

! ONE-bURN MISSION | 

Pumping Power - Coolant Return 


M 

161.76 

- 


- 

- 

m 

Bn 

181.78 

B 

166.08 

B 

Batti^iry Power - Presaurent 
Heat Eaehancer 


WM 

- 

■ 

B 

246 

2112 



- 

B 

B 

m 

Total Power Coneumptlon 

m 

m 

161.76 

B 

B 

248 

2112 

m 

m 

161.78 

- 

u 

165.06 

B 










Pumping Power - Coolant Return 

\/ 

\ / 

216.04 




_ 

\/ 

m 

215.04 


244.8 

- 

Battery Power - Preaaurant 
Heat Exehaneer 

X 


- 

- 


864 

2304 

X 

H 

- 

B 

B 

- 

Total Power Coneumptlon 

m 

ail 

216.04 

- 

- 

864 

2304 

/ ^ 


216.04 

- 

244.8 

- 










Pumping Pv/wer - Coolant Return 


BB 

273. CO 



B 

B 


m 

273.60 

- 

317.3 

B 

Battery Power - Pressurant 
Ik at Exi^banger 

gi 

Ki 

B 

B 

B 

2016 

2640 


ill 

- 

B 


** 

Total Power Consuntptlon 


A 

273.6 

- 

- 

2016 

2640 


/A 

273.6 

- 

317.3 

- 


RL10A-3-3A engine. Concepts E, H, N and P do not require electrical power. 

For the RLIO Category I engine, Table 4-45 shows that the highest power requirement 
is for the option that uses helium pressurization for the LH 2 tank and cryostored helium 
(Option H2). This is followed by the option using autogeneous LH 2 tank pressurization 
and cryostored helium. Option O, pumping turbopump thermal subcooling flow and 
propellant duct cooling flow back into the tanks is the next highest in electrical power 
constunption. This consumption is slightly greater than that required for Option D, 
pumping only turbopump subcooled flow back into the tank. Options A, B, K and L were 
not designed for the RLIO Category I engine. Options E, H, N and P do not use 
electrical power. 

4.6 MISSION PROFILE FLEXIBILITY 

Mission profile flexibility assessments were made. For systems using settling, added 
start sequence time will be required to accomplish settling. This will have an impact 
on the existing mission profiles for the Centaur D-IS. Main engine firing with capillary 
devices can be initiated more quickly than with settling thrusters. 

For other missions, capillary devices could have a limiting influence on minimum bum 
time and maximum time between bums. Minimum bum time, with a start basket 
capillary device Is the required time (in excess of the start sequence) to settle and refill 
the start baskets under main engine thrust. Capillary device thermal conditioning 
requirements are directly determined, from the maximum time between burns and are a 
factor in determining the capillary device volume. Large capillary device volumes are 
incompatible with short final bums because liquid level in the tank will be insufficient 
to allow the capillary device to refill. The considerations discussed in this paragraph 
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have no bearing on Centaur D-IS mission profile flexibility and are presented only for 
general information and extrapolation to other missions and vehicles. 

4. 7 C ONCLUSIOMS AND RECOMMENDATIONS 

For the Centaur D-lS. feed systems using capillary devices have the greatest mission 
profile flexibility. Feed systems having the lowest hardware weight are those using 
propellant settling and thermally suboooled boost pumps with coolant pumped back into 
the tank (Concept B). Feed systems havii^ the lowest pi^load weight penalty are 
Concepts A and B. Concept A uses propellant settling and pressure fed boost pumps. 

For missions with significantly greater burns than those required by the Centaur D>1S 
vehicle, capillary device payload weight penalty will be less than that of propellant 
settling. Concepts with no electrical power consumption are Concepts A,E,L,N and P. 
The highest reliability concept is Concept H which utilizes pressure fed turbopumps. 

The study shows that several areas are worthy of investigation depending on the 
direction taken by new vehicle design requirements. If high reliability is the major 
criteria, then pressure fed vehicles should be studied. If low payload penalty and low 
power consumption are most significant then the baseline Centaur D-IS using propellant 
settling and boost pumps (Concept Af is best. If low hardware weight is most important 
and missions of two burns or less are required payload penalty is then also lowest) 
then propellant settling with thermal subcooling and coolant returned to the tank should 
be selected. For missions greater than five bums (approximately 10 burns), capillary 
devices are attractive Rising thermally subcooled boost pumps with coolant returned 
to the tank, Concept K). 

Ihe study has shown that the baseline system is attractive and, under a specific set of 
assumptions, pressure fed turbopumps, thermally subcooled boost pumps with coolant 
pumped back into the tank and capillary devices are worthy of additional study. 
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APPENDIX A 


REFILLING PROGRAM DOCUMENTATION 


A.l PROGRAM CAPABILITY 

Program REFILL predicts start basket refilling conditions as a function of vehicle and 
mission conditionSt capillary device characteristics and engine requirements. A flow 
chart of the program is given in Figure 2-2. The program determines capillary device 
liquid level as a function of time. The analysis includes the effects of dynamic pressure » 
screen wicking, multiple screen barriers, window (standpipe) screens that can be of a 
different mesh than the main screens, time dependent liquid collection, and variable 
acceleration level due to changes in vehicle mass. Outflow from the tank is included 
either as an input or calculated based on feed system pressures. Vapor pullthrough 
height is accounted for as a function of tank outflow rate. Options are included for 
simulating spilling, wetted or dry standpipes and wicking between screen/plate barriers. 

The following pages present a listing of the program and a sample input and output. 

A. 2 PROGRAM LISTING 

The following pages give a brief description of the main program and each program 
subroutine followed by a listing of the corresponding program element. 

Listings are presented for the following program elements: 

DRIVER PROGRAM 

REFILL 

SUBROUTINES 

DVCOL 

FLOW 

IMPOUT 

INPUT 

OUTFLOW 

STAND 

SWET 

TIME 

WRIT 

DIAGNOS 
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PROGRAM REFILL 


llie main program serves as a driver for the flow subroutines after initializing 
boundary conditions for each time step. Flags for spilling, screen wetting, standpipe 
calculations, screen impingement, and outflow calculations are used to direct control 
to the appropriate subroutines. Vehicle acceleration level is calculated. Collected 
liquid height outside the basket is determined from tables. Impingement calculations 
for determining the dynamic pressure due to collected fluid velocity outside the start 
basket are performed. Pullthroug^ tables are used to determine if the liquid level in 
the basket has dropped below the point at which vapor will enter the channels. 

Liquid levels inside and outside the basket are calculated. Initial conditims for the 
first time step are determined and printed. 


PROGRAM REFILL ( INPUT . TAPESsiNPUT , OUTPUT . TAPE6»0UTPUT ) 

PROGRAM COMPUTES START BASKET REFILLING UNDER VARIABLE DRIVING HEAD. 
VARIABLE WETTING.VARIABLE DYNAMIC PRESSURE. MULTIPLE COMPARTMENTS. AND 
MULTIPLE SCREEN BARRIERS 

CALCULATED 4UANTITIES-INTERMEDIATE 

AL : AREA FOR LIQUID INFLOW. FT2 
ALI : AREA ADJACENT TO LIQUID INSIDE THE BASKET. FT2 
ALO = AREA ADJACENT TO LIQUID OUTSIDE THE BASKET, FT2 
AV = AREA FOR VAPOR OUTFLOW ACROSS UNWETTED SCREEN. FT2 
AVW= AREA FOR VAPOR OUTFLOW THROUGH WETTED SCREEN ADJACENT TO VAPOR, FT2 
B= SPACING BETWEEN SCREEN BARRIERS, FT. IN SUBROUTINE SWET 
CV0L(N)= COLLECTED VOLUME AT TC(N), FT3 
DPIN= BASKET INTERNAL PRESSURE MINUS ULLAGE PRESSURE, PSF 
DPS = MAXIMUM SURFACE TENSION DELTA P.PSF 

DPS2=SURFACE TENSION PRESSURE DIFFERENTIAL FOR TOP SCREEN, PSF 
G=AMBIENT ACCELERATION, FT/SEC2 
HI=LIQUID LEVEL INSIDE BASKET, FT 

HLEV^COMPUTED LIQUID LEVEL WHEN START BASKET AND TANK LEVEL ARE EQUAL, FT 
HO=LIQUID LEVEL OUTSIDE BASKET, FT 

HS - HEAD THAT CAN BE SUPPORTED BY SURFACE TENSION, FT 
IFBAD =FLA6 FOR INTERPOLATION CALCULATIONS. IFBAD*! FOR GOOD EXIT, =2 BAD 
QL = LIQUID VOLUME INFLOW RATE, FT3/SEC 
T = TIME, SEC 

TVOL =T0TAL FLIUD VOLUME IN THE TANK DURING ONE TIME STEP,FT3, 

VCL = COLLECTED VOLUME AT A GIVEN TIME STEP.FT3 
VCL1= COLLECTED VOLUME FROM CURRENT TIME ST0P.FT3 
VCL2= COLLECTED VOLUME FROM LAST TIME STEP.FT3 

VIL=L1QU1D VOLUME INSIDE THE START BASKET AT START OF TIME STEP.FT3 
VOL=LIQUID VOLUME OUTSIDE START BASKET AT START OF TIME STEP.FT3 
VOREF = VOLUME FLOW INTO THE BASKET DURING THE LAST ITERATION, FT3 
VV = VAPOR INFLOW VELOCITY ACROSS UNWETTED SCREEN, FT/SEC 
WT= WETTED HEIGHT OF SCREEN 

C0MM0N/FL0W/AL1,AL2,AV1,AV2,C0NVERG,C0NVI.C0NV2,DPIN1,DPIN2,H0LI, 

1H0L2, NITER 

:OMMON/IMPINGE/A,AREA(20),ATOT,B,NDR,NIM,NW,VOLO(20),Y,Z 
COMMON/IMPOUT/CVOL,DPS,HDPIN,HPULL,VOL 

COMMON/INFILL/CVI(20),DBP,DBP2,HCOL(20),HPT(lO),J,K,MPROP,MR,MVEH, 
lN,NG,NS,NSP,QOUTa0),TC(20),THRST,TIM(20),TITLE(7),VILI,VTOT(20) 
COMMON/INPUT/DT,G,HT(20),HTOT,I,KC,NOUT,NST,OA,QP,RHOL,SIGMA, 

1VOLK20) 

C0riM0N/REFILL/AL,AV,AVW,DPIN,DPS2,HI,HILl,HIL2,H0,HS,HX,JUMP,QL. 
1Q0,T,VCL,VIL,VILMAX,V0,V0L0UT,V0REF,VV,WT 
COMMON/TIME/HTOL.TMAX 

COMMON/WRIT/ADL,AVL,AVLS,AVS,AVT,AVWS,HDYN,VL,VT,VVL,VVLS,VVS,VVW, 

IVVWS 

REAL KC.MPROP.MR.MVEH 
C JUMP IS A VARIABLE USED TO RESTART NEW CASES MHBlO/77 

i/ ASSIGN 10 TO JUMP 
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10 CALI INPUT 

ADLsAL3AV:AVL=AVT3AVUsC0NV1sC0NV2sDPZNsDPXN1sDPIN2sHDYNsHLEVs 
IQL sQOsTsVCLsVlLnAXsVlsVOsVOLOUTsVOReP^VTBVVsVVLsVVUsMTsO 
URITE(6.20) (TITLE(L).Lsl,7) 

20 FORMAT(1H1»7A10) 

IP(NS.LE.l) GO TO SO 

AsNSMANO.05 

B3NSKBM0.9S 

YSNSXYW0.95 

Z:NSKZK0.9B 

C INITIALIZED BASKET NEIONT NMMMMMMMMMNMMMMMMMNKMHNMMiiMKMNMMMMNMilMMMnHB10/77 
SO CALL TABL(VILI.HT>VOLia).HT(l).ia>l.I»IFBAD) 

VILsVUI 

CALL TABL(T.HO.TC(l),HCOL(l).l,l.l.N.IPBAD) 

CALL TABL (HO . VCLl . HT( 1 ) . VOLO( I ) ,1 . 1 »1 . 1 > ZFBAD) 

VOLsVCLl 

ir(NSP.EQ.O. 0 R.HI.lt. HO) 00 TO AO 
VILI=V1L1*VCL1 

CALL TABL(VILI>HLEV.VTOT(l)»HT(l)aa.iaaFBAD) 

CALL TABL(HLEV,VOL,HT(l).VOLO(l).l.ia>iaFBAD) 

CALL TABL(HLEV,VIL.HTa).VOLia)»iaaa»IFBAD) 

HO^HLEV 

HI=HLEV 

AO URITE(A.50) T.HO.HI 

50 FORMAT(>« INITIAL CONDITIONS. TIME sxpB.AnSEC, TANK LIQUID LEVEL 
I =«FS.AKFT, BASKET LIQUID LEVEL s»F8.A>«FT«) 

CALL TABL(HC0L(K).CV0L.HT(1),V0L0(1).1.1.1»I.IFBAD) 

DPS:8IGMAM2 . 72NS . 0 ABE05/DBP 
1F(DBP2.EQ.0.)DBP2:DBP 

DPS2=S10nAX2 . 72*3 . 0ABE05/DBP2 MHBlO/77 

IF(NDR.EQ.l) DPS2*0 
00 TO 120 
100 T=T*DT 
VCL2SVCL1 

CALL TABL(T.H0.TC(1).HC0L(1)»1»1.1.N»IFBAD) 

CALL TABLCHO. VCLl. HT(1).V0L0(1) .1.1.1. I. IFBAD) 

VCL=VCLl-VCL2 
TVOL=VOL*VIL 

CALCULATION TO DETERMINE LIQUID LEVEL OUTSIDE AND INSIDE BASKET IFMHBlO/77 
LEVELS ARE EQUAL MKK«()(KKM)on(K)(iHO(«(»K)(i(MXinn(KXKioo(MKKtooo.ioo(t(K«O(K)O(MHB10/77 
CALL TABL(TV0L.HLEV.VT0T(1).HT(1).1.1.1.1.1FBAD) 

CALL TABL(HLEV.VIL1.HT(1).V0LI(1).1.1.1.I.IFBAD) 

CALL TABL(HLEV.V0L1.HT(1).V0L0(1).1.1.1.I.1FBAD) 

CALL T ABL ( HTOT . VTOP . NT( I ) . VTOT (1 ) . 1 . 1 . 1 . 1 . IFBAD) 

VOLNAX=VOL‘VOLl 

IF(VOLMAX.LT.VOL-VTOP) VOLMAX=VOL-VTOP 
IF(VOREF.GT.VOLMAX) VOREF=VOLMAX 
VILMAX=VIL1-VIL 

C CALCULATIONS IF LIQUID HAS BEEN COLLECTED OUTSIDE BASKET xxxxx*xx«MHB10/77 
IF(VOREF.GT,VILMAX+VOLOUT) VOREF=VILMAX+VOLOUT 
VILsVIL+VOREF-VOLOUT 
VOLsVOL+VCL. 

IF(VOL.LT.O.)VOL=0. 

IF(VIL.LT.0.)VIL=0. 

CALL TABL(V0L.H0.V0L0(1).NT(1). 1.1. 1.1. IFBAD) 

CALL TABL ( VIL.HI. VOLK l).HT(l). 1.1.1. I. IFBAD) 

CALL TABL (QP.HPULL.Q0UT(1).HPT(1).1.1.1.J. IFBAD) 
IF((NOUT.NE.O).A.(HI.LT.HPULL))GO TO 3007 
IF(N1TER.GT.10)CALL DIAGNOS 
CALL TIME 

IF(T.OE.TMAX) 60 TO 10 
120 IF(NIM.EQ.O)GO TO 103 

CALL TABL(T.VJ.TIM(l).CVl(l).l.l.l.NlM.irAAD) 

F*2. 

BK 

C VO IS THE IMPINGEMENT VELOCITY CALCULATED FROM SYMONsS CORRELATION******** 
V0 :(-FkBK*A 4SQRT((FKBK*A)*M24‘A.*(RH0L*VJ**2)*(F*BK*B«RH0L*BK)))/(2 
l.*(F*BK*B4RH0L*BK)) 

103 IF(NG.EQ.l) GO TO 150 
MPROP=MPROP-QP*RHOL*DT 

G = 32.17*TI,RST/(MVEH4MPR0P*(14MR)) RlfRODUCIBILITY OF 'iu- 

A -3 OiaClNAL JAGB IS POOR 



ISO HS:32.17KDPS/(Gi(RHOL) 

HX=HS 

! F( NST . NE . 0 ) HX:32 . 1 7NDPS2/ (GNRHOL ) 

IFdlO.OT.HTODHOsHTOT 

C SELECTION OF UETTING CALCULATIONS N)IMMNNNNNKXNNMMMKi(KXKKMMNilNKMRMKnHB10/77 
TF(NVI.EG.ni<IT=HO 
IF(NU.EQ.2)WT=HT0T 
IF(NU.EQ.3)CALL SUET 

^ OUTER LIQUID BELOU TOP OF BASKET AND EQUAL TO WETTED HEIGHT XMXMKxnHBlO/77 
IF<HO,EQ.MT.A.HO.LT.HTOT)GO TO 300 

C OUTER LIQUID AT TOP OF BASKET AND EQUAL TO WETTED HEIGHT MMXkKMMHBlO/77 
IF(HO.EQ.UT.A.HO.EQ.HTOT)GO TO 3A0 

C WETTED HEIGHT GREATER THAN LIQUID LEVEL AND BELOW TOP OF BASKETxnmMHBIO/77 
IF(UT.6T.H0.A.H0.LT.HT0T)G0 TO 350 

C WETTED HEIGHT GREATER THAN LIQUID LEVEL AND EQUAL TO BASKET HEIGHTnHBIO/77 
IF(UT.GT.HO.A.WT.GE.HTOT)GO TO 3A5 
WRITE(S.200)HO,WT.HTOT 

200 rORHATCM FRROR IN COMPUTING HEIGHTS. HO^MFIO. 7, N WT=XF10.7.X HTOT^M 

iriO.7) 

GO TO 10 

C CASEl xxxx HO.LT.HTOT HO.EQ.WT XRXR XRXX 

300 CALL TABL(HO.ALO.HTa).AREA(l)»l.l.l.I.IFBAD) 

CALL TABL(HI.ALI.HT(l).AREA(l).l.I.l.I.IF6AO) 

IF(NOUT.NE.O)GO TO 350 

AL=ALO-ALI 

AV=AT0T-AL0 

IF(NIM.NE.O) CALL DVCOL 
IF(NIM.NE.O)GO TO 100 
!F(HO.GT.HI) call FLOW 
IF(HO.LE.HI) VV=0.0 
QL=AV«VV«OA 
0riN=AKVV+BxVVXX2 
HDPIN=DP1NX32.2/(RH0LNG) 

IF(HOPlN.GT.(HO-HI)/lO.)GO TO 350 
V0REF=QLXDT 

1F(V0REF.GE.VILMAX)G0 TO 100 
CALL OVCOL 
GO TO 100 
3A0 AVW=0.0 
3A5 AV=0.0 

350 IF(NIM.NE.O) CALL DVCOL 
IF(NIM.NE.O)GO TO IOC 
CALL FLOW 
GO TO 100 

3007 WR1TE(6,3017)HPULL 

3017 FORMAT(x PULLTHROUGH OCCURRED, HPULL=XFO. 5) 

GO TO 10 
END 


■^f ... 
T" ■■ • 
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SUBROUTINE DVCOL 


Subroutine DVCOL calculates flow when liquid level inside and outside the start basket 
are unequal at the end of the time step and net refilling occurs. Screen impingement 
calculations are included in the pressure balances. Standpipe areas are incliuled in the 
calculation. Convergence calculations are performed. 


SUBROUTINE DVCOL 

CALCULATIONS UHEN LIQUID LEVEL INSIDE AND OUTSIDE THE START BASKET ARE 
UNEQUAL AT THE END OF THE TIME STEP AND REFILLINO OCCURS MKMMMMNNKnHBlO/T? 
DVCOL ^SUBROUTINE TO COMPUTE FINAL HEAD UHEN VOREF. LT. VILNAX 
AVI : INITIAL UNUETTED VAPOR FLOU AREA.FT2 
AV2 = FIMAl UNUETTED VAPOR FLOU AREA.FT2 

CONVI: CONVERGENCE VARIABLE FOR CONTINUITY EQUATION-CURRENT ITERATION 
C0NV2- CONVERGENCE VARIABLE FOR CONTINUITY EQUATIOH-LAST ITERATION 
CONVERGE CONVERGENCE VARIABLE THAT MUST BE UITHIN ZERO 40R-RANGE 
DPlNls STORAGE FOR CURRENT DELTA P ITERATION 
DPIN2S STORAGE FOR CURRENT DELTA P ITERATION 

DPIN : PRESSURE DIFFERENCE OF GAS INSIDE AND OUTSIDE BASKET KNMKMK 
FOR FLOU AND DVCOL »DPIN MUST BE POSITIVE AND THE GAS PRESSURE INSIDE 
THE BASKET EXCEEDS THE PRESSURE OUTSIDE THE BASKET xKXKXNNKNimHKNH 
HALl s INITIAL DISTANCE FOR LIQUID INFLOU.FT 
HALE 2 FINAL DISTANCE FOR LIQUID INFLOU.FT 
HAVL12 INITIAL DISTANCE FOR VAPOR FLOU INTO LIQUID. FT 
HAVL2S FINAL DISTANCE FOR VAPOR FLOU INTO LIQUIO.FT 
HAVUIs INITIAL DISTANCE FOR VAPOR FLOW ACROSS UETTED SCREEN, FT 
HAVU22 FINAL DISTANCE FOR VAPOR FLOU ACROSS UETTED SCREEN. FT 
NITER 2ITERATION COUNTER 
VORl= STORAGE FOR VOREF 

UT1.UT2:UETTIHS PARAMETERS AT START AND END OF TIME STEP 
COMMON/FLOU/ALl. ALE. AVI. AV2.C0NVERG. CONVI. C0NV2.DPIN1.DPIN2.H0L1, 

1H0L2. NITER 

C0Mt10N/DVC0L/AST.A2.B2. RANGE 

COKMON/IMPINGE/A. AREA(20).ATOT.B.NDR.NIM.NW.VOLO(20).Y,2 
CONNON/IMPOUT/CVOL.DPS.HDPIN.HPULL.VOL 

COMMON/INPUT/OT.G.HT(20).HTOT.I.KC.NOUT.NST.OA.QP.RHOL.SI6MA, 

1V0LK20) 

C0MM0N/REFILL/AL.AV.AVU.DPIN.DPS2.HI.HILI.HIL2,H0.HS.HX.JUMP,QL, 

IQO.T.VCL. VIL.VILMAX. VO. VOt OUT. VOREF. VV.WT 
COKMON/STAND/AVLR . AVR . AVUR 

CONMON/URIT/ADL.AVL.AVLS.AVS.AVT.AVUS.HDYN.VL.VT.VVL.VVLS.VVS.VVU, 

IVVUS 
REAL KC 
URITE(6.5) 

5 FORMAT(N SUBROUTINE DVCOLN) 

NITER=0 

IF(MOUT.NE.O) DPINsO.OOl 
10 VORlsVOREF 

IF(NOUT.EQ.O)VOLOUT=0. 

VILX2=VIL+V0REF-V0L0UT 
V0LX2=V0L-V0REF 

200 CALL TABLCVIL.HILl.VOLKD.HUn.l.l.l.I.IFBAD) 

CALL TABL(VOL.HOL1.VOLO(1).HT(U.1.1.1,I.IFBAD) 

CALL TABL(VILX2.HIL2.V0LI(1).HT(1).1.1.1.I,IFBAD) 

CALL TABL(V0LX2.H0L2.V0L0(l).HT(l).l.l.l.i,IFBA0) 

IF(HIL2.GE.HT0T)HI12=HT0T 
IF(H0L2.GE.HT0T)H0L2=HT0T 
IF(H0L1.GE.HT0T)H0L1=HT0T 

CALL TABL(H0L1.AVLT1.HT(1).AREA(1).I.1.1.I.IFBAD) 

CALL TABL(H0L2.AVLT2.HT(1),AREA<1).1,1.1,I,1FBAD> 

CALL TABL (HILl. AVTl.HT(l). AREAL 1), 1.1,1. 1. IFBAD) 

CALL TABL(H1L2.AVT2.HT(1),AREA(1).1.1,1.I.IFBA0) 

C AREA AND LIQUID LEVEL CALCULATIONS NMKMMKIH(NMNKMNMNIiNIUINMNMMMNMKMMMHB10/77 
AVTs(AVTl*AVT2>/2. 
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H0SH0L2 

C SELECTION OF UETTINO CALCULATIONS MMKMMNMMNMMNKMMMMMKKMKMkMMMKMMMMnHSlO/77 
1F(NU.EQ.1)UT:H0 
IF(NU.EQ.2)UT:HT0T 
IF(NU.EQ.3)CALL SUET 
UT2=WT 
H0:H0L1 

IF(NU.EQ.1)UT:H0 
IF(NU.E0.2)UT:HTOT 
IF(NU.EQ.3)CALL SUET 
UT1:UT 

CALL TABL(UTl»AVUTl.HT(l),AREA(i)a.l.l.I»IFBAO) 

CALL TABL(UT2,AVUT2,HT(1).AREA(1).1.1»1.I»1FBAD) 

AVlsATOT-AVUTl 
AV2=AT0T-AVUT2 
20 CONTINUE 

HDPIN:DPINm32.2/(RH0Lm6) 

IF(NIII.NE.O) GO TO 1 
HAVLlsHILl+HS 
HAVL2SHIL2+HS 
GO TO 11 

1 CONTINUE 

ROUTINE CALCULATES BASKET REFILLING. SIHILAR TO DVCOL* LIQUID LEVELnHBlO/77 
INSIDE AND OUTSIDE THE BASK"T UNEQUAL AFTER A TIKE STEP - NET MHBlO/77 
REFILLING OCCURS. VAPOR OUTFLOU TO ALLOU LIQUID TO ENTER>»««««»»««««MHB10/77 
DYNAMIC PRESSURE CALCULATED FROM IMPINGEMENT VELOCITY k*(xi«»(m»o«»*««»«umHB10/77 
URITE(6.6) 

6 FORMAT(» IMPINGE*) 

DYNAMIC PRESSURE CALCULATED FROM IMPINGEMENT VELOCITY *>(**>Oi******MHB10/77 
DPDYN=A«V0+B‘<V0*«2 
HDYN=0P0YN>«32.2/(RH0L*G> 

HAVLlsHILHHS + HDYN 
HAVL2SHIL2+HS+HDYN 
11 IFCHAVLl.CT.HOLDHAVLlsHOLl 
IF(HAVL2.GT.H0L2)HAVL2=H0L2 
HAL1=H0LI-HDP1N 
HAL2SH0L2-HDPIN 

CALL TABL(HAVL1.AVL1.HT(1).AREA(1).1.1.1.I.IFBAD) 

CALL TABL(HAVL2.AVL2.HT(1).AREA(1).1.1.1.I.IFBAD) 

CALL TABL(HAL1.AL1.HT(1).AREA(1).1.1.1.I.1FBA0) 

CALL TABL(HAL2. AL2.HT<1).AREA(1).1.1.1.1.IFBAD) 

AVW1=AVUT1-AVLTI 
AVW2=AVWT2-AVLT2 
AL1=AL1-AVT1 
IF(ALl.LT.O) AL1=0 
AL2SAL2-AVT2 
IF(AL2.LT.0) AL2=0 
AVLI=AVLT1-AVL1 
AVL2=AVLT2-AVL2 
IF(AVL1.LT.0. )AVL1=0. 

IF(AVL2.LT.0. )AVL2=0. 

IF(NIM.EQ.O) GO TO 2 
HADL1=H1L1*HDPIN-HDYN 
HADL2=HIL2+HDPIN-N0YN 
IFCHAOLl.GT.HILDHADLlsHILl 
IF(HADL2.GT.HIL2)HADL2:HIL1 
IF(HADL1.LT.0)HADL1=0. 

IF(HADL2.LT.0)HADL2:0. 

CALk IABL(HADL1.ADL1.HT(1).AREA(1).1.1.1.I.1FBAD) 

CALL TABL(HA0L2.ADL2.HT(1).AREA(1).1.1.1.I»IFBAD) 

ADL3=AVTI-ADL1 
ADIA=AVT2-ADL2 
ADL1:ADL3/2.*AVLT1-AVT1 
ADL2=ADLA/2.*AVLT2-AVT2 
IF(ADL1.LT.0.)ADL1:0. 

IF(ADL2.LT.0. )ADL2s0. 

ADL=(ADL1*ADL2)/2. 

2 CONTINUE 
AV:(AVl*AV2)/2. 

AVU:(AVUUAVU2)/2. 

AVL=(AVLl*AVL2)/2. 

AL =(ALl*AL2)/2. 

IF(OPIN.LE.OO)VVsVVSsO.O . „ 
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VAPOR FLOW ACROSS DRY SCREEN.OUT OF BASKET MKKXXKMXKXXKKXKMMMMKMMKMHBIO/?? 
VVS =VAPOR VELOCITY FOR UNWETTED STANDPIPE xxxxxxxxxxxxMHBlO/77 

IF(OPIN.GT.O.O.AND.NST.NE.O)VVS=(-A2+SQRT(A2XX2+4XB2XDPIN))/(2XB2) 
IF(DPIN.6T.O.O)VV=(-A+SQRT(AJ««+A>«BxDPIH))/(2XB) 

IF(NIM.NE.O) GO TO 3 

IF(NST.EQ.O)GO TO 83 NHBlO/77 

IF(DPIN-DPS2.LE.O.)GO TO 82 

VVUS=VAPOR VELOCITY FOR WETTED STANDPIPE XXXXXXXXXXXXMHBlO/77 

VVWS:(-A2+SQRT(A2XX2+A.XB2X(DPIN-DPS2)))/(2.«B2) MHBlO/77 

VVLS=VAPOR VELOCITY FOR STANDPIPE SURROUNDED BY LIQUID xxxxkxxkmHBI 0/77 
VVLS=(-A2+SQRT(A2«»2+2.J<B2X(DPIN-DPS2)))/(2.XB2) MHBlO/77 

GO TO 83 MHBlO/77 

82 VVWS=0. 

84 VVLS=0. 

S3 CONTINUE MHBlO/77 

IF(DPIN-DPS.LE.O.)GO TO 25 

VAPOR FLOW ACROSS WETTED SCREEN INTO VAPOR, OUT OF BASKETXxxxxxxxxxMHBlO/77 

VVW=(-A*SQRT(AXX2+4.»BX(DPIN-DPS)))/C2.XB) 

VVL=(-A*SQRT(Axx2+2.xBx(DPIN-DPS)))/(2.xB) 

GO TO 26 

25 VVW=0. 

VVL=0. 

26 CONTINUE 
GO TO 12 

3 CONTINUE 

1F(NST.EQ.O)GO TO 87 
IF(DPIN-DPS2.LE.O.)GO TO 86 

VVWS=VAPOR VELOCITY FOR WETTED STANDPIPE xxxxxxxxxx)OfMHB10/77 

VVWS=(-A2+SQRT(A2xx2+4.xB2x(DP1N-DPS2)))/(2.xB 2) MHBlO/77 

GO TO 87 
86 VVWS=0.0 


87 CONTINUE 

IF(DPIN-DPS.LE.O.)GO TO 23 

VAPOR FLOW ACROSS WETTED SCREEN INTO VAPOR, OUT OF BASKETXXXXXXXXXXMHBl 0/77 
VVW=(-AtSQRT(AXX2+4.J<BX(DPIN-DPS)))/(2.xB) 

GO TO 24 

23 VVW=0.0 

24 CONTINUE 

1F(DPIN-DPS-DPDYN.LE.0.)G0 TO 27 

VAPOR FLOW INTO LIQUID, OUT OF BASKET xxkxxxxxxxxxxxkxxxxxxxxxxxxxxMHBIO/ 77 
DYNAMIC PRESSURE USED TO CALCULATE LIQUir> AND VAPOR FLOW xxxxxxxxxMHBlO/77 


IF(NST.EQ.0)6O TO 85 MHBlO/77 

IF(DPIN-DPS2-DPDYN.LE.O.)GO TO 88 

VVLS=VAPOR VELOCITY FOR STANDPIPE SURROUNDED BY LIQUID xxxxxxxxMHBl 0/77 

VVLS=(-A2*SQRT(A2XX2+2.xB2x(DP1N-DPS2-DPDYN)))/(2.»B2) MHBIO/77 

GO TO 85 
88 VVLS=0.0 

85 CONTINUE MHBlO/77 

VVL=(-A+SQRT(AXX2+2.«BX(DPIN-DPS-DPDYN)))/(2.>iB) 

GO TO 31 

27 VVL=0. 

31 CONTINUE 
12 CONTINUE 

H0PIN=2.XHDPIN 

IF(H0L1-HIL1+H0L2-HIL2-H0PIN.LE.0. )GO TO 28 

LIQUID FLOW INTO VAPOR, INTO BASKET xxxkxxxxxxxxxxxxxxxxxxxxxxxxxxxmHBIO/ 77 
VL=(-Y+SQRT(YXX2+ 2x( H0L1-HIL1+H0L2-HIL2-2. KHDPIN)XRH0U<G/32 . 2 ) ) 

1/(2. XZ) 

LIQUID FLOW INTO LIQUID. INTO BASKET xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxMHBlO/77 
VT=(-Y+SQRT(Yxx2+2.xZX(H0L1-HIL1+H0L2-HIL2-2.xHDPIN)XRH0LxG/32.2)) 

1/(2. XZ) 

GO TO 29 

28 VL=0. 

VT = 0. 

29 CONTINUE 

xxxx CONVERGENCE CALCULATIONS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxKXxxxx 
DPIN1=DPIN 

if(nst.eq.o.and.n:m.eq,o) go to 93 

IF(NST.EQ.O.AND.NIM.NE.O) GO TO 103 
HI=(HILl+HIL2)/2. 

CALL STAND 
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IF(NIM.NE.O) GO TO ^ 

CONVER6-AVRNVV'^AVS)(VVS^AVUR)(VVUUVWS)<VVUS'^AVLR)<VVl.^AVLSXVVLS-ALKVLi1HB10/77 
1-AVT?(VT HHBlO/77 

C0SVERG::C0NVERG)<2. 

QVO=AVR«VV+AVS«VVS+AVWRXVVW+AVWS«VVWS+AVLSi<VVLS+AVLRJCVVL 
GO TO 99 MHBlO/77 

93 CONTINUE MHBlO/77 

C0NVERG=(AVl+AV2)i<VV+(AVUl+AVW2)*<VVW+(AVLl+AVL2)J<VVL-(ALl+AL2)MVL- 
1AVT><VT 

OVO=AVxVVrAVW}<VVM>AVLXVVL 

99 CONTINUE MHBlO/77 

Q0 = 0. 

IF(NOUT.NE.O)CAU OUTFLOW 
QIH=VLX(ALl+AL2)><0A/2.+VT>(AVT>«0A 
QL=QIN 
GO TO 13 
A CONTINUE 

C0NVER6= AVRxVV’^AVSxVVS+AVWRXVVW+AVWS)<VVWS+AVLSmLS+AVLRXVVL-ALMVMHBU/77 
1L-ADL?<V0-AVT^VT MHBlO/77 

C0NVfc'RG=C0KVERGJ^2. 

QV0=AVRKVV+AVS>tVVS+AVWR«VVW+AVW5»VVWS+AVLSXVVLS+AVLR«VVL 
GO TO 9S 
103 CONTINUE 

CCNVERG=(AV1+AV2)XVV+(AVW1+AVW2)>CVVW+(AVL1+AVL2)XVVL-(AL1+AL2)MVL 
1-(ADL1+ADL2))<V0-2.5<AVT5<VT 
QV0 = AVXVV + AVW*^VVW+AVL«VVI 
9« CONTINUE 

IF(NOUT.EQ.O)QO=0. 

IF(NOUT.NE.O)CALL OUTFLOW 

QIN = (ADLl+ADL2)><OA5<VO+(ALl+AL2)J<OAXVO+AVTJ<VT><OA 
QIN=QIN/2 
13 CONTINUE 

IFCNITER.GT.DGO TO 73 
IF(QO.EQ.O*)GO TO 73 

C DETERMINES IF NET BASKET OUTFLOW EXISTS 
IF(QO.GT.QIN)GO TO 71 
73 CONTINUE 

C0NVERG=C0NVERG^0A/2+Q0 

CONV1=CONVERG 

IF(CONVERG.EQ.O..A.VL.EQ.O.)GO TO 35 
ir(ABS(CONVERG),LT.RANGE)GO TO AO 

35 CONTINUE 

IF. NITER. EQ.O)GO TO 30 9 

IF(WT.LT.HTOT.OR.NDR.EQ.l) GO TO 36 

HCH=(H0Ll+H0L2-HILl-HIL2)/2. 

IFCHX.GE.HCIDGO TO 37 

36 CONTINUE 

IF(NITER.GE.75)GO TO 300 
IF(C0NV1.EQ.C0NV2)G0 TO 30 
NITER=NITER+1 

IF(NITER.GT.10)CALL DIAGNOS 

i)PIN = (DPIN2^CONVl-DPINlxCONV2)/(CONVl-CONV2) 

HDPIN = D?IN>t'52.2/RH0L/G 
IF(NW.EQ.2,AND.NDR.NE.l) GO TO A3 
IF(WT.GE.HTOT.AND.NDR.NE.l) GO TO A3 
IF(AVS.NE.O) GO TO A1 
A3 CONTINUE 

H0PIN=HDPIN><2. 

IF(HOPIN*GT.(HOL1-HILHHOL2-HIL2))HOPIN=HOL1-HIL1+HOL2-HIL2-(HOL1- 

1HIL1+HOL2-HIL2-HX)/(NITER»10000) 

HDPIN=H0PIN/2. 

IF(HDPIH.LE.HX)HDPIN=HX+(HOLI-HIL1+HOL2*- HIL2-HX)/(NITER«10000 ) 

GO TO A2 

A1 IFCDPIN.LE.O. )HDPIN = 1./(NITER>^10) 

A2 CONTINUE 

DPIN = HDPINXRH0LJ<G/32.2 
C0NV2=C0NV1 
DPIN2=DPIN1 
GO TO 20 
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30 DPIN2=DPIN1 
C0NV2=C0HV1 
NITER=NITER*1 

lF(NirER.GT.10)CALl DIAGNOS 
DPIM=OPIM/2. 

HDPIN=0PIN»32.2/(RH0L»<6) 

IF(NITER.EQ.1.A.WT.GE.HTOT.A.HX.NE.O.)H?PINSHX 
DPIN=HDPlHKRH0LxG/32.2 
60 TO 20 
AO CONTINUE 

C KKM OUTFLOU RATE CALCULATION KliXMKiniKKMMXXKXMXKNKKKXitKXXKXKKMMKK 
IF(NIN.NE.O) GO TO 7 
QL = (ALl*AL2)/2.XVLX0AmxAVT«0A 
VOREF=OLXDT 
GO TO lA 
7 CONTINUE 

QL=(AL1+AL2)/2.»VLX0A+(ADL1*ADL2)/2.XW0X0A+VTXAVTX0A 
VOREF=QLXDT 

IF(NOUT.NE.O)CALL OUTFLOU 
lA CONTINUE 

IF(ABS(VOREF-VORl).LE.RANGE)GO TO 72 
V0REF=(W0REF+V0Rl)/2, 

IF(NITER.GT.IO)CALL DIAGNOS 
CO TO 10 

37 IF(NOUT.NE.O)GO TO 71 
URITE(6.38) 

38 FORMAKX HS IS GREATER THAN H0L1-HIL1*H0L2-HIL2. NO SOLUTION FO 
IR HDPINX) 

VOREF=0. 

IF(VOLOUT.HE.O.)GO TO 72 

IF((HOL1.LT.HTOT).A.(VOL+VIL.LT.CVOL))GO TO 72 

GO TO JUMP, TRANSFER CONTROL TO MAIN PROGRAM-STATEMENT 10, IN ORDERMHBlO/77 
TO RESTART A NEU CASE 
60 TO JUMP, (20) 

301 FORMAKX TOO MANY ITERATIONS IN SUBROUTINE DVCOL x) 

300 URITE(6,301) 

GO TO JUMP, (20) 

71 CALL IMPOUT 

72 RETURN 
END 
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SUBROUTINE FLOW 


Subroutine FLOW calculates flow when the liquid level inside and outside the basket are 
equal at the end of the time step. (Generally* subroutine FLOW would only be used if 
collection was complete, the screens surrounded by v^x>r were dry* or the flow rate 
across screens was relatively hig^ considering the time st^. ) Net refilling occurs. 
Djnnamic pressure (screen impingement) effects are not included. Standpipe areas are 
included in the calculation. Convergence calculations are performed. 

SUBROUTINE FLOM 

CALCULATIONS MHEN LIQUID LEVEL INSIDE AND OUTSIDE THE START BASKET ARE 
EQUAL AT THE END OF THE TIME STEP AND REFILLING OCCURS 
C0KV1= CONVERGENCE VARIABLE FOR CONTINUITY EQUATION- CURRENT ITERATION 
C0NV2= CONVERGENCE VARIABLE FOR CONTINUITY EQUATION- LAST ITERATION 
CONVERG=CONVER6ENCE VARIABLE FOR CONTINUITY EQUATION. ZERO REQUIRED 
DPINl= STORAGE FOR CURRENT DELTA P ITERATION 
DPIN2= STORAGE FOR LAST DELTA P ITERATION 

FOR FLOW AND DVCOL.DPIN MUST BE POSITIVE AND THE GAS PRESSURE INSIDE 
THE BASKET EXCEEDS THE PRESSURE OUTSIDE THE BASKET 
HAL = HEIGHT U'HERE LIQUID INFLOW OCCURS. FT 

HAVL= HEIGHT WHERE VAPOR FLOWS ACROSS SCREEN SURROUNDED BY LIQUID. FT 
HAVW= HEIGHT WHERE VAPOR FLOW OCCURS ACROSS WETTED SCREEN. FT 
NITER =ITERATIOH COUNTER 
COMMON/DVCOL/AST . A2 . B2 > RANGE 

COMMON/FLOW/AL 1 . AL2 . AVI . AV2 . CONVERG. CONVl . C0NV2 . DPINl . DPIN2 . HOL 1 . 

1H0L2. NITER 

COMMON/IMPINGE/A. AREA(20).ATOT.B.NDR.NIM.NW.VOLO(20).Y.Z 
COMMON/IMPOUT/CVOL , DPS . HDPIN. HPULL . VOL 

COMMON/INPUT/DT.G,HT(20),HTOT,I.KC.NOUT.NST.OA.QP.RHOL.SI(5MA. 

1VOLK20) 

C0MM0N/REFILL/AL.AV.AVW,DPIN,DPS2.HI.HIL1.HIL2.H0.HS.HX.JUMP.QL. 
IQO.T.VCL.VIL.VILMAX.VO.VOLOUT.VOREF.VV.WT 
COMMON/ST AND/ AVLR, AVR , AVWR 

COMMGN/WRIT/ADL.AVL.AVLS.AVS.AVT.AVWS.HDYN.VL.VT.VVL.VVLS.VVS.VVW. 

IVVWS 
REAL KC 
WRITE(6,5) 

5 FORMATOf SUBROUTINE FLOW^O 

NITER=0 

IF(NOUT.NE.0)DPIN=0.001 

IF(DPIN.EQ.O. )OPIN=0.5K(HO-HI)/32.2«G»RHOL 
HDPIN=DPIN‘<32.2/(RH0L)(G) 

C AREA AND LIQUID LEVEL CALCULATIONS xxxxxxxxxxxiotxxxxxxxxxioeNXKXXKKMHBlO/TT 
10 CALL TABL(WT,AVWT,HT(1),AREA(1).1,1.1,I.IFBAD) 

CALL TABL<H0,AVLT,HT(1),AREA(1),a,1,1,I,IFBAD) 

CALL TABL(HI,AVT.HTa),AREA(l).l,l.l.I.IFBAD) 

AV=AT0T-AVWT 
AVW=AVWT-AVLT 
20 CONTINUE 
HAVL=PI+HS 

IF(HAVL.GT.H0)HAVL=H0 

HAL=HO-HDPIN 

CALL TABL(HAVL.AVL.HT(I),AREA(1).1.1.1.I.IFBAD) 

CALL TABL(HAL,AL.HT(1).AREA(1).1.1.1.I.IFBAD) 

al=al-avt 

AVL=AVLT-AVL 

1F(HS.GE.HDPIN)AVL=0. 

IFTAL.LT.O. )AL=0. 

HDP1N=DPIN><32.2/(RH0L«G) 

IFTDPIN.LE.O. )GO TO ?1 
IFCAV.EQ.O. )G0 to 21 

IF(NST.EQ.O)GO TO 80 MHBlO/77 

C VVS =VAPOR VELOCITY FOR UNWETTED STANDPIPE NXXiHnmxxKlixMHBl 0/77 

VVS=(-A2+SQRT(A2XX2+A.XB2X(DPIN)))/(2.XB2) • MHBlO/77 

80 CONTINUE MHBlO/77 

C VAPOR FLOW ACROSS DRY SCREEN. OUT OF BASK'"'^ XNXX>nn»*»»*NNNXN>tJn»NX}<KMHB10/77 
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VV:(-A^SQRT(AXMm.MBN(DPIN)))/(2.»B> 

GO TO 22 

21 VV:0. 

vvs=o. 

22 CONTINUE 

IF(NST.EQ.O)GO TO BS HHBlO/77 

IF(DPIN-DPS2.LE.O.)GO TO 82 

C VVWS=VAPOR VELOCITY FOR WETTED STANDPIPE iiANin<)OMMn«xxMHB10/77 

VVUS=(-A2«SQRT(A2K)(2^A.mB 2K(DPIN-DPS2)))/(2.MB2) MHBlO/77 

C VVLS=VAPOR VELOCITY FOR STANDPIPE SURROUNDED BY LIQUID xkkkkkmkmHBIO/TT 

VVLSsI-AR^SQRKAEXKEI-R.KBRKCDPIN-OPSRDI/LR.kBR) nHBlO/77 

GO TO 83 riHBlO/77 

82 VVUS=0. 

84 VVLS=0. 

83 CONTINUE nHBlO/77 

IF(DPIN-DPS.LE.O.)GO TO 25 

C VAPOR FLOW ACROSS WETTED SCREEN INTO VAPOR. OUT OF BASKETKKXMXKXXXxriHBlO/TT 
VVW=(-A+SQRT(A»<»2*4.ABX(DPIN-0PS)))/<2.XB) 

IF(AVW.EQ.0.)VVW=O. 

C VAPOR FLOW INTO LIQUID, OUT OF BASKET KKXXKKXKMXKMKXxxionntintKKKKKiiKnHBlO/TT 

VVL=<“A*SQRTCAXX2+2.«BX(DPIN-0PS)))/(2.NB) 

GO TO 26 

25 VVU=0. 

VVL=0. 

26 CONTINUE 

IF(HO-HI-HDPIN.LE.O.)GO TO 28 

C LIQUID FLOW INTO VAPOR, INTO BASKET xxxkxxxxxxmkkkkxxxxkkxxxxxxxxxxI1HB10/77 

VL=(-Y+SQRT(Yxx2+2.XZ»((H0-HI-HDPIN)xRH0LXG/32.2))/(2.x2) 

C LIQUID FLOW INTO LIQUID, INTO BASKET xxxxxkxxkxxxkxxxxxxxkxxxxxxxxxKHB10/77 

VT=(-Y+SQRT<Yxx2+6.x2X(H0-HI-HDPIN)xrH0LXG/32.2))/(2.x2> 

GO TO 29 

28 VL=0. 

VT=0. 

29 CONTINUE 

C CONVERGENCE CALCULATIONS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


DPIN1=DPIN 

IF(NST.EQ.0)G3 TO 93 nHBlO/77 

CALL STAND 

CONVERG=AVRKVV4^AVSKVVS4^AVWRXVVW>AVWSXVVWS^AVLRNVVL+AVLSKVVLS-ALxVLMHB10/77 
1-AVTXVT MHBlO/77 

QV0=AVRKVV+AVSKVVS«AVWRXVVW4^AVWSXVVWS-»AVLSXVVLS4^AVLRXVVL 

GO TO 99 MHBlO/77 

93 CONTINUE MHBlO/77 

CONVERG=AVxyV+AVWXVVW+AVLXVVL-ALXVL-VTXAVT 
QV0=AVXVV^AVWXVVW4AVLXVVL 

99 CONTINUE HHBlO/77 


Q0 = 0. 

IF(NOUT.NE.O)CALL OUTFLOW 
qin=alxvl«oa+oaxvtxavt 

QL=QIN 

IF(NITER.GT.l)GO TO 72 
IF(QO.EQ.O.)GO TO 72 

C DETERMINES IF NET BASKET REFILL EXISTS xxxxxxxxxxxxxxxxxxxxxxxxxxMHBlO/77 
ir(QO.GT.QIN)GO TO 71 

C NO VAPOR OUTFLOW xxxxxxxxxxxkxkkxkxxxkxxxxxxxxxxxxxxxxxxxxxxxxxxxkmhb10/77 
72 CONTINUE 

C0NVER6=C0NVER6X0A4-Q0 

CONVl=CONVERG 

IF(ABS(CONVERG).LE.RANGE)GO TO 40 
IF(NITER.EQ.O)GO TO 30 
IF(WT,GE.HTOT.A.HX.GE.HO-HI)GO TO 37 
IF(NITER.GE.75)G0 TO 300 
IF(C0NV1.EQ.C0NV2)60 TO 30 
NITER=RITER+1 

IF(NITER.GT.10)CALL DIAGNOS 
DPIN=(DPIN2XC0NV1-DPIN1XC0NV2)/(C0NV1-C0NV2) 

HDPIN=DPINX32.2/CRH0LXG) 

IF(NW.EQ.2.AND.NDR.NE.l) GO TO 43 
IF(WT.GE.HTOT.AND.NDR.NE.l) GO TO 43 
IF(AV.NE.0.)6O TO 41 
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43 CONTINUE 

IF(HOP1N.IE.HX)HDPIN=HX4^(HQ-HI*HX)/(NITERM100000) 

IF(HDPlN.6T.(HO-Hl))HDPINs(HO-Hl)-(HO-HI-HX)/(NlTERxiOOOO) 

CO TO 42 

41 IF(DPIN.LE.0.)HDPIN:1./(NITERK10) 

42 CONTINUE 

IF(NITER.EQ.25.AND.ABS(HX-HDPIN).LT.HX/100) HDP1N:2MHDPIN 
0PIN=HDP1NKRH0LKG/32.2 
C0NV2:C0NV1 
0P1N2=DP1H1 
CO TO 20 
30 DP!N2=DPIN1 
C0NV2=C0NV1 
NITER=NITER+1 

IF(NITER.CT.10)CAU DIAGNOS 
DPlN=DPIN/2. 

HDPIN=DPIN»32.2/(RH0LXG) 

IF(NITER.EQ.1.A.UT.GE.HTOT.A.HX.NE.O.)HDPIN:HX 
DPIN=HDPIN»RH0L><G/32.2 
GO TO 20 
40 CONTINUE 

C KKX OUTFLOM RATE CALCULATION KKintxxKiOoniKMMKKXMMXKMXKinnnoMiKXXKK 
60 QL=ALXVL>«0A+AVTXVTX0A 
V0REF=QLXDT 

IF(NITER.GT.10)CALL DIAGNOS 
IF(VOREF.LT.VILMAX)CO TO 77 
GO TO 270 
77 CALL DVCOL 
GO TO 270 

37 IF(NOUT.NE.O)60 TO 71 
URITE(6>38) 

38 FO'^MAT(x HS IS GREATER THAN HO-HI, NO SOLUTION FOR HDPIHX) 

’ JREF=0. 

IF((HC.LT.HTOT).A.(VOL+VIL.LT.CVOL))GO TO 270 
IF(VOLOUT.NE.O. )GO TO 270 

GO TO JUMP. TRANSFER CONTROL TO MAIN PROGRAM-STATEMENT 10, IN ORDERMHBlO/77 
TO RESTART A NEU CASE xxxxxxxxxxxxxkxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxmHB10/77 
GO TO JUMP, (20) 

300 URITE(6,301) 

301 F0RMAT(x TOO MANY ITERATIONS IN SUBROUTINE FLOW X) 

GO TO JUMP, (20) 

71 CALL IMPOUT 
270 RETURN 
END 
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SUBROUTINE IMPOUT 


Subroutine IMPOUT performs flow oaloulations when liquid levels inside and outside the 
basket are unequal after a time step and net outflow from the basket ooours* (The basket 
is being drained rather than refilled. ) A separate subroutine is needed for net outflow 
because the pressure differences inside and outside the basket are reversed compared to 
that when net refilling is occurring. Pressure of the vapor inside the basket is greater 
than the pressure outside the basket. Dynamic pressure (screen Impingement) effects 
are not included. Standpipe areas are included in the calculations. Convergence 
calculations are performed. 

SUBROUTINE IMPOUT 

ROUTINE CALCULATES BASKET REFILLING. SIMILAR TO OVCOL- LIQUID LEVELMHBlO/77 
INSIDE AND OUTSIDE THE BASKET UNEQUAL AFTER A TIME STEP - NET MHBlO/77 
OUTFLOW OCCURS. CREATING VAPOR INFLOW TO REPLACE LIQUID LEAVING xkkmHBIO/ 77 
DYNAMIC PRESSURE USED TO CALCULATE LIQUID AND VAPOR FLOW KKXK)Oon(KMHD10/77 
FOR IMPOUT A POSITIVE VALUE OF DPIN AND HDPIN MEANS THAT CAS PRESSURE 
INSIDE THE BASKET IS LESS THAN GAS PRESSURE OUTSIDE THE BASKET-THIS 
IS OPPOSITE TO THE SIGN CONVENTION USED IN DVCOL AND FLOW. 

DVCOL=SUBROUTINE TO COMPUTE FINAL HEAD WHEN VOREF.LT.VILMAX 
AVI = INITIAL UNWETTED VAPOR FLOW AREA.FT2 
AV2 = FINAL UNWETTED VAPOR FLOW AREA.FT2 

C0NV1= CONVERGENCE VARIABLE FOR CONTINUITY EQUATION-CURRENT ITERATION 

C0NV2= CONVERGENCE VARIABLE FOR CONTINUITY EQUATION-LAST ITERATION 

CONVERGE CONVERGENCE VARIABLE THAT MUST BE WITHIN ZERO 4-OR-RANGE 

DPIN1= STORAGE FOR CURRENT DELTA P ITERATION 

DPIN2= STORAGE FOR CURRENT DELTA P ITERATION 

DPIN = PRESSURE DIFFERENCE OF GAS INSIDE AND OUTSIDE BASKET 

MALI = INITIAL DISTANCE FOR LIQUID INFLOW, FT 

HAL2 = FINAL DISTANCE FOR LIQUID INFLOW, FT 

HAVL1= INITIAL DISTANCE FOR VAPOR FLOW INTO LIQUID, FT 

HAVL2= FINAL DISTANCE FOR VAPOR FLOW INTO LIQUID, FT 

HAVW1= INITIAL DISTANCE FOR VAPOR FLOW ACROSS WETTED SCREEN, FT 

HAVW2= FINAL DISTANCE FOR VAPOR FLOW ACROSS WETTED SCREEN, FT 

NITER =ITERATION COUNTER 

VORl= STORAGE FOR VOREF 

UT1,WT2=WETTIN6 PARAMETERS AT START AND END OF TIME STEP 
COMNON/DVCOL/AST , A2 , B2 , RANGE 

COr.MON/FLOW/ALl,AL2,AVl,AV2,CONVERG,CONVl,CONV2,DPINl,DPlN2,HOLl. 

1H0L2, NITER 

COMMON/IMPINGE/A, AREA(20>,ATOT,B,NDR,NIM,NW,VOLO(20),Y,Z 
COMMON/IMPOUT/CVOL, DPS, HDPIN, HPULL, VOL 

COMMON/INPUT/DT,G,HT(20),HTOT,I.KC,NOUT,NST,OA,QP,RHOL,SIGMA, 

1VOLK20) 

COMMON/REFILL/AL,AV,AVW,DPIN,DPS2,HI,HIL1,HIL2.HO,HS,HX,JUMP,QL, 
1Q0,T,VCL,VIL, VILMAX, VO, VOLOUT, VOREF, VV.WT 
COMMON/STAND/ AVLR , AVR , AVWR 

CONNON/WRIT/ADL,AVL,AVLS,AVS,AVT,AVWS,HDYN,VL,VT,VVL,VVLS,VVS,VVW, 

IVVWS 

COMNON/INFILL/CVI(20),DBP,DBP2,HCOL(20),HPT(10),J,K,MPROP,MR,MVEH. 

lN,NG,NS,NSP,QOUTa0),TC(20),THRST,TIM(20),TITLE(7),VILI,VTOT(20) 

REAL KC,MPROP,MR,MVEH 
WRITE(6,5) 

5 F0RMAT(« SUBROUTINE IMPOUT *) 

NITER=0 
DPIN=0.001 
10 VOR1=VOREF 

IF(NOUT.EQ.O)VOLOUT=0. 

VILX2=VIL+VOREF-VOLOUT 

V0LX2=V0L-V0REF 

200 CALL TADL(VIL,HILl,VOLI(l),HT(l),l,l,l,I,IFBAD) 

CALL T ABL ( VOL , HOL 1 . VOLO< 1) , HT( 1 ) , 1 , 1 , 1 , I , IFBAD) 

CALL T ABL ( VILX2, HIL2, VOL I(1),HT(1), 1,1, 1,1, IFBAD) 

CALL TABL(QP,HPULL,Q0UT(1),HPT(1),1,1,1,J, IFBAD) 
IF((NOUT.NE.O).A.(HIL2.LE.HPULL))GO TO 302 
CALL TABL(V0LX2,H0L2,V0L0<1),HT<1), 1,1, 1,1, IFBAD) 
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C DVNANIC PRESSURE INCLUDED IN PLOW CALCULATIONS KMNMMMKMKMMKioniXMNMnHBlO/?? 
DPDYN=AKV0*B»«V0KX2 
HDYNsDPDYNxSZ.R/IRHOLXC) 

IF(HIL2.OE.HT0T)HIL2sHT0T 

IF(HOL2.0E.HTOT)HOL2=NTOT 

IF(HOLl.GE.HTOT)HOLlsHTOT 

C AREA AND LIQUID LEVEL CALCULATIONS KMNMMMNKMKMKMMMKKKMMKiiMMMXNMXKKnHBlO/?? 
20 IF(DPIN.EQ.O.)DPIN=HS/(NITER*U 
HDP1N=DPIN«32.2/(RH0LJ«6) 

HAL1:HOL1«HOPIN 

NAL2:H0L2^HDPIN 

IF(HDPIN.6T.O.)HAL1:HOL1 

IF(HDPIN.GT.O.)HAL2:HOL2 

CALL TABL(HALl.AVLTl.HT(l),AREA(l).ia»l»I*lFBAD) 

CALL TABL(HAL2.AVLT2,HTa).AREA(l).l.l.l»I.lFBAD) 

CALL TABL(HOLl,AVLJl>HT(l),AREA(l)a*ia>l>IFBAD) 

CALL TABL(H0L2,AVLJ2.HT(1),AREA(1).1»1.1»I.IFBAD) 

CALL TABL(HlLl»AVTl,HTa)»AREA(l),l.l»l*I>IFBAD) 

CALL TABL(HIL2»AVT2.HTa).AREA(l).l.l»l.I.IFBAD) 

AVT=(AVTl+AVT2)/2. 

H0=H0L2 

1F(H0.GT.HT0T)H0=HT0T 

IF(NU.EQ.1)UT:M0 

IF(NW.EQ.2)WTsHT0T 

if(nn.eq.3)call suet 

UT2=WT 

HO^HOLl 

IF(HO.GT.HTOT)HO=HTOT 
IF(NW.EQ.l)WTsHO 
IF(NW.EQ.2)UT=HT0T 
IF(NU.EQ.3)CALL SUET 
UT1=UT 

CALL TABL(UT1,AVUT1,HT(1).AREA(1)»1»1»1»I.IFBAD) 

CALL TABL(UT2.AVUT2.HT(1)>AREA(1),1>1.1>I.IFBAD) 

AV1=AT0T-AVUT1 

AV2=AT0T-AVUT2 

AVW1=AVUT1-AVLJ1 

AVU2=AVUT2-AVLJ2 

AL1=AVLT1-AVT1 

IF(ALl.LT.O) AL1=0 

AL2=AVLT2-AVT2 

IF(AL2.LT.O) AL2=0 

HDL1=H1LI-HDP1N-HDYN 

HDL2=HIL2-HDPIN-HDYN 

1F(HDL1.LT.0.)HDL1=0. 

IF(HDL2.LT.0.)HDL2=0. 

CALL TABL(HILl.AVLl>HT(l),AREAa^.l»I*l>I>IFBAD) 

CALL TABL(HIL2,AVL2.HT(l).AREA(l),iaa»I»IFBAD) 

CALL TABL(HDLl.ADLl.HT(l).AREA(l).ia.l>I.IFBAD) 

CALL TABL(HDL2»ADL2.HT(l),AREA(l)>l.ia>I»lFBAD} 

ADL3=AVL1-ADL1 

ADLA=AVL2-ADL2 

ADL1=AVLT1-ADL1 

ADL2=AVLT2-ADL2 

IF(ADL1.LT.O.)ADL1=0. 

IF(ADL2.LT.O. )ADL2=0. 

ADL = ( ADL1+ADL2 )/2 . + ( ADL3+ADLA )/A . 

AV=(AVl+AV2)/2. 

AVU=(AVUl+AVU2)/2. 

AL =(ALl+AL2)/2. 

H0PIN:DPINX32.2/(RH0LKG) 

IF(DPIN.LE.O.)GO TO 21 

C VAPOR FLOU ACROSS DRY SCREEN, OUT OF BASKET 

IF(NST.EQ.O)GO TO BO MHBlO/77 

C VVS =VAPOR VELOCITY FOR UNUETTED STANDPIPE XHKimnnimximMHBlO/77 

VVS=(-A2+SQRT(A2«»2+A.XB2K(DPIN)))/(2.»B2) MHBlO/77 

\/vS=-VVS MHBlO/77 

80 CONTINUE MHBlO/77 

VV=(-A+SQRT(A«X2*A.XB«(DPIN. )/(2.xB) 

VV=-VV 
GO TO 22 
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c 

c 

c 

c 

c 


21 VV^O. 

vvs*o. 

22 CONTINUE 
IF(NST.EQ.O)GO TO 83 
IF(DPIN>DPS2.LE.O.)60 TO 82 
VVWS:VAFOR VELOCITY FOR WETTED STANDPIPE 
VVMSs(-A2+SQRT(A2«»<2+4.»<B2>«(DPlN-DPS2)))/(2.J<B2) 

VVWSs-VVWS 
GO TO 83 

82 VVWS=0. 

83 CONTINUE 
IF(DPIN-DPS.LE.O.)60 TO 25 

VAPOR FLOW ACROSS WETTED SCREEN INTO VAPOR, OUT OF BASKETkxmkkmmkmnmHBIO/?? 
VVW=(-A+S0RT(A«»2*^.»»BN(DPIN-DPS)))/(2.«B> 

VVW=-VVW 
GO TO 26 


MHBlO/77 

KMKMMNKNKNKKtlHBlO/77 

MHBlO/77 

MHBlO/77 

MHBlO/77 

MHBlO/77 

MHBlO/77 


25 VVW=0. 

26 CONTINUE 
HOPINsR.KNDPIN 

IF(HOL1-HIL1+HOL2-HIL2^HOPIN.LE.O.)GO TO 28 

LIQUID FLOW INTO VAPOR, INTO BASKET MKitKXMXKKKKKXiOiNKKKKXKKKKXinontKnHBlO/TT 
VL = (-Y+SQRT(YXX2+ 2K(H0H-HIL1*H0L2-HIL2+2 .«HDPIN)xRH0L»«G/32.2)) 

1/(2. XZ) 

LIQUID FLOW INTO LIQUID, INTO BASKET xxxxxxxxxxxxxxxxxxxxxxxxkxxxkxmHB10/77 
VT=(-Y+SQRT<Yxx2*2.x2x<H0H-NIL1+H0L2-HIL2+2.xHDPIN)xRH0LXG/32.2)) 

1/(2. XZ) 

GO TO 29 
28 VL=0. 


VTso. 


KKXMXKXMXXKXXXXXXXXXXXKXXXXMXXKXXXXX 


29 CONTINUE 

xxxx CONVERGENCE CALCULATIONS 
DPIN1=DP1N 

IF(NST.EQ.O)GO TO 93 
HI=(HILl+HIL2)/2. 

CALL STAND 

converg=avrxvv+avsxvvs+avwrxvvu*avwsxvvws-alxvl-adlxvo+qo-vtxavt 

GO TO 99 
93 CONTINUE 

converg=avxvv+avwxvvw-alxvl-adlxvo+qo-vtxavt 

99 CONTINUE 

QL=(AL1+AL2)/2.XVLXOA+(ADL1+ADL2)/2.XVOXOA*AVTXVTXOA 

CONVERG=(CONVERG-QO)XOA+QO 

CONV1=CONVERG 

IF(CONVER6.EQ.O..A.VL.EQ.O.)GO TO 35 
IF(ABS(CONVERG).LT.RANGE)GO TO 60 


MHBlO/77 


MHBlO/7 

MMBlO/77 

MHBlO/77 

MHBlO/77 


35 CONTINUE 

IF(NITER.EQ.O)GO TO 30 
IF(NITER.GE.75)GO TO 300 
IF(CONV1.EQ.CONV2)GO TO 30 
NITER=HITER+1 

IF(NITER.CT.10)CALL DIA6N0S 

DPIN=(DPIN2XCONV1-DPIN1XCONV2)/(CONV1-CONV2) 

C0NV2=C0NV1 

DPIN2=DPIN1 

DX=DPS2 

IF(NST.EQ.O)DX=DPS 

IF(DPIN.LT.O. .A.DPIN.LT.-DX)DPIN=-DX+DX/(NITERX10) 
HDPIN=DP1NX32.2/RH0L/6 
GO TO 20 
30 DPIN2=DPIN1 
C0NV2=C0NV1 
NITER=NITER+1 

IF(NITER.GT.10)CALL DIAGNOS 
DPIN=DPIN/2. 

HDPIN=DPINX32.2/(RH0LXG) 

GO TO 20. 

40 CONTINUE 
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C OUTFLOW RATE CALCULATION AKKNItAXKAAAKXAAUAKAAHAKKKmiAllHAFUKII 

QLsCALl+ALaT/Z.AVLuOA+CADLl+ADLRI/R.xVO^OA+AVTAVTAOA 
V0REF=QLADT 

IF(NOUT.HE.O)CALL OUTFLOW 
IF(ABS(V0REF-V0RI).LE.RAHGE)00 TO 72 
VOREF=(VOREF+VORl)/2. 

IF(NITER.GT.10)CALL DIAGNOS 
GO TO 10 

300 URITE(6.301) 

301 FORMAT(A TOO MANY ITERATIONS IN SUBROUTINE IMPOUTA) 

GO TO JUMP >< 20) 

502 MRITE(6,303)I1PULL _ , 

303 FORMAT(A LIQUID LEVEL BELOW PULITHROUGH HEIGHT, HT=AF9. 5) 

GO TO JUMP. (20) 

72 RETURN 
END 



SUBROUTINE INPUT 


All Inputs are defined and formatted in Subroutine INPUT. Namelist input is used. 
Comment cards define the inputs. 

SUBROUTINE INPUT 
C INPUT QUANTITIES 
C 

C ALL TABLES MUST BE MONOTONIC FUNCTIONS kkmxmkkkknknknknmknknknkkk NHB3/7S 

C A : VISCOUS CONSTANT IN PRESSURE LOSS EQUATION FOR VAPOR 

C AREA(I): BASKET SURFACE AREA BELOU HT(1).FT2 

C AST« STANDPIPE SURFACE AREA MHBlO/77 

C ATOTs total START BASKET AREA.FT2 

C AWIDTH=WIDTH OF SECTION, FT 

C A2:VISC0US CONSTANT FOR TOP SCREEN PRESSURE LOSS 

C B = INERTIAL CONSTANT IN PRESSURE LOSS EXPRESSION FOR VAPOR 

C B2=INERTIAL CONSTANT FOR TOP SCREEN PRESSURE LOSS 

C CC^CORRECTION FACTOR FOR ORIFICE ROUNDED ENTRANCE, MAD PIG 3.22 

C CVl(N) = IMPINGEMENT VELOCITY FOR TIM(N)VS CVKN) TABLE, FT/SEC 

C D: spacing BETWEEN SCREEN BARRIERS, FT 
C DBP: SCREEN BUBBLE POINT DIAMETER, MICRONS 

C DBP2 = BUBBLE POINT FOR TOP SCREEN, MICRONS 
C DORsOUTLET DIAMETER, FT 

C DPRESSs TANK PRESSURE DIFFERENCE ABOVE INTERFACE VAPOR PRESSURE, PSF 
C DPX s PRESSURE CORRECTION TERM FOR WICKING,PSF 
C DT = TIME STEP FOR FLOW CALCULATIONS, SEC 

C ENsNUMBER OF HOLES IN DIRECTION PERPENDICULAR TO UICKING FLOW 

C G=AMBIENT ACCELERATION FOR CONSTANT G CASE, FT/SEC2 
C HCOL(N)s COLLECTED HEIGHT FOR TC(N) VS HCOL(N) TABLE, FT 

C HPT(N) PULLTHR0U6H HEIGHT FOR QOUT(N) VS HPT(N) TABLE, FT 

C HT(I)= BASKET HEIGHT CORRESPONDING TO VOLI(I),VOLO(I),VTOT(I) 

C HTOL=UNFILLED HEIGHT FRACTION 

C HTOT s MAXIMUM HEIGHT OF BASKET, FT 

C I=NUMBER OF ENTRIES IN INPUT TABLES. HT(I) ,VOLl(I) ,VOLO(I) ,VTOT(l) ,AREA(I) 
C J = NUMBER OF ENTRIES IN QOUT VS HPT TABLE 
C K = INDEX OF MAXIMUM HCOL VALUE 

C KC=CORRECTION FACTOR FOR ABRUPT CONTRACTION, WAD FIG. 3. 21 
C MPROPsPROPELLANT MASS, LBM 

C MRsMIXTURE RATIO, MR=6 FOR LH2, MR=l/6 FOR L02 
C MUL = LIQUID VISCOSITY, LBM/FT-SEC 

C MVEHsVEHICLE MASS, LBM 

C NsNUMBER OF ENTRIES IN INPUT TABLES TC(N),AND CVOL(N) 

C NDRsSTANDPIPE SCREEN WETTING FLAG, NDRsOi SCREEN WET, NDR.NE.O* SCREEN DRY 

C NG^G LEVEL FLAG, NG=0 FOR VARIABLE G CASE, NG=I FOR CONSTANT G CASE 

C NIMsFLAG,IF.NE.O. IMPIMGEMENT TABLES ARE INPUT AND IMPINGEMENT CALCULATED 

C NOUT=FLAG FOR OUTFLOW, NO OUTFLOW IF NOUT sO 
C NS= NUMBER OF SECTIONS OF SCREEN 

C NSP s SPILLING FLAG,NSP=0-NO SPILLING. NSP.NE.O, SPILLING 
C IF SPILLING DOES NOT OCCUR INITIAL BASKET LEVEL MUST BE LESS THAN TANK 

C H5T=FLAG FOR STANDPIPE, NS=0-NO STANDPIPE MHBlO/77 

C NTSCRs TYPE OF SCREEN BARRIER, NTSCR=l ,P/S-S/P,NSCRs2,P/S-P/S 
C NW s flag to determine METHOD OF COMPUTING SCREEN WETTING 
C NW:1,WT:H0 NW=2,WT=HT0T NUs3,WT IS COMPUTED ACCORDING TO UICKING 
C OA = SCREEN SURFACE OPEN AREA FRACTION 

C PsPOROSITY, SCREEN VOID FRACTION IN THE DIRECTION OF WICKING FLOW 
C PTs plate THICKNESS, ft 

C QOUT(N): FLOW RATE USED IN OUTFLOW RATE VS PULLTHROUGH HEIGHT, FT3/SEC 

C QP = CONSTANT OUTFLOW RATELUSED WHEN N0UT=5) , FT3/SEC 
C RANGE = ITERATION TOLERANCE 

C RHOLs LIQUID DENSITY, LBM/FT3 

C SIGMAS SURFACE TENSION, LBF/FT 

C ST s SCREEN THICKNESS, FT 

C TC(N)sCOLLECTION TIME AT HCOL(H),SEC 

C THETAS AH6LE BETWEEN SCREEN BARRIER AND HORIZONTAL 
C THRSTSVEHICLE THRUST, LBF 

C TlM(N)s TIME AT IMPINGEMENT , CVKN) , SEC 
C TMAX s maximum TIME FOR CALCULATIONS, SEC 
C VILIsINITIAL LEVEL IN START BASKET, FT3 

C VOLKDsVOLUME INSIDE BASKET AT HT(I), FT3 
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V0L0(1)SV0LUHE OUTSIDE THE BASKET AT HT(I), FTS 
VT0T(I)s total tank VOLUME AT LEVEL HT<I) 

Y sVISCOUS CONSTANT FOR LIQUID FLOU ACROSS SCREENS 
Z ^INERTIAL CONSTANT FOR LIQUID FLOM ACROSS SCREEN 
COMMON/DVCOk/AST , A2 » D2 , RANGE 

COMMON/IMPINGE/A»AREA(ZO)»ATOT»B.NDR.NIM.NU»VOLO(ZO)>Y,Z 
COMM0N/INFILL/CVI(20>,OBP.D&P2.HCOL(2O),HPT(10)»J,K.MPROP»MR»nVEH. 
lM,NG,NS.NSP.QOUTaO),TC(20)»THRST,TIM(20),TITLE(7),VILl,VTOTt2Q) 
COMMON/INPUT/DT , 6. NT < 20 ) , HTOT » I . KC > NOUT , NST » OA , QP » RHOL » SIGMA , 
1VOLK20) 

COMMON/OUTFLOU/CC . DOR . DPRESS 

COMMON/SUET/AUIDTH.D.DPX»EN.MUL>NTSCR»P.PT.ST»THETA 
COMMON/T IHE/HTOL , TMAX 
REAL KC.MPROP.NR.nUL.MVEH 

NAMELIST/INDATA/A, AREA, AST. AT0T»AUI0TH>A2,B»B2»CC,CVI»D,DBP»DBP2. 
lOOR , OPRESS , OPX , DT , EN , G . HCOL , HPT , NT , HTOL , HTOT , I , J , K , KC .MPROP .MR, 
2MUL.MVEH,N,NDR,NG.NIM,N0UT,NS,NSP,NST.NTSCR,NW,0A,P,PT,Q0UT.QP, 
SRANGE. RHOL .SIGMA. ST, TC. THETA, THRST. TIM. TITLE. TMAX, VILI.VOLI.VOLO. 
AVTOT.Y.Z 
REA0(5.IN0ATA) 

IF(E0F(5>) 10.20 
10 CALL EXIT 
20 URITE(6. INDATA) 

RETURN 

END 
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SUBROUTINE OUTFLOW 


Subroutine OUTFLOW oalculatee the outflow rate from the start basket based on tank 
pressure, outlet pressure, basket geometry and outlet geometry. 


SUBROUTINE OUTFLOU 

C CALCULATES OUTFLOU RATE BASED ON BASKET PRESSURE AND OUTLET GEOnETRYXKKMK 
COr.MON/INPUT/DT,G,HT(20),HTOT,I,KC,NOUT,NST,OA,OP,RHOL,SIGMA, 

1VOLK20) 

COMMON/OUTFLOW/CC, DOR , DPRESS 

C0m0N/REFILL/AL.AV,AVU,DPINr0PS2,HI,HILl,HIL2»H0.HS,HX, JUMP.QL, 
lQO,T,VCL,VIL,VILHAX.VO.VOLOUr.VOREF,VV,UT 
REAL KC 

IF(NCUT.EG.5)60 TO 20 

DPF=RH0L»«G/32.2>«(HlLl*HlL2)/2.+DPRESS-DPIN 
U'RITE(6,30)HIL1.HIL2.DPRESS,DPIN,DPF 
30 FORMATOt HILl = >iF10 .3, X HIL2 = HFin . 3,x DPRESS=«Fl0.3,x DPIN=xF10.5, 
lx DPF=xF10.3) 

1F(DPF.LE.0. )GO TO 10 

,C DETERMINES START BASKET OUTFLOU RATE xmxmm»mkmnmmmmmkmmmm»mmkmmmmmmHB10/77 
VOUT=5QRT(2.x32.2XDPF/(KCXCCXRHOl)) 

GO TO 1> 

10 D?F=0. 

VOUT=0. 

15 CONTINUE 

Q0=(V0UTX3.1A159xD0Rxx2)/A. 

GO TO 21 

C CONSTANT OUTFLOU RATE - INPUT XMMKKMMKMMKMMKKXi(«(MK«(M)OO(MM)OOOO()O(MMMHB10/77 

20 Q0=QP 

21 CONTINUE 
VOLOUT=QOxDT 
VILMAX-VILMAX+VOLOUT 
RETURN 

END 
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SUBROUTINE STAND 


Subroutine STAND is used when a standpipe on top of the basket is employed. Standpipe 
flow areas are computed. The &reas depend on screen wetting, basket pressure 
differences and liquid levels inside and outside the basket. Calculations determine 
vapor flow area of the standpipe for flow into or out of the start basket. Vapor flow 
area for flow across dry screen, wet screen and screen surrounded by liquid is 
computed. 


(A flag is used to determine whether the screen is wet or dry independent of the 
calculations of SWET or flags for NW. This was used to simulate conditions that 
occurred during testing. ) 

SUBROUTINE STAND 

STAND COMPUTES STANDPIPE AREA FOR VAPOR AND LIQUID FLOU mmmmnkmmm 
AVR =AREA FOR UNUETTED VAPOR FLOW OUTSIDE STANDPIPE KUKMXMKtmxxwMHBlO/?? 
AV5 =AREA FOR UKMETTED VAPOR FLOW ACROSS STANDPIPE »on<«>i)noo*xx»(MHB10/77 
AVWRsAREA FOR WETTED VAPOR FLOW OUTSIDE STANDPIPE mmnnnnknkknmMHB10/77 
AVWS = AREA FOR WETTED VAPOR FLOW ACROSS STANDPIPE X»»<««»O»«»«««»*MIIB10/77 
AVLR = AREA FOR VAPOR FLOW INTO LIQUID. OUTSIDE STANDPIPE»‘»>“‘»>»K«x»»»'<MHBlO/77 
AVLS=AREA FOR VAPOR FLOW INTO LIQUID ACROSS STANDPIPE«xxxxxxxxxxxMHB10/77 
COMMON/IMPOUT/CVOL.DPS.HDPIN.HPULL.VOL 

COMMON/INPUT/DT,C.HT(20).HTOT,I,KC,NOUT,NST,OA,QP,RHCL,SIGMA, 

1VOLK20) 

C0MM0N/DVC0L/AST,A2,B2, RANGE 

COM;10H/IMPINGE/A.AREA(20).ATOT.B,NOR,NIM,NW.VOLO(20).Y.Z 
COMr*.ON/REFILL/AL.AV,AVW,DPIN,DPS2,HI,HILl.HIL2,H0.HS,HX, JUMP.QL. 

IQO.T .VCL.VIL.VILMAX.VO.VOLOUT.VOREF.VV.WT 
CCMMON/STAND/AVLR.AVR.AVWR 

COMMON/WRIT/ADL.AVL, WLS.AVS.AVT.AVWS.HDYN.VL.VT.VVL.VVLS.VVS.VVW, 

IVVWS 

HA:HI+HX 

CALL TADL(HA,AM,HT(1),AREA(1),1,1,1,I,IFBAD) 


If (AV.LE.AST)GO TO 

MHBlO/77 

1F(AV.CT.A5T)AVS=AST 

MHBlO/77 

AVR=AV-AST 

MHBlO/77 

Av;.;s=o. 

MHBlO/77 

AVLS=0. 

MHBlO/77 

AVl'R = AVW 
AVLR=AVL 

MHBlO/77 

GO TO 97 
AVWTO = AV-^AVU 

IF(AVWTO.LE.AST)GO TO 95 

MHBlO/77 

IF(AVWTO.GT.A$T)AVS=AV 

MHBlO/77 

AVR=0. 

MHBlO/77 

AVWS=AST-AVS 

MHBlO/77 

AVIS=0. 

MKBlO/77 

AVUR=AVW-AVWS 

MHBlO/77 

AVLR=AVL 

MHBlO/77 

GO TO 97 
AVWX=AV+AVW^AVl 
IF(AVUX,LE.A5T)G0 TO 96 

MHBlO/77 

IF(AVUIX.GT .ASTJAV^^AV 

MHBlO/77 

AVUS=AVU 

MHBlO/77 

AVR=0. 

MHBlO/77 

AVWR=0. 

MHBlO/77 

AVIS=AST-AVS-AVMS 

IF(HO.GE.HTOT.A.HX.LT.HS.A.HO.OT.HI+HX)AVLS=ATOT-AM 

IF(AVLS.GT.A5T)AVLS=AST 

MHBlO/77 

AVLR=AVL-AVLS 
CO TO 97 

MHB 10/77 
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9« AVSsAV MKBlO/77 

AVR=0. NHBlO/77 

AVUS^AVU 

AVWR=0. MHBlO/77 

I AVLS^AVL MHBlO/77 

' 1F(H0.66.HTOT.A.HX.IT,HS.A.HO.CT.HI+HX)AVLS=ATOT-AM 

ZF(AVLS.OT.AST)AVLS=AST ■ 

AVLRsO. MHBlO/77 

f 97 IF(NDR.EQ.l) AVS^AST 


RETURK 

END 


A-21 





i 


SUBROXJTINE SWEl 

Subroutine SWET computes wtoklng along the main screen using equations developed in 
''Centaur Propellant Thermal Conditioning Study/' by M. H. Blatt, R. L* Pleasant and 
R^ Erickson, NASA CR- 135032, CASD-NAS-76-026, NAS3-19693, July 1976. Plate/ 
screen-screen/plate and plate/screen«*plate/screen, the two ccmfiguratlons that were 
found to be most promising, are programmed for wicking between the inner barriers. 
An Iteration on time and wicking distance travelled is performed until a solution Is 
reached that satisfies the time step taken. 


The subroutine yields the distance wicked ahead of the level of liquid outside the start 
basket. (Options are available that do not use subroutine SWET to compute screen 
wetting. One option wets the entire screen instantaneously; the other option keeps 
the screen wetting at the liquid level with no wicking. ) 


c 


c 


c 


c 


c 


SUBROUTINE SWET 

CALCULATES SCREEN WETTING CONDITIONS AS A RESULT OF WICKING 
CONMON/INPUT/DT,6,HT(20),HTOT,I,KC,NOUT.NST,OA,QP,RHOL,SIGriA, 

1VOLK20) 

C0MM0N/REFILL/AL,AV,AVW,DPIN,DPS2,HI,HIL1.HIL2,HC,HS,HX,JU«P,QL, 
1Q0.T.VCL,VIL,\;ILMAX.V0»V0L0UT,V0REF,VV,WT 
COmON/SWET/ AWI DTH , D , DPX , EN , MUL , NTSCR , P , PT , ST , THETA 
REAL KC,MUL 
S1 = 0 
TI1=0. 

S=H0/5. 

AS=2)iSTJ<P 

PLATE/SCREEN-SCREEN/PLATE WICKING CALCULATION 

IF(NT5CR*EQ.2)G0 TO 10 

AF=D 

5 TI=AF/(AF-AS)'^:2<i.)<riULKSIGMA/RH0L«)(2/DKX3/GKX2x32-2/(SIN( THETA )XJ€2) 
lJ<(-RH0LJfG/32.2>^D>i5)^SIN( THETA )/2./SIGNA-AL0G(l.“RH0L«G/32.2^»DxS>«SIN 
l(ThETA)/2./SIGNA -DPX>«D/2 . )>^( 1 . -DPXxD/2 . ) ) 

GO TO 30 

PL ATE/SCREEN-PL ATE/SCREEN WICKING CALCULATIONS>^KX)^J(X^)t)(XX)(>()«JiXJ<i<i(XMHB10/77 
10 AF=D+PT«OA 

15 TIA = AF/(AF-AS))<192.^nULJ^SIGMA«(AMIDTH+EN)(PT))(«3 

TIB=(AWIDTHX5<3x(2.?<D+PT)i<K3XRH0L«»2>^6xJ<2^(SIN(THETA))K>J2/32.2) 

TIC=(-RH0LxG/32.2K(2.>JD+PT)«AWIDTH>''SIN(THETA)KS) 

TID = A.XSIGMA)f(AWIDTH + EN>iPT) 

TIP = (1.-DPX^AWIDTH?<(2.>^D+PT)) 

TIR = RH0L^<G><5'-SIN(THETA)J«AWIDTH)<(2.^<D+PT)/32.2 

IF(TIR.LE.TIP) go to 26 

S=S/10. 


00 TO 15 

26 1IU=AL0G((TIP-TIR)/TID) 

TI = (TIA)<TIB^<T1C-TIP^TIU)/TID 
30 WRITE(6,150)TI,S 

150 F0RriAT(K TI= ?^F12.4 ,k S = KF12.A) 

CONVERGENCE CALCULATIONS xxx«KXXX)<xXKXX)()O<xxKKXXi<K»KX)<)<XJ<M>tKXXKX)6KMHB10/77 

IF(TI~T.LT.T/10.)GO TO 100 

S2=S1 

S1=S 

TI2^:Tn 

TIUTI 

IF(TI2.NE.O. )G0 TO 60 
IFCTI.LT.O. )G0 TO 50 


S=Sl/2. 

IF(NTSCR.EQ.2) 15,5 


50 S = S19<2 

IFCNTSCR.EQ.2) 15,5 
60 SL0PE=(TI1-TI2)/(S1-S2) 

S=(TI1-T)/SL0PE+S1 
IF(NTSCR.EQ.2) 15,5 

END OF CONVERGENCE CALCULATIONS KX»xxxxxx>nixxxxxxx)oixxxxxxxxxxxxxxMHB10/77 


100 WT=H0+S 
RETURN 
END 
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SUBROUTINE TIME 


Subroutine TIME inorements the value of time and oheoks to see if the end time has been 
reached or the basket is Ml. 


SUBROUTINE TIME 

INCREMENTS TIME, DETERMINES UHETHER BASKET IS FULL JMtXKJomxxxKXXXXKMHBlO/?? 
DELT=TRUNCATED VALUE OF DT DUE TO BASKET FILLING 
VTIL2=V0LUME OF LIQUID IN THE BASKET IF FILLING IS UNCONSTRAINED 
COMMON/INPUT/OT,G,HT(20),HTOT,I,KC,NOUT,NST,OA,QP,RHOL,SIGMA. 

1V0LK20) 

CGMM0N/REFILL/AL,AV,AVU,DPIN,DPS2,HI,HIL1.HIL2,H0.HS.HX,JUMP,QL. 
1Q0.T,VCL,VIL,VILMAX,V0,V0L0UT,V0REF,VV,WT 
COMNON/TIME/HTOL , TMAX 
REAL KC 

IF(MT.EQ.HTOT) HMAX=HTOT-HX 

IF(UT.LT.HTOT)HMAX=HTOT 

VTIL2=VIL+VQREF-Q0XDT 

CALL TABL(VTIL2,HM.V0Htl),HTa>,l,l,l,I,IFBAD) 

IF(HM.LT.HmX)GO TO 240 

DIAGNOSTIC MHBll/77 

CALL TABL(HMAX,VMAX,HT(1).V0LI(1),1,1,1,I,IFBAD) 

WRITE(6,«)VTIL2,VMAX,V0REF,Q0,DT MHBll/77 

8 FORMATOt VTIL2 =»F10.3i<FT3, VMAX=XF10.3»FT3, VOREF=»F10.3XFT3,QO=J<MHBll/77 
1F10.3XFT3/SEC, DT=»F10.3>«SE(K) MHBll/77 

DELT=DT»<VTIL2-VMAX)/(V0REF-Q0*DT) MHBll/77 

URITE(6,10) 

10 FORMAT(20XX BASKET FULL*) 

CALL WRIT 

HI=HMAX 

T=T+DELT 

WRITE(6,110)T,HO,HI 

110 FORMAT!* TIME AT END OF TIME STEP = *F8.4*SEC, TANK LIQUID LEVEL 
1=*F8.4*FT, BASKET LIQUID LEVEL =*F8.4*FT*) 

60 TO JUMP, TRANSFER CONTROL TO MAIN PROGRAM-STATEMENT 10, IN ORDERMHBlO/77 
TO RESTART A NEW CASE ********x*k*k*x***********x***k****m****k**>>MHB10/77 
60 TO JUMP, (20) 

240 CALL WRIT 

IF((AB5(HMAX-HM)).6T.(HXXHT0L))60 TO 20 
WRITE(6,21)HT0L 

21 FORMAT(20X.XFILLING LEVEL WITHIN HS TIMES *F8.5*FTX) 

60 TO JUMP, (20) 

20 CONTINUE 
RETURN 
END 



SUBROUTINE WRIT 


Subroutine WRIT formats and directs the output for eaoh time step after time zero. 


SUBROUTINE WRIT 

C WRITES CONDITIONS AFTER EACH TIME STEP 

COmON/INPUT/DT,G,HT(20),HTOTrIfKC,NOUT,NST,OA,QP»RHOL,5IGMA. 

1VOLK20) 

CCr:M0N/REFIU/AL.AV»AVW.DPIN»DPS2,HItHILI»HIL2,H0.HS.HX, JUMP.QLt 
1Q0,T,VCL,VIL.VILMAX,V0,V0L0UT,V0REF,VV,WT 
CCmON/TIME/HTOL , TMAX 

C0r;M0H/WRIT/ADL,AVL,AVLS,AV5,AVT,AVWS,HDYN,Vl,VT,VVL,VVLS,VVS,VVM, 

ivvws 

REAL KC 
WRITE(6,11)QL 

11 FORMATCK INFLOW RATE DURING TIME STEP =^FS . ^J<CUBIC FEET/SECX) 
WRITE(6,20)VV,VVW,VVL,VL,VT 
WRITE(6,30)AV,AVW.AVL,AL,AVT 

20 FORMAT 0< VV = ^F8 . 5>tFT/SEC^7XMVVW=J^ FB . 5«FT/SECJ<5XXVVL=J<F8 . S^^FT/SEC^S 
lX^VL = )fF8.5i^FT/5EC«5X^ VT = kF8.5) 

30 FORMAT CX AV=xF8 . A»FT2>^5XKAVW=«F8 . AJ^FT2)<5Xi<AVL=JtF8 . A>fFT2X5X>tAL=>^F8 . 
1A>^FT25<5XJ( AVT = mF8.A) 

WRITE(6,31)VO»ADL,HDYN 

31 FORMATCX V0=xFS.5i< ADL=>^F8,5K HDYN=^F8.5) 

70 WRITE(6,71)VOREF.VCL,QO,VOLOUT 

71 FORMATO^ VOREF =5<F10.A,K VCL = XF10.A,K QO =i«F10.A,x VOLOUT =XF10. 
lA) 

WRITE(6.72) DPINrHX 

72 FORMAK^ DPIN = >^F10 .5»LBF/FT2 HX= J^F10.5»FT») 

WRITE(6,73) AVS.AVLS^AVWS 

73 FORMAT (K AVS = ^^F8 . A«FT2H5XhAVLS=^^F8 . AJ<FT2«5XXAVWS=xF8 . A><FT2X) 
WRITE(6,7A) VVS,VVLS»VVWS 

7 A FORMAT CK VVS = «FS . 5XFT/SECX5X«VVLS=XF8 . 5HFT/5ECX5X^VVWS = »F8 , 5«FPS^< ) 
WRITE(6.110) T,HO,HI 

110 FORMATO^ TIME AT END OF TIME STEP = KF8.AK5EC, TANK LIQUID LEVEL 
l=^ra.A^FT, BASKET LIQUID LEVEL =«F8 . AJ<FT?f ) 

IF(T.GE,TMAX)GO TO 200 
DO SO L=1,A 
50 WRITE(6,A0) 

AO FORMAT(IHO) 

WRITE(6,10) T.HO.HI 

10 FORMATCK INITIAL CONDITIONS, TIME =^^F8.AHSEC, TANK LIQUID LEVEL 
1 =^F8.AXFT, BASKET LIQUID LEVEL =XF8.AXFTX) 

200 RETURN 
END 



SUBROUTINE DIAGNOS 


Subroutine DIAGNOS is used for debugging and following the convergence pattern in the 
main subroutines. A printout of significant values in each Iteration is formatted and 
directed from this subroutine. (This subroutine is not normally required but is included 
for completeness of documentation. To delete the subroutine all CALL DIAGNOS 
references should be removed as well as SUBROUTINE DIAGNOS. ) 


SUBROUTINE DIAGNOS 

C DIAGNOS USED FOR PRINT OUT FOR DEBUGGING 
COMMON/DVCOL/AST,A2,B2, RANGE 

COnNON/FLOM/ALl,AL2,AVl,AV2,CONVERG,CONVl,CONV2,DPINl,DPIN2,HOLl, 
1H0L2, NITER 

CONMON/IMPINGE/A,AREA{20),ATOT,B,NDR,NIM.NM,VOLO(20).Y,Z 
COMMON/IMPOUT/CVOL , DPS.HDPIN.HPULL , VOL 

COMMON/INPUT/DT,G,HT(20),HTOT,I,KC,NOUT,NST,OA,QP,RHOL,SIGMA, 

1VOLK20) 

C0MM0N/REFILL/AL,AV,AVW,DPIN,DPS2,HI.HIL1,HIL2,H0,HS,HX.JUMP,()L, 
1Q0.T.VCL.VIL,VILMAX,V0,V0L0UT,V0REF,VV,WT 
COMriON/WRIT/APL,AVL,AVLS,AVS,AVT,AVMS,HDYN,VL,VT,VVL,VVLS,VVS,VVW, 
IVVU'S 
REAL KC 

WRITE(6,311)C0NVERG.DPIN,DPIN1,DPIN2,C0NVI,C0NV2 

311 FORMATOf CONVERG =kF 10.5« DPIN =*F10.5>f DPINl =KF10.5» DPIN2 =KF10 
1.5»e CONVl =xF10.5» C0NV2 =XF10.5) 

WRITE(6,312)HDPIN,HX.H0L1,H0L2.HIL1,HIL2 

312 FORMATOf HDPIN:>«F1 0 . HX=«F10 . 5><H0L1 =XF10 . 5x H0L2=*F10 . 5X HILl=xF 
110. 5*« HIL2:«F10.5) 

WRITE(6,212)H0,HI 

212 FORMAHx H0 = XF10.5»< HI=«F10.5) 

WRITE(6,513)NITER.VV,VVW,VVL,VL,VT 

313 FORMAT(J< ITERATION NUMBER J«I3« VV=XF10.5XVVW=»»F10.5* VVL=XF10.5x 
1 VL=xF10.5,5XXVT=xF8.5) 

WRITE(6,30)AV,AVU, AVL,AL.AVT 

30 FORMAT (X AV=XF8.5>FT2X5XXAVW=XF8.5XFT2><5XXAVL=XF8.5XFT2X5XXAL=XF8. 
15*(FT2x5XxaVT = <F8.5) 

WRITE(6,31'>)AV1,AV2,AL1,AL2,0L,Q0 

31A FORMATTX AV1= XF10.5X AV2= XFIO.SX AL1= XFIO . 5XAL2=XF10 . 5x QL= »F1 
10. 5x Q0 =xF 1C.5) 
lv'RITE(6,73) AVS,AVLS,AVWS 

73 FORMAT (X AVS=XF8 . 5xFT2x5XxAVLS=xF 8 . 5xFT2x5XXAVUS=XF8 . 5XFT2X ) 
K'.RITE(6,74) VVS,VVLS,VV:.'S 

7A FORMAT (X VVS=xF8 . 5XFT/SECX5XXVVLS=XF8 . 5XFT/SECX5XX VVWS=xF8 . 5XFPSX) 
WRITE(6.315)V0L,VIL,V0REF 

315 F 0 RMAT(x VOL =xF10.5x VIL =xF10.5x V0REF= XF10.5) 

WR1TE(6,316) 

316 FORMATCIH ) 

RETURN 

END 
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A. 3 PROGRAM REFILL SAMPLE INPUT 


The program input Is In namelist format. Definitions of input variables are given in 
the initial portion of subroutine INPUT. The input shown on the following pages is 
for the Centaur D-IS LO 2 tank refilling during the fourth burn of the five bum 
low earth orbit mission. 


PROGRAM INPUT 


9INDATA 

A*0.I63, 

AREA«0. 0,2.016,6. 581, 8.088, 9. 023,9.397,11.0^1, 14. 145,24.601, 24. 444, 
24.445, 

AST»0.043, 

AT0T=24.444, 

AUIDTH=30.97, 

A2*0.I63, 

8 = 0 . 11 , 

82=0.11, 

CC=1.0, 

0*0.0015, 

«BP=38.0. 

DBP2=65.0. 

noii;=o.208, 

DPRESS=5000.0, 

DPX=0.25, 

DT=0.025, 

EH=557.46, 

G=69.?4, 

HC0L=2-295, 2.3,2.31, 

HPT =0.0,0. 032, 

HT=0. 0,0. 02083,0. 1 083, 0.144, 0.226,0.258,0. 403, 0.759,0. 959, 1.1425, 
2.3042, 

HTOL=0.05, 

HT0T=1.1425, 

T=11, 

J=2, 

K=3, 

KC=0.6. 

MUL=0.000108, 

N=3, 

NG=1, 

N0UT=5, 

NS=2, 

NST=1 , 

NTSCR=1, 

N»=2, 

0A=0.51, 

P=0.635, 

PT»0. 00208, 

00=0.835, 
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80UT«0.0»».83S, 

RAIIGE*0.0004, 

RH0L-6J.5, 

SI BHA»0. 000782, 

ST=0. 00103, 

TC«0.0,0.01.50.0, 

T»ETA=I.571, 

nTLEsJlHCEHTABR 0-1 S L02 BASKET 50X250, 

THAX»17.0, 

VILI*0.t35, 

MOLI=0.0, 0.0325,0.701 ,0.732, 1.45, 1.79, 3.31, 0.91 ,7.57,7.5701,7.5702, 
m0>0. 3073, 0.304, 0.4480, 0.7908, 1.1006, 1.2406, 2. 3384, 8.341 1,14. 8523, 
21.0389,94.942, 

VT0T»0. 3073, 0.3965, 1.1496, 1.5288,2.5500, 3.0306,5.6484, 15.2511, 22.4223, 
29.209,102.5126, 

Y»3.38, 

Z»15.01, 

SEND 


A. 4 PROGRAM REFILL SAMPLE OUTPUT 


T he output shown on the following pages corresponds to the case input shown in A. 3. 
Results for the 4th bum of the five-bum low earth orbit mission are shown for 
Centaur D-IS LO 2 start basket refilling. 


Initial conditions are printed first. The first two time steps and the last time step 
are shown. For the first time step convergence was obtained in less than ten 
iterations. The second time step shows the printout driven from subroi^ineDIAGNOS 
if greater than ten iterations are required. This feature was used extensively in the 
debugging state of program development in order to eliminate ccmvergence problems 
and programming errors. (Subroutine DIAGNOS and associated subroutine call 
statements can be eliminated from the program to reduce computer and printer costs. 
Calculations performed will be identical with or without subroutine DIAGNOS. ) 


The calculations show that the capillary device refills into the standpipe region in 
2. 5521 seconds. This is well below the allowable refilling time of 16. 9 seconds for 
this bum. 
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